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ABsnucr 

STOCBASIIC  ANALYSE  CF  FAOUIIES  HARIKT^ 

C}QNVE^mOP^  WEAPCIE 

A  procedure  is  developed  ixtimis  a  probdnlity  of  kill  of  a  hardened  &cility 

taking  into  account  two  types  of  uiKxrlffiitt^  weapon  deliveo^  aocui»^  and  stmctural 
diaracteristics  or  intelligence  uncertainties,  llie  kill  criteria  are  based  on  the  structural 
response  of  die  fiadlity  exceeding  predetermined  limits  which  rqxesent  the  achievement 
of  die  attack  objective.  Perfect  knowledge  is  rardy  known  about  die  structural 
characterisdcsofa  target  once  a  conflict  is  initiated.  Analyst  tasked  to  perform  pre-attack 
wesqxms  analysis  and  post-attack  weqxxis  effectiveness  must  be  dile  to  rqxxt  to  thdr 
superiors  realistic  probabilities  of  achieving  die  objective  of  an  airbome  strike  on  atarget 
Current  mediodologies  do  incorporate  weapon  delivery  accuracy,  however  they  overlodc 
uncenaindes  in  target  structural  diaracterisdcs  which  can  make  a  dramatic  difitoice  in 
a  probability  of  kill  predicdon. 

A  mxilinear,  ncmdimensional,  cxxitinuous  l^steretic  beam  modd  is  devdoped  to 
represent  a  section  of  a  hardened  fecili^  subjected  to  caivendonal  weqx)ns  effects.  The 
model  returns  reqxxise  calculations  across  die  hdght  of  the  section  as  required  to  provide 
informaticm  for  determining  the  kid  state  of  a  hardened  target  The 

nondimensicxialization  allows  fen*  ease  of  parameter  iiput  and  serves  the  stodiasdc 
analysis  well  where  structural  diaracteristics  are  continuously  changed 

Existing  enseal  models  which  generate  conventional  weapon  blast  pressure  time 
histexies  as  a  function  of  the  TNT  equivalent  throw  wdght  are  modified  to  become  a 
function  of  ^lace  as  well  as  time.  A  new  model  is  generated  which  returns  the  pressure 
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at  aiQr  point  up  a  ^waU  section  as  a  iinction  of  time;  i9»c«  and  angle  of  hidden  This 
type  of  load  representation  is  required  to  feed  tiie  continuous  beam  model  referenced 
^xive.  The  cond>ination  of  tiie  beam  model  and  tiie  load  model  is  termed  the  NONUN 
model. 

Robust  statistical  models  or  re^xnse  sur&ces  (RS)  are  derived  from  NONUN 
results  calculated  using  typical  combinations  of  weqxm  throw  weight  and  range,  taget 
wall  hei^its,  dqjtfis  and  concrete  compressive  strenglfas.  Tlie  use  of  a  dedgn  of 
e^qseriments  (DCB)  or  experimental  design  approach  to  the  RS  development  ensured  the 
RS  closely  replicated  the  iipm  data  across  the  parameter  space  of  interest  Theresultwas 
a  multidimensi(xial  RS  returned  tiie  structural  reqxnse  given  a  set  of  the  five 
parameters  stated  above. 

Two  Monte  Carlo  simulation  programs  are  developed  which  take  into  account 
structural  diaractoistics  as  randexn  inpiti  variables  in  addition  to  tiie  traditional  wetpon 
delivery  accuracy.  The  first  program  utilizes  the  NCX^IUN  model  in  the  siinulaticms 
whoeas  the  seemd  program  rqilaced  tiie  full  NONUN  model  witii  several  RS  models 
valid  ovCT  key  parameter  ranges.  Use  of  the  RS  rqilaoenient  models  allowed  tiie 
simulaticxi  to  be  run  in  less  than  1%  of  the  time  required  to  run  tiie  simulation  with 
NONUN.  This  time  saving  is  crucial  to  tiie  use  of  this  tool  in  a  d|ynamic  wartime 
enviremment 

This  work  demonstrates  that  uncertainties  in  the  structural  characteristics  of  a 
target  may  significantly  effect  its  reqxxise  to  conventional  weqxxis  and  the  determination 
of  the  resulting  probability  of  kill  given  an  attadL  The  use  of  robust  re^xxise  surfiices 


to  replace  ooini^  analytic  procedures  demonstrates  diat  trnidy  calculation  of 
probabilities  of  damage  may  be  generated  in  ^ite  of  using  a  simulation  technique  such 
as  Monte  Carlo.  The  medxxiology  presented  will  accommodate  studies  to  single  out 
die  most  critical  random  structural  variables  Old  dieir  ranges  allowing  the  proper  empha^ 
to  be  placed  on  variable  significance  in  a  targeting  analy^  and  data  gathering  exercise. 
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Uroduction 

In  the  past  twenty-five  years  diere  has  been  a  resurgence  of  intere^  in  the  effects 
of  c(xiventi(ml  wesqxns  (xi  structures  (Kraudiammer,1987).  This  is  due  in  part  to  die 
many  uncertainties  vdiich  diaracterize  the  loading  and  re^xmse  process  of  a  structure  in 
a  convendoial  weqxxi  environment  plus  die  realization  that  diere  is  singly  no  Mire  in 
the  use  of  nuclear  werpxis.  The  occurrence  of  EXESERT  STY3RM  dnlher  served  to 
illustrate  the  effecdvoiess  of  advanced  axivendonal  wefqxxis  against  targets  which  only 
two  decades  ago  could  (xily  be  held  at  risk  widi  nuclear  weqxxis.  The  introduction  of 
precisi(xi  guided  muniticxis  and  tough,  penetrating  warheads  has  caused  a  re-diinking  of 
the  "near  miss"  mentality  diat  has  dominated  "survivability  designs"  to  date. 

OxicurTently,  the  mediods  by  which  "hostile"  structures  are  targeted  and 
we^neered  are  also  being  redefined.  In  targeting,  uncertainties  in  delivery  accuracy  are 
accoiqianied  by  uncertainties  in  the  structural  characteristics  of  the  target,  such  as  wall 
thickness  (internal  and  external),  pocent  of  steel  reinforcemoit,  concrete  strengdi,  etc.. 
The  target  structural  charactoistics  uncoUunties  above  will  be  referred  to  as  target 
intelligence  uncertainties  for  the  remaindo’  of  diis  dissertation.  The  ability  to  perform 
meaningful  pre-  and  post-attack  Battle  Damage  Assessment  (BDA)  is  driven  ty  the 
amount  of  information  known  about  the  target  coi^led  with  the  analyst's  c^)abiliy  to 
perform  weapcxis  effects  calculations.  Qnrent  probabilistic  structural  re^xmse  methods 
do  not  inccxpcxate  target  structural  uncotainties  into  their  probability  of  kill  calculations. 

The  primary  ccxicem  of  commanders  during  a  conflict  is  whether  diey 


2 

accon^lisheddieiriiiilhaiyolyecdvewifhtfiefin^  For  air  ooimiianders  this 

translates  into  answering  tilie  following  question;  Did  an  attadc  achieve  the  level  of 
damage  desired  or  do  more  pilots  need  to  be  sent  in  harm's  w^  to  re-strike  the  target? 
Coihbat  decisirxi  makers  require  personnel  to  perfonn  pre-  and  post-attack  damage 
analysis  and  provide  realistic  assessment  given  whatever  level  of  target  intelligence  is 
available,  fo  a  oxnbat  environment  or  in  prqjaradon  for  a  contingent^  operation,  damage 
predictions  are  required  in  nev  real  time.  Apu^isontoincxaporatestnicturaldaniage 
predictions  and  graphical  rqxesentations  into  cmfoet  simulators  and  trainers  foerefore 
exacefoating  foe  real  time  requirement  Realistic  bounds  are  required  for  foe  stnjctural 
diaracteristics  presented  in  Table  42.  The  mean  values  and  probability  distribution  of 
eadi  variable  will  be  based  on  known  local  design  and  construction  practice.  For 
example,  it  was  found  during  on-site  inflections  following  DESERT  STORM  that 
similarly  designed  airoait  shelters  built  in  Iraq  and  Kuwait  had.veiy  different  quality 
control  on  ccxicrete  mix  and  placement  This  was  evident  in  the  refxmse  of  foidc 
concrete  burster  liters  to  contact  detcxiations  and  penetrating  wefxms.  There  is  mudi 
work  underwE^  within  the  DoD  communi^  to  quantify  foe  sensitivify  to  variations  in  the 
design  characteristics  in  Table  42  when  sutgected  to  blast  fragment  loadings  as  well  as 
projectile  penetrations.  A  procedure  to  qiply  a  measure  of  effectiveness  or  probabilify 
of  damage  to  tiiis  process  is  lacking  in  the  communify  at  this  time.  This  work  will 
present  a  mefoodology  to  ocnnbine  foe  weqxms  effects  and  target  intelligence 
uncatainties.  This  mefoodology  will  be  illustrated  for  a  single  case,  that  of  an  above 
grcxind  hardened  structure  targeted  wifo  non-precision  wefxxis,  however  fiplicatioi  to 
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adier  scenarios  ate  rdevant  and  will  be  discussed. 

A  conventional  weapon  blast  produces  vay  high  local  pressures  of  ve^^  short 
duration.  This  is  in  contrast  to  the  much  hitter  and  much  longer  pressure  loading  of  a 
nuclear  blast  In  die  converidonal  arena,  structures  are  designed  to  widistand  a  ^)ecified 
weqxxi  damage  mechanism  based  on  a  threat  assessment  for  that  &dlity.  Forexamfde, 
structures  >^ch  house  hi^  value  assets  or  critical  mission  hmctions  will  be  derigned  to 
withstand  the  rnost  potent  weqxxi  dseat  Ihe  rnost  dseatenmg  conventional  weqxns 
existing  tod^  are  precision  guided,  penetrating  weqxxis.  Structures  designed  to  resi^ 
such  weqxms  will  not  be  studied  her^  Structures  i^ch  house  less  than  mission 
critical  functions  are  traditionally  constructed  above  ground  and  designed  to  withstand  the 
primary  damage  mechanisms  of  unguided  munitions,  i^ch  are  airbL  .  and  fiagments. 
Unguided  munitians  may  be  delivered  by  tactical  fighter-bomber  aircraft  (ie.  F-1 1 1,  F- 
ISE,  F-1 17)  CM*  strategic  bombers  (ie.  B52,  Bl).  The  unguided  munitions  of  interest  are 
delivered  in  groiqis  of  up  to  ei^  in  a  "stidc"  for  tactical  aircraft  and  much  larger 
numbers  for  strategic  aircraft.  The  ccnnbination  of  tiirow  weight  of  eiqilosive  and  location 
of  detonation  determine  the  magnitude  and  duration  of  a  conventional  loading  i^di  may 
produce  pressures  on  die  order  of  30,000  psi  decaying  esqionentially  to  atmospheric 
pressure  oyer  a  poiod  of  1  to  4  milliseconds(msec).  From  a  survivability  stand^int, 
structures  designed  to  withstand  standoff  loadings  of  this  type  are  typically  called 
"semihardened”  and  are  traditionally  constructed  of  reinforced  concrete  fixxn  12  to  30 
indies  in  thickness.  In  contrast,  structures  designed  to  witiistand  direct  hits  fixxn 
conventiaial  wetqxms  are  built  oon^letely  or  partially  underground,  widi  tfaidcer  walls 


Attack  scenario 

Duration  of  attack 
Target  priorities 
Sin^e  target 

Blanket  bon^l»ng^inuhq>ie  targets 
Delivery  method 
Airborne 
Sur&oe 
Naval 

Nkilt^le  attadcs 
Weiqxntype 

Fh^cal  dimensions 

LengdiAviddi  ratio 
Casing  diidcness 
'Ihrow  wei^  (TNT  equivalent) 

Ouidanoe  system 
Detonation  method 
In^iact 
Timed 
Penetration 
Delayed  fuse 
Pressure-Time  histoy 

Pressure  and  casing  fragment  synergism 
Spadal  representation 
Arrival  aid  duration  representation 
Initial  condition  rqjresentatioi 


Table  1.1  Loading  Uncertainties 
and  a  burster  barrier  cover. 

To  fully  define  an  attack,  and  die  subsequent  damage  mechanism  seen  fay  a 
structure,  one  must  start  >vidi  die  intent  of  die  attack.  Once  a  plan  of  is 
hypothesized,  the  probability  of  survival/damage  of  a  structure  can  be  developed  by  using 
ccxipita*  models  to  predict  die  reqxxise.  I^iedicting  die  re^xxise  of  a  known  structure 
to  die  direat  werqxms  of  interest  is  a  random  process  due  to  the  uncertainties  in  the 
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structural  material  properties  as  wdl  as  inherent  uncertainties  in  the  anulatkn 
mediodologyr.  Ihe  random  uncertainties  assodated  with  the  attadc/kxKling  aid 
structiv^c^pcnse  events,  given  known  structural  chaacleristics,  ae  shown  in  1.1 
and  \2. 

The  mryofity  of  woric  dealing  with  the  reqxnse  of  individual  structures  to 
conventional  wer^xxis  oonaders  deterministic  and  often  very  simplistic  moddk 
^)ecifically,  spatial  variations  in  tiie  loading  are  ignored  and  the  structure  b  represented 


Ihysical  properties 

Material  strength 

Nonlinear  te^xmse 

Strength  enhancement  due  to  high  stress,  idiort  duration  loading 
Compoate  action 

Blast  arresting  device  physical  characteristics 
Berms 


Revetments 
Sacrificial  slabs 

Multiple  hits  and  strengdi  denudation 
Spalling^scabbing  on  fiont  wall 
Blast  attenuaticxi  and  ^jailing  on  inside  wall 
Local  breadung 
Wolonanship 
Sq^iistication  of  m^l 

Single  degree  of  fiuedc^SDOF) 

Miltidegree  of  fieedc»n(MIX>F) 

Ccxitinuous 
Finite  difference 
Finite  element 
Sqjport  conditions 
Mode  of  Mure 
Shear 


Flexure 


Table  12  Structure  Uncertainties 
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by  a  angle  d^ree  of  freedom  modd.  Ihe  weapons  effects  otgective  of  diis  dissert 
is  to  provide  an  analysis  technique  wiiich  takes  into  account  die  stodiasdc  loading  aid 
resistance  characteristics  of  die  system  while  using  one  dimensional  continuous  load  and 
structure  representations.  The  target  intelligenoe  uncertainties  will  be  incorporated  into 
die  structure  rqxesentation  as  well. 

Oiqiter  2  reviews  the  current  deterministic  and  stochastic  methods  of  design  aid 
andyasofhardened  as  wdl  as  traditional  ftcUides.  Research  on-gcang  in  diese  areas  will 
be  examined,  dit^iter  3  describes  die  development  of  the  re^xxise  statistic  generator 
code  to  include  die  load  and  structural  modds  and  the  solution  tednique.  Chester  4 
describes  die  basis  and  in^lementation  of  die  stochastic  response  anafysis.  Stochastic 
analyses  are  coi^iared  for  various  real-wnrld  delivery  situations  and  stnictural  variable 
distributions.  Qa^iter  5  provides  reconmiendatians  and  conchiaons. 
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CHAPIfRl 

BackgiTNiid  nd  Lilentim 

2.1  DelmiriirisIkDesi^iaDdAiialj^oflVolBCtiveStnKtiicsAgriiBt^^^^ 

Hie  mediodology  used  tfie  Armed  Forces  in  desigp  and  analysis  of  semi- 
hodened  stnictures  is  ouUined  in  die  US  Air  Fcice  maiBials, ’Trotectkn  Fitm  Fkxni^ 
Weqx)ns"(Chnvford,1971),  "Stnictures  to  Resist  die  Effects  of  Accidental 
Eiqdosioiis"(TMS-1300,1990),  Tratecdve  Consliuction  Deagp  MBnual''(  Drake  et 
al.,1989)  and  die  Army  manual  "Fundamentals  of  Protective  Design  For  Conventional 
Weapons'’CIM5-855,1986). 

Drake  etal,1989  is  the  current  audioriQr  for  the  US.  Air  Force.  According  to  diis 
manual,  die  diice  cases  of  interest  for  sur&cestiiictures  being  considered  are:  (l)direct 
hits  resulting  in  extreme  damage,  (2)  close-in  bursts  causing  extremdly  high  and 
nonuniform  overpressures  and  fiagment  loadings,  and  (3)  bursts,  occurring  fer  erou^ 
away  for  the  pressure  loadings  to  be  r^iprmdmately  uuform.  With  the  accurate  of 
modanwa9X)ns,(l)  and  (2)  are  becoming  die  events  of  interest  They  are  also  the  areas 
about  which  the  least  is  known. 

Ihe  stqs  in  the  analysis  and  design  of  hardened  structures  (Drdce  et  al,1989)  are 
given  in  Table  2.1.  In  the  ronainder  of  diis  section  die  current  methodologies  used  to 
perform  those  stqis  will  be  discussed 

2.1.1  Loadhig  CbnrideralioiiB 

The  eiT^cal  forinulas  refered  to  in  Stqil  above,  were  developed  to  determine 
the  minimum  wall  thickness  required  to  prevent  breaching  fiom  contact  or  near  contact 


Stq>  1  -  As  a  pcdUminafy  aialy^  use  bieadiing  curves  (described  bdow) 
to  esdinate  the  minimum  standoff  retirements  of  the  wall  section  to  resist 
breaching. 

Stq)  2  *  Perform  a  structural  analysis  assuming  the  wet^xxi  burst  produces 
a  uniform  load,  provided  die  total  ^jfdied  impulse  to  die  stnictiBie  is  maintained. 

The  equivalent  uniform  load  is  in  a  sin^e  degree  of  fieedom  (SDCXO 
analysis  resulting  in  the  maximum  structural  re^xmse  at  a  ^lecified  loodon. 

Step  3  -  Perform  a  fiagmentation  analysis  to  determine  if  the  design 
fragment,  based  on  die  desi^i  or  threat  weapon  case  characteristics,  will 
perforate  die  wall. 

Table  2.1  Stqs  in  the  Desigp  and  Analysis  of  Hardened  Struchses 
bursts.  These  formulas  are  based  on  National  Defense  Research  Comniittee(NIX^  tests 
(White,1946),  sipplemented  by  more  recent  test  data  (McVay,1988).  These  curves  can 
be  found  in  recem  manuals  (Drake  etal,1989*JDAHS,I993).  Examples  ofthese  curves  are 
shown  m  Figures  2.1.fflid  11  and  were  taken  out  of  Hyde,1988.  These  curves  are  for 
cased  and  uncased  oqiiosives  and  are  based,  for  the  most  part,  on  scaled  model  tests.  It 
has  been  reported  that  strain  rate  effects  (McVay,1988),  cannot  be  modeled  or  scaled 
ejqilicitly,  causing  model  structures  to  widistand  slightly  less  damage  dian  otherwise 
equivalent  full-scale  structures.  This  foct  may  result  in  the  under  prediction  of  real  world 
damage  based  (xi  these  curves.  Breach  and  lower  levels  of  damage  m^  be  read  fixxn 
these  curves  using  the  scaled  wall  thickness,  and  the  scaled  range  For 

cased  weapons  the  scaled  range  is  multipled  by  where  C  is  the  case  weight 


Figure  2.1  NI^C  Curves  for  a  Cased  Ejqplosive  (Hyde,1988) 


Figure  22  NEttlC  Curves  for  an  Iftioeed  E7q)losive  (ffyd^l988) 
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IheequBtiGn 


t 

ffl/3 


-b 


w*c 


(2.1) 


was  used  to  generate  the  curves  of  Figises  2.1and  22.  Tfd>le  22  gives  die  values  of 
a,  b,  and  c  used  to  generate  those  curves. 

Descriptions  of  near  field  bursts  have  traditionally  dtst^^rded  die  nonunifonnity 
of  die  blast  loading  as  well  as  weqxxi  caang  effects  (Q:awfocd,1971;  TN&1300,1990). 
In  fact,  it  has  been  shown  diat  the  blast  loading  is  very  concentrated  about  the  centrcHd 
of  the  weqxxL  Ihe  effects  of  the  wer^pon  case  on  the  load  attributed  to  a  werq^  has 
been  widely  ddioted.  One  method  of  including  such  effects  is  to  decrease  die  blast 
loading  to  account  for  die  energy  absorbed  in  case  break-up.  It  is  correct  that  the  blast 
pressure  seen  by  the  structure  fixnn  a  given  ^  cased  eiqilosive  is  less  than  if  the  charge 


a 

b 

c 

No  Damage 
Cased  Weqxm 

0.43 

03 

0.3 

1  Beach 

1  Cased  We^xxi 

0.23 

0.3 

0.3 

1  No  Damage 
Bare  Charge 

0.3 

0.62 

0 

Breadi 

Bare  Charge 

0.13 

0.62 

0 

Table  22  NDRC  Ec^iation  Coefficents  (Ifyde,1988) 
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wereuncased  In^inelhodhonwever,theini|MlseddiveraibyAefiagn^ 

Ibis  impulse  has  been  shown  in  sane  cases  to  be  very  destmcdve  (I^der,1983).  Ihe 
current  method  (Drake  et  al,19S9d)AHS,1993)  does  not  reduce  the  blast  loading  due  to 
case  bfeak-^>  but  does  not  esqdicitly  aocouit  for  die  fiagment  loading.  Only  recently 
halve  researchers  becone  very  interested  in  the  condmied  or  synei^stic  effects  of  blast 
and  firagments  fixxn  cased  esgilosives  (Hader,1987;  Marchandl986;  Kropatsche(h,1983; 
Koos,1987).  Studies  to  eiqdictly  model  fiagpaent  loading  and  synergistic  effects  have 
been  ^xxisoed  1^  various  DoD  agencies  beginning  in  1988  (Marchand,1988;  Sues  aid 
Twisdale,1993). 

2.U  Stnictural  Response  Cbnsideiations 

The  SDOF  commonfy  selected  is  die  mici^int  deflection  or  support  rotation  of 
a  slab  representing  the  portion  of  the  protective  structure.  Ihe  wall  dimensions  typical 
in  protective  ccxistruction  in  coijunctiai  widi  observed  re^xxisejcharacteristics  has  led 
to  the  commcxi  modeling  of  wall  secdcns  as  cxie-way  slabs  (Biggs,1964;  IM-5-8SS- 
1,1986).  Equivalent  SDOF  parameter  rqxesentaticxis  of  such  structures  are  provided  in 
tabular  form  in  several  publications  (Bigg^l964;  Crawford  et  al,1971;  IM-5-855-1,1986). 

Ihe  Justificaticm  fen*  use  of  the  uniform  load  and  SDOF  models  hinges  on  die 
uncertainties  surrounding  the  loading  and  material  re^xxises  not  warranting  elaborate 
and/or  expensive  analysis  techniques  {Biggs,1964).  Additionally,  in  protective 
coistructioi  design,  the  main  design  driver  has  been  maximum  displacement,  which  has 
beoi  shown  to  be  modeled  reasonably  well  widi  the  SDOF  rqxesentation  for  above 
ground  and  shallow  buried  structures  (Cohhaip  et  al,198S;  Kraudiammer,1986).  Ihe 
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short&lls  of  the  SDC^  metfiod  diis  dissertation  addresses  are: 

1.  Re^xxise  diaracteristics  generated  for  the  designated  degree  of  fieedmn  only. 

2.  ^)atial  load  variations,  to  include  fiagment  effects,  can  not  be  analyzed. 

3.  9iear  fixmi  the  sqsports  cfn  not  be  analyzed  directly. 

4.  Si^iificaiit  participation  ofhi^ier  inodes  ofvibration,ty|ncal  under  inpilsive 
loadings,  canned  be  evaluated. 

A^/hen  imder  ncffl’ field  loading,  tiiear  fidhire  has  been  tiiown  to  be  prevalent  (Ross 
and  Krawinkler,198S;  Krauthammer, 198Q  as  well  as  die  early  achievement  of  the  ultimate 
mcHnent  ciqiacity  at  points  away  from  the  midspan.  the  addition  of  fiagment 
inq^acts  and  maximum  nonuniform  pressures  near  the  bottrxn  of  die  wall,  shear  Mures 
away  firom  the  siqjports  have  been  observed.  This  so  called  punching  shear  or  vertical 
shear  fiiilure  has  been  called  an  early  time  Mure  (Ross  and  Kiawiiilda,1985;  Van  Der 
Veen  and  Blaauwendiaaad,1983;  KrauthanBna*,lS^.  The  Mure  occurs  so  early  diat 
there  qspears  to  be  no  flexural  contribution.  Thus  it  has  been  proposed  that  sqiarate  shear 
and  flexural  criteria  can  be  used  to  predict  Mure  once  the  reqxinse  has  been  ompited 
(Ross  and  Rosengren,1985;  Krauthammer,1986).  If  the  structure  does  not  develrp  a 
critical  shear  re^xmse,  it  can  fail  lata-  in  flexure.  The  ultimate  mcxnent  may  occur  away 
fiom  the  midspan  if  the  pessures  are  sufficioitly  hi^  and  are  coicentrated  in  cm  area, 
normally  the  bottom  half  of  the  wall.  In  these  cases  yielding  dien  propagates  toward  the 
mid^ian  (Ross  and  Rosaigrai,1985;  Krauthammer,1986). 

12  Probabilistic  Metiiods  in  Design  and  Analysis  of  IVttiBctive  Stnirtures  Against 
A^ie^ions  Effects 
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111  Gfeoerai 

Probabilistic  analysis  of  protective  structures  was  first  used  for  hardened  &cilities, 
or  facilities  required  to  resist  nuclear  loacb.  In  1968  a  program  called  FAST  in  was 
completed  whidi  modeled  various  elements  of  die  hardened  or  strategic  ^'stems  Mlure 
problem  (R0wan,1977).  Failure  Analysis  by  Statistical  Techniques  (FAST)  was  developed 
to  sijtfy  nuclear  attack  survivability,  using  Mmte  Carlo  simulaticm.  >^fithin  ^  program 
simulations  are  run  with  probable  attadc  scenarios  and  possible  ftihire  modes  and 
probability  statistics  are  conpjted.  The  code  underwent  several  revisions  (Rowan, 1977) 
until  it  was  iqilaoed  Ity  Probability  Assesanent  of  Strategic  Systenos  (PASS)  code  adopted 
by  the  Air  Fcxce  in  1988  (Kung  et  al,1988).  The  new  code,  which  covers  a  wider  variety 
of  scenarios  and  failure  medianisms  than  FAST,  mirras  current  technologies  but  uses  the 
same  Mmte  Carlo  approach  to  derive  re^xxise  and  survivahility  statistics. 

Several  air  base  attack  simulation  models  use  apr(habilistic.£pproQch  for  assessing 
overall  damage  effects  fiom  a  (xodetermined  attack  scenario.  One  sudi  program  is 
entitled  TSARINA,  and  was  develq3ed  by  the  RAND  Crapraafioa  The  main  objective 
of  this  code  is  to  goierate  air  base  attack  outccanes  that  exhibit  realistic  levels  of  damage 
across  Mcmte  Carlo  trials.  It  was  not  designed  to  accurately  estimate  die  damage  to  any 
particular  target  for  a  particular  impact  point  and  wespon.  Damage  assessmem  is 
acconplished  using  a  "cookie-cutto^*  method  based  on  we^xm  effectiveiKSs  data  takoi 
fixim  the  Joint  Munitions  Effectiveness  N^uals/  Air  to  Surface  (JMEM^AS)  manual 
whidi  will  be  discussed  later  in  the  di^iter.  Ftx*  eadi  kind  of  point  inpact  weprai  the 
code  assigns  an  effective  miss  distance  (EMD)  accrading  to  each  target  type. 
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Semihardened  Qpe  &cilities  would  be  one  sudi  target  type.  If  a  fecility  target  fiiUs 
within  die  radius  of  a  weqxxis  EMD  after  an  attadc,  die  fiaction  of  die  facility  coveted 
leads  to  die  determinadoi  of  the  fiacd(xial  damage  to  die  fecility.  Cumulative  damage 
fitMn  multiple  point  in^iocts  m^  be  attributed  to  a  &cility.  Failure  criteria  for  different 
types  of  ftcUities  are  based  on  reaching  a  certain  fractional  damage  level. 

The  EfifecdvenessA^ulnerability  Assessment  in  Three  E>imensions  (EVA3D),  is  a 
model  developed  under  funding  from  the  Air  Force  to  assess  the  vulnerability  of  hard 
targets  (Bessette,1988).  The  results  of  diis  code  could  feed  an  overall  air  base  or  theater 
ccHiflict  simuladon  model  such  as  TSAR  mentioned  thove.  EVA3D  is  a  Monte  Carlo 
simulation  model  \^ch  sedcs  to  provide  a  realistic  assessment  of  damage  to  a  structure 
and  its  ccmtents.  Either  a  structural  kill  or  a  funcdcxial  kill  may  occur.  The  ability  to 
predict  a  fiincticxial  kill  is  much  more  dqpendent  on  the  fordknowledge  or  intelligence 
available  on  that  facility.  Within  die  code  a  variety  of  predictive  "tools"  are  available. 
The  user  has  die  ability  to  £^ly  a  fcnee  in  one  of  four  ways.  They  are:  (1)  laser- 
guided  bomb  delivery,  (2)  electro-qitical  guided  bomb  delivery,  (3)  unguided  ch*  sdcdc 
bomb  delivery  and  (4)  Pk  regicxi  anai}^is.  Cases  1  and  2  both  deal  widi  precisicai  guided 
muniticms  and  are  not  relevant  to  this  work.  Case  4  is  available  to  qiecifically  provide 
output  data  that  would  be  consistent  with  the  data  recpiirements  of  the  TSARASARINA 
model.  The  stick  bomb  delivoy  method  is  the  cme  of  choice  for  this  effort  and  will  be 
described  in  the  context  of  its  EVA3D  ^licadcm  as  background  EVA3D  uses  a 
MSBASIC  program  develc^ied  by  the  JMEM  community  fcff  evaluating  runway 
survivability  from  an  attack  by  an  aircraft  carrying  a  stick:  of  "dumb"  bcnnbs.  A  stick 
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delivay  is  a  mediod  of  bombing  in  ^cfa  two  or  more  bombs  are  rdeased  at  a 
predetermined  interval  frcxn  one  aircraft  as  a  result  of  a  single  actuation  of  the  bondv 
release  mechanism  on  die  aircraft.  Stidc  patterns  or  how  the  weapons  are  delivered 
relative  to  a  reference  point  are  generated  imng  the  JMEM^AS  data  diat  is  fed  into  the 
Sdddbomb  code.  The  stidc  patterns  are  a  function  o^  ddivay  profile,  die  weapons 
loadout,  intervelcxneter  setting,  and  weqxxi  aerodynamic/stability  characteristics.  These 
variables  are  iipits  into  the  Stiddxxnb  code.  After  the  Stickbomb  code  returns  a  stick 
pattern  relative  to  a  reference  point  the  EVA3D  code  determines  the  actual  impact  points 
by  implying  a  prooedLire  diat  accounts  for  aiming  error  and  ballistic  error.  Each  bond) 
also  is  evaluated  for  fuse  reliability,  that  is  one  bcxid)  in  so  many  is  conddered  a  "dud.” 
Fcm*  the  purposes  of  diis  work  these  "fitters"  will  not  be  qiplied  and  die  reference  point 
dted  above  will  beccane  the  aim  point  or  center  of  the  target  A  samfde  stick  pattern  is 
shown  in  Figure  2.3. 

Giou  and  Chang,1987  arrests  using  Markov  chains  to  model  die  detonation 
locations  of  wespons  during  a  sequoKe  of  attacks.  This  way  the  actual  load  seal  by  a 
structure  may  be  aiproximated  since  the  loading  is  directly  proportional  to  die  structures 
distance  fium  the  detcxiaticxi  point  of  the  weqxm. 

Wong,1985,  extended  the  fir^-orda,  seccnid-mcHnent  mediod  in  attonpting  to 
quantify  the  uncotainties  in  the  represoitaticxi  of  dynamic  soil-structure  intaaction(SSI). 
Dynamic  SSI  refers  to  the  coipling  between  the  reqxxise  of  a  buried  structure  and  the 
loading  exoted  by  the  medium  on  the  structure  afta  the  dehxiadon  of  a  conventional 
we^)cxi  eitha  above  or  below  ground.  The  parameters  modeled  as  uncertain  are  the 


16 


Figure  2.3  l^icil  Sdck  Pattern  About  An  Aim  Point 

nominal  overpressure  loading,  die  loading  moduli,  die  unloading  moduli,  die  pressure  rise 
time  and  the  soil  density. 

Researdi  funded  by  the  Air  Force  has  result^  in  work  on  the  feasibility  and 
a|^licati(»is  of  stodiastic  methods  in  design  and  analysis  of  protective  structures 
(Twisdale,et  al.,1988;  Ross,  et  al.,1988).  Both  reports  recommend  stochastic  methods. 
Twisdale,  et  al.,1988,  focused  on  develqwg  a  reliability  based  design  for  protective 
structures.  Iheir  main  ^^iroach  is  based  cm  defining  limit  states  fa*  probable  failure 
modes.  Th^  recomniaid  the  use  of  a  First-Order,  Second-Mcxnent  method  (Ang  and 
Tang,1984)  or  Monte  Carlo  simulaticxi  as  discussed  earlier.  Ross  et  al,1988,  recomniend 
the  use  of  Monte  Carlo  simulation  for  assessing  die  influence  of  uncertainties  in  design. 
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Voy  recent  enophasis  in  tfie  conventional  weqxsns  efifects  comnuiity  has  been 
placed  on  increasing  the  ability  to  predict  fiagment  loadings  on  structures  (Hader,  1983; 
Koos,  1988;  Mardiaid,  1988;  Sues  and  Twis(fale,1993).  Ihe  work  in  diis  area  centers 
around  die  ability  to  predict  a  design  fiagment  and  the  loads  generated  fiom  fiagment 
sfjplications.  In  protective  construction  dedgp,  a  deagn  fiagment  is  spedfied  for  i^hkfa 
a  structural  element  must  resist  penetration,  perforation  and/or  qxdling.  In  the  same 
stucty  advances  have  been  made  in  die  area  of  fiagpient  structure  interaction  to  indude 
predicting  die  loads  on  the  structure  due  to  penetrations  versus  ricochets.  \^th  current 
medKxis,  given  a  weqxm  location  relative  to  die  structure,  the  number  of  fragments  that 
penetrate  and  ricodiet  at  each  section  can  be  predicted.  In  turn  the  load  imparted  to  the 
structure  m^  be  predicted  for  either  penetration,  perforation,  cratering,  and  ricochet 
Predictive  tools  under  development  oaiently  lack  sufiQdent  empirical  data  to  allow 
reasonable  verification  of  the  results. 

122  ^iplication  of  Rdudbility  Methods 

The  use  of  reliability  methods  in  engineering  design  has  significantly  increased  in 
recentyears.  Because  ofthe  large  uncertainties  in  most  engineering  problems  the  measure 
of  safety  has  shifted  fium  singly  flying  a  safety  fkrtor  to  devdoping  probability  based 
design  factos  i^di  account  for  the  probability  of  failure  or  survival  in  a  systematic  and 
accqited  fashi(xi  (Ang  &  Tang,1984).  Reliability  problons  aun  to  determine  a 
probability  of  failure  (  Pf )  for  general  oigineoing  problem  Mure  is  easily  and 
traditioially  defined.  In  the  area  of  design  and  analysis  of  protective  structures  die 
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ddSnhkn  of  Mure  is  not  so  well  defined  lliese  stniclures  are  desigtied  to  widisla^ 
certain  level  of  plastic  reqx»se  dierefore  generating  a  nonlinear  proUeuL  Where  that 
nonlinearity  crosses  the  level  diat  tenders  the  structural  dement  under  analysis  uraisdrfe 
is  debatable.  Ihetefoie,  in  lieu  of  determining  a  specific  Pf  it  will  be  more  useful  to 
devdop  a  cuinulative  distributkn  fiaiction  (Q^  for  each  tc^Mnse  parameter  that  inigbt 
rqxesent  a  fiiihse  mode  that  can  be  used  for  later  devdopment  of  one  or  more  P/ s.  It 
has  been  shown  that  these  same  reliability  methods  can  be  used  to  develop  the  desired 
CDFs(Wu,etal,1989). 

The  structural  reliability  methods  alluded  to  dxwe  were  devdoped  by  I&sofer  and 
Lind,1974,  Rackwitz  and  Feisler,1978,  HdientHdiler  and  Rackwit2^1981,  Tve(h,1983  and 
Kfadsen  et  al,  1986.  The  initial  stq)  in  these  methods  is  to  define  a  limit  state  function 
or  performance  function.  A  typical  petfonnance  fiinction  tdoes  the  form: 

g(20  =  V- '  (22) 

in  Adiidi  the  's  rqxesent  a  set  of  basic  state  variables  that  definethe  performance 
dther  e?q)licitly  ch-  inplicitly.  In  general  foe  X/  's  m^  be  correlated  randnn  variables . 
Fcm*  the  purposes  of  this  discussion  th^  are  assumed  to  be  uncorrelated  and  mutually 
independent. 

The  limit  state  may  then  be  defined  as  ^20  -  0,  whidi  for  reliability  purposes 
states  that  if  g(X)  <  0  then  the  system  is  in  a  "safe  state"  and  if  gQQ  >  0  the  tystem 
is  in  a  "Mure  stale."  The  probability  dendty  fimction  (Fl^  of  foe  state  varid)les  is 
defined  asfx(x)  -  (xi  x„ ).  The  cumulative  distribution  function  (COT) 

is  further  defined  as 
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P(t^g)‘FJig)^  f  400*  (2.3) 

This  muh^le  integral  is  in  general  \tty  difScuk  to  evaluate  sinoe  tfie  Auction  gQQ  is 
typically  nonlinear..  The  stnictural  rdialnlity  ooinnaiii^  ovoxxxne  this  difficulty  by 
adopting  and  devdoping  a  variety  of  approxinoate  mediods  for  defining  ptobdxlity  of 
safety  or  in  more  general  terms  the  COF  described  above.  Oie  alternative  is  a  Monte 
Carlo  simulation  ^ch  requires  a  loge  number  of  oomputer  nns.  This  is  fiitfaer 
con^licated  g(X)  is  not  an  explicit  function  as  will  be  the  case  in  the  stnictural 
analysis  of  diis  work.  Odier  {^proximate  methods  are  summarued  by  stating  that  they 
rely  on  die  determination  of  a  most  probable  fiiilure  pmnt  (MPFP)  to  evaluate  the  integral 
above.  The  MPFP  defined  as  the  point  on  die  &ilure  sur&ce  generated  by  nailtqde 
evaluatiais  of  ^x).  that  gives  die  minimum  distance  to  the  origin  of  a  plot  of  that 
surface.  Given  a  plot  of  the  g-function,  at  the  MPFP  that  fioicfion  is  fiather  linearized 
if  it  is  not  Itydefiniticxi  a  linear  function.  Indiecase^heregisnotaneiqilidtfiaiction, 
a  polynotnial  or  reqxMTse  surface  rqxesentationofg  must  be  determined  by  multiple  riBis 
of  the  implicit  code  that  represents  g.  Prior  to  determining  die  MPFP  die  generally  non¬ 
normal  dependent  random  variables,  ,  are  transformed  into  a  set  of  independent, 
standard  normal  variables  Uj . 

NASA  began  a  10-year  program  in  1984  called  PSAM  (probabilistic  structural 
analysis  methods)  for  the  purpose  of  developing  probabilistic  mediods  and  structural 
analysis  codes  for  the  oxiponents  of  cunent  and  future  reusable  ^laoe  propulsion 
tystons.  Southwest  Research  Institute  (Svcdll)  of  San  Antonio  is  one  of  die  key 
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oGOInclors  in  tins  effixt  Uider  this  prognnn  a  prolxMlisdc  ftute  dement  code  called 
NESSl^  (Numerical  Evaluation  of  Stochastic  Structures  Iftider  StressXWii  and 
WiisGhin&l987B)  was  devdoped  by  S^^.  NESSUS  is  designed  ^)ecificdiy  for 
predicting  structural  le^xxise  caused  by  uncertain  basic  variaMes  such  as  loads,  material 
properties,  geonietiy  and  boundaiy  oondhioos. 

SwiU  has  pid)lished  lamierous  papers  uring  what  diey  call  a  iKW,  &5t  probabitity 
analyris  for  tdiability  (Wii  and  \K^rsdiing;1987A).  Ihey  extend  the  FFl  method  of 
Rackwhz  and  Fiessler  ,1978  and  Chen  and  Lind,1983.  Iheir  new  FPI  method  has  been 
demonstrated  to  provide  fist  and  accurate  estimates  of  probability  of  fidlure  for 
performance  functions  normally  associated  with  etigineeringapi^cations.  Ihtyintroduce 
an  improved  schenK  for  constructing  the  equivalent  normal  distributions  and  a  quadratic 
ifjproximation  of  the  original  limit  state  with  log  transformation  options. 

223  Joiit  MmhiotB  Effech  Mnnab 

Ihe  form  of  tiie  weifxxi  delivery  accuracy  adopted  by  the  U.S.  Dqnrtment  of 
Defense  (DoD)  is  outlined  in  the  Jbint  hfunitiom  Effectiveness  Mcnud  (JME3$.  Ihe 
process  of  estimating  force  requirements  stems  from  tiie  estimation  of  die  effectiveness 
of  the  werqxxis  systems  die  DoD  maintains.  These  same  methods  may  be  used  to 
estimate  dw  probabilities  associated  with  die  effectiveness  of  an  adversaries'  weqxxis. 
The  three  conditions  which  lead  to  die  determination  of  the  weqxxi  to  be  used  for  a 
specific  attack  scenario  are: 

1.  Capabilities  and  Characteristics  of  We^xxis  and  Fuses 

2.  Target  Characteristics  and  Vulnerability 
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3.  Weapcn-system  EflSxdvai^ 

Of  tiiese  tlvee  tfie  third  will  be  of  most  impoctanoe  in  this  work.  It  is  assuned 
diat  the  first  two  have  bem  ansidered  by  die  strat^jst  and  have  lead  to  the  sdectkn 
of  a  certain  weapon  or  type  of  weapoa  Iherefore,  the  structural  analysis  will  be  run 
against  a  ^ledfic  weqxn  assuming  it  will  be  delivered  in  the  most  advantageous 
mediod.  Weapon  ddivety  accuracy  is  die  prime  measure  of  weapon-tystem 
effectiveness.  lUs  parameter  will  be  the  prinvay  input  infinaation  and  it  can  be 
influenced  tty  many  unpredictaUe  assodated  varialUes.  Ihe  unoertainties  associated 
widi  die  manuals  estiniates  are  based  on  data  derived  fiom  past  esqperience,  test  results, 
and  dieoredcal  calculations.  The  probalnlity  of  hit  of  a  typical  conventional  weiqxn  is 
characterized  in  die  JMEM  by  the  term  Circular  Enor  ProbaUe  or  CEP.  For  a  large 
luiniier  of  bomb  drops  the  inpact  penis  describe  a  noimaldistribu^  TheprobahUity 
structure  of  die  accuracy  of  a  weapon  may  be  described  in  one  or  tyvo  dimensions.  When 
die  target  on  be  described  by  one  dunenrion  so  may  the  weapon  accurtKty.  An  examine 
would  be  the  bombing  of  a  highway  uhere  die  accuracy  may  be  characterized  by  the 
inpact  distance  to  the  left  or  ri^  of  the  roadway  centerline.  Ihe  jSliEKf  states  diat  die 
crosstrack  and  alongtrack  values  of  bond)  impacts  are  indpendent  and  normally 
distributed  events.  As  suck  die  probfdnlity  of  bodi  events  is  die  product  of  eadi 
individual  event  Ihe  joint  probability  demity  ftmetion  is  shown  as: 


Ino/iy 


-K»  -  l»j)*  /  ♦  0  -  lijJ*  / 


(2.4) 
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\Vliere  die  above  iinctkii  describes  tbe  probability  sinicliiv  0^ 

rectangular  oooidiiiates,  it  has  become  oonvenienC  to  describe  it  in  circular  coordinates. 

Maldng  die  appropriate  transfoimatiom  die  distribution  fiaction  above  na^ 


P  =  1  - 


(2.5) 


with 


=  X*  +  y*  Ojj  =  Oy  *  o  (2.6) 

To  find  die  vahie  of  r  diat  wiU  give  a  dide  containing  SO  percent  of  the  inipactS) 
P  is  set  equal  to  0.5.  Ibis  results  in: 

0 J  =  1  -  (2-2) 

or 

^(-r*/2o*)  ^  (2.8) 

Taking  the  natural  log  and  rearranging,  yields: 

r  *  Ov®(^  =  1.17740  =  CEB  (2-^) 

Ibis  is  defined  as  die  CEP  or  die  radius  of  die  circle  in  ivhich  SO  percent  of  die  impacts 
willoccur.  Wesqxins  and  associated  weapons  systems  or  platforms  are  assi^iedCEPs  for 
strategic  planning  piaposes. 

In  most  cases  die  distribution  is  noncircular  or  Ox  »  Oy .  Hie  delivery  accuracy 
is  dierefore  defined  in  terms  of  two  normal  variates.  Ihey  are  die  range  error,  R,  and 
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die  deflectkn  error,  D.  Hie  range  error  'R*  rqdaoes  the  varial^  'x'  and  die  deflection 
error  *0  rqilaoes  die  variaUe  y.  Range  error  probdile  (REP)  and  deflection  error 
probidile  (I£P)  values  are  derived  in  the  same  manner  as  the  CEP  and  cai  be  found  in 
the  JMEM  for  calculadrig  delivery  aocurades.  The  values  of  REP  and  IXP  are  given  in 
die  JMEM  and  can  be  used  to  back  out  the  first  and  second  moments  of  the  variaUes 
using  die  following  relations  developed  in  die  manial: 

REP  =  0.67450,  DEP  =  0.67450^  (2.10) 

Here  Or  and  Oq  are  the  ddivery-accuracy  measiaes  in  range  and  deflection, 
reflectively.  To  illustrate  diis,  let  die  aimpmnt  be  die  target  and  let  the  fli^  padi  be 
over  die  target  If  two  parallel  lines  oe  (kawn,one  on  either  side  of  the  flight  padi,  at 
a  distance  CO*,  an  infinite  strqi  would  be  formed  that  contained  50%  of  the  impacts  of 
the  individually  aimed  and  delivered  wefxm.  Likewise;,  an  infinite  stifi  similariy 
constructed,  at  90°  to  the  first  str^,  having  2*REP  as  its  over^  widdi,  would  also 
contain  50%  of  the  infiacts.  The  intersection  of  these  two  strips  about  the  target  or 
ainfxiint  will  contain  25%  of  die  infiacts. 

23  Experimental  Test  Des^ 

Eiqieriniental  test  design  or  design  of  eiqierinients  (EX^  is  a  mediod  of  tailoring 
eiqierimaits  in  order  to  extract  the  most  meaningful  information.  The  method  is  very 
useful  when  there  are  many  variables  that  contribute  to  die  design  of  a  single  eiqieriment 
Cte  run  of  a  coinputer  program  may  be  iriterpreted  as  a  single  esqieriiiient  The  mediod 
pronfits  the  user  to  evaluate  the  outixit  requirements  desired  and  ifi-fiont  make  dedsions 
that  can  result  in  fewer  runs,  fiply  sinfiler  analysis  tedmiques  and  yidd  more 
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infonnadon  than  undesigned  Iben^jut  is  generally  nxxeidiaUe  than  the 

enai^oyinent  of  traditional  experimentation  techniques  ^ch  involve  varying  process 
ii^juts  in  an  unsystematic  way.  This  method  often  includes  letting  the  results  of  oiie  nm 
determine  die  ii^iuts  to  the  next  run  i^ch  may  tdl  you  what  you  wait  to  know  or  may 
in  effect  be  very  misleading.  In  contrast,  experimental  design  requires  careful 
considerations  about  all  a^iects  of  yoiar  oqieriment  before  you  make  a  single  run.  This 
mediod  can  save  a  great  deal  of  time  dut  may  be  used  in  imgwtwrf  nns  diat  often  lead 
the  analyst  down  a  long  and  ui^iroductive  road.  >^/hat  assumptions  one  makes  before 
onbaildng  on  a  systonatic  eiqierimental  process  can  be  critical  to  die  qipropriateness  and 
the  validity  of  die  results. 

The  stages  of  die  eiqierimental  design  process  are  as  Mows. 

a)  define  the  problem 

b)  ^lecify  the  model 

c)  select  a  design 

d)  running  the  eigieriment 

e)  analyze  the  data 

0  interpret  the  results 

The  following  is  a  brief  descri|Xuxi  of  each  stage 
a)  Define  die  Problem 

Defining  the  problem  requires  the  oqierimenter  to  bectxne  femilia  with  what  he 
is  trying  to  achieve.  While  defining  the  problon  an  objective  must  be  determined.  The 
factcxs  or  variables  which  will  interact  in  a  specified  process  to  produce  or  accomplish 
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tfiat  objective  are  detemuned  Finaify,  the  response  or  responses  of  die  process  that 
measure  ^^iiedier  the  otgective  has  been  met  are  establi^ied  If  a  representation  of  the 
process  over  ai  entire  operating  range  is  desired  the  response  sur&ce  metfiodtdogy  (RSM) 
is  reccmnended. 

b)  l^iecify  die  Mbdd 

Proper  modeling  of  the  process  in  order  to  stuefy  die  interactions  and  contributions 
ed*  each  &ctor  to  die  observed  reqxnse  is  the  key  of  the  method.  Ihe  model  wiU  be  a 
generic  equation  that  will  be  used  to  predict  die  outcome  of  the  process  after  initial  data 
has  been  collected  and  analyzed  The  model  may  be  eidier  linear  or  more  cooaplex. 
Throu^  an  analysis  of  a  limited  set  of  actual  process  runs  die  ooefiBdents  of  the  model 
equation  are  delennined  The  narherofactual  process  runs  required  is  determined  based 
on  the  aelectod  deaiin  method 

c)  Select  a  Dcai^ 

The  design  is  the  niediod  of  determiiiiiig  the  actual  runs  required  of  an  esqieriinent 
to  achieve  the  demed  oliyective.  The  design  and  its  qiecific  fonn  are  determine  on  the 
overall  objective,  die  nundier  and  Qpe  of  factors  and  die  nder  of  the  terms  in  die  model. 
Classical  and  optimal  are  the  two  general  cat^ories  of  designs.  The  following  are  a  few 
types  of  classical  designs: 

>  Plankett-Burman  -  Fractiaial  Factorial 

>  Full  Factorial  •  Box-Bdmken 

•  Central  Composite  >  Orthogcaial  Arr^ 

These  classical  designs  are  generally  inflexible  on  die  number  of  tuns  and  variable 
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settings  reqpnred.  Optinml  (fesigns  are  less  oonvendonal  and  use  oomplex  algoritfams  to 
determine  a  reduction  of  runs  and  variable  settings.  Although  fewer  runs  are 
aooonqjlished,  the  interpretation  of  the  results  is  often  inore  difficult  One  of  the  most 
common  optimal  design  is  the  D-optimal  desigp.  The  D  in  D-optimal  refers  to  a 
"determinant”  manipulation.  The  maximizing  of  die  value  of  die  determinaitt  of  die 
co^cients  is  equivalent  to  minimizing  die  uncertainty  in  die  coefiBcient  values.  RSM 
erqperiments  generally  require  classical  methods  due  to  their  complex  ooefiBdents. 

d)  Run  the  Expaimmt 

Once  the  number  of  runs  is  determined  along  with  die  desired  variable  values 
associated  with  each  run  a  worksheet  is  developed  and  die  results  are  cataloged. 

e)  Analyze  the  Residts 

Regression  techniques  are  used  to  fit  the  data  to  the  model.  Once  the  initial 
results  are  fit  the  model  is  evaluated  and  refined  until  the  model  results  are  reasonable. 

f)  Interpret  die  Residts 

With  a  satisfactory  model  it  can  be  used  to  learn  more  about  the  basic  process. 
In  this  case  the  model  will  be  used  to  determine  die  probability  structure  of  the  reqxHise. 
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CHAPIER3 

ReqMMBe  l^atistics  Granation 

3.1  Production 

Ihis  du^pter  describes  die  development  of  load  and  stnjctural  models  which 
together  generate  response  statistics  for  typical  hardened  &cilities  subjected  to 
conventional  wet^xxis  effects.  The  load  model  is  derived  fitxn  coded  blast  environment 
prediction  formulas  and  gn^ihs  found  in  the  literature.  The  model  presented  extends 
current  methods  tty  inaxporadng  die  ability  to  jxedict  loads  on  structures  in  both  time 
and  ^lace.  A  typical  wall  uiit  is  modeled  as  a  condnuous,  hysteretic  beam  in  order  to 
accept  the  ^xttial/tai^xiral  load  model  iipit  and  produce  reqxxise  statistics  across  its 
hei^  The  loading  parameters  as  wdl  as  die  structure  equations  of  state  are 
nondimensionalized  within  the  code  to  allow  investigations  of  various  structural 
cc»ifigurati(His  with  a  minimum  of  ii^xit  The  equatims  of  state  of  die  structure  are 
solved  by  the  method  of  weighted  residuals,  cxi  an  itoadve  basis.  The  combined  load  and 
structure  models  is  termed  the  Re^ionse  Statistics  Generator  (RSG)  mode.  RSG  model 
results  are  conpared  to  full  and  sub-scaled  test  results  ffcxn  die  litaature  (Pyle  and 
Baber, 1991). 

The  following  statemoits/assun^ois  are  made  as  a  pre&ce  to  die  rmiainder  of 
this  ch^er: 

1.  The  attack  scenario  posits  the  use  of  unguided  munitirxis  against  the  target 
These  unguided  munitiois  are  ncMiiially  dqiloyed  in  "sdcks”  fiom  four  to  eight  weqxms 
po*  airoraft 
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2.  A  wall  of  the  above  ground  hardened  structure  considered  in  this  study, 
subjected  to  blast  and  fiagment  loadings  frran  a  conventional  wes^xm  detonated  in  its 
vicinity,  will  respond  as  a  one-way  slab  (Biggs,1964;  Coltharp,1985;  Hyde,1988).  This 
assun^cxi  is  based  on  die  &ct  diat  a  typical  hardened  structure  wall  has  a  length  to 
hei^  ratio  which  allows  for  this  characterization.  'Therefore,  typical  wall  units  will  be 
described  as  a  series  of  beam  models.  'The  beam  model  will  rqiresent  a  section  of  the 
structure  described  as  the  region  of  interest  (ROI)  in  this  work  (Figure  3  4). 

3.  The  we^XHis  to  be  investigated  \rill  be  general  purpose  btxnbs  which  will  be 
assumed  to  detonate  at  a  near  votical  oientation  or  perpendicular  to  the  ground  sur&oe 
with  the  nose  toudiing.  Casings  effects  will  not  be  considered  esqilicitly  but  will  be 
included  by  using  the  total  wet^xm  wei^  (eiqilosive  plus  casing  weight)  for  calculating 
blast  loadings. 

3.2  Loading  Model 

The  blast  oivironment  from  a  surface  detonatiai  of  a  typical  omventicxial  we^xxi 
is  depicted  in  Figure  3.1.  In  this  Figure,  taken  from  TMS-85S,1986,  the  commcm 
approximation  that  the  load  is  tg^lied  unifcxmly  at  the  structure  is  dq)icted  when  in  fact 
it  varies  rapidly  over  the  hei^  of  the  wall.  In  order  to  accurately  rqxesent  die  loading 
of  a  strucmre  from  the  near  field  detonation  of  a  ctxwendcxial  wetqxxi,  a  mo^l  is 
required  that  can  represoit  both  the  spatial  and  temporal  diaractoistics  of  this  fdienomena. 
The  free-field  pressure-time  histcxy  from  ccxiventiixial  weqxxi  debxiation  is  shown  in 
Figure  3.2,  also  freon  Th/B-855.  The  pogram  ConWep,  a  conventional  wespois  effects 
code  develqied  by  the  Amty  Corp  of  Engineers,  Waterwitys  Eiqperiment  Statitoi 
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Figure  32  Typical  Free-Field  Pressure-Time  ffistoy  From  a  ConvendKxial  Weapm 

EtetcHiaticxi  (TM5-855,1986) 
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(Hyde,1988),  uses  the  physics  described  later  in  this  seetkn  to  predict  die  blast 
environment  from  a  weqxm  detonation.  The  incident  pressure-time  histoiy  seen  at  the 
base,  die  mid-heig^  and  top  of  the  ROi  are  shown  in  Figures  3.3  dvoug^  3.5.  Clearly 
die  blast  seen  die  wall  is  not  uniform.  This  variation  is  further  exaggerated  by  die 
hitler  reflection  and  magnification  of  die  blast  at  the  lower  angles  of  inddenoe  for  the 
lower  p(xti<xis  of  the  wall.  It  is  argued  diat  die  variations  are  not  critical  for  calculating 
die  bending  response  of  the  section,  hofweva  tfaqr  seem  to  be  critical  frr  predicting  shear 
fidlures  away  from  the  si^ipcxts  as  was  discussed  in  Section  2.1. 

The  principal  parameters  diat  can  be  taken  from  Figure  3.1  are: 

(1)  Time  ofarrival,afier  detonation,  ofdie  blast  wave  to  die  structure.  (t„msec) 

(2)  Positive  pressure  t^iase  duration  of  blast  afier  arrival  at  structure.  (fo,nisec) 
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(3)  Peak  positive  incident  pressure  seen  by  the  structure.  (Pa^psi) 

(4)  Positive  incident  inpilse  delivered  to  the  structure,  (i^  nisec*psi) 

These  parameters  are  presented  in  grt^cal  and  ecpiation  form  in  Kingery,1966. 
Figure  33  is  an  extract  from  Kingety,1966  and  was  constructed  uang  the  results  of 
e7q)enniental  tests  accrenplished  widi  hemispherical  sur&ce  bursts  of  payloads  ranging 
fixxn  S  to  SOO  tons.  Qibe  root  scaling  and  altitude  corrections  were  made  to  bring  die 
results  to  standard  sea  level  conditions  for  a  1  pound  TNT  charge.  The  curves  were  then 
generated  using  standard  curve  fitting  techniques.  The  required  iqxit  for  pulling 
information  frrxn  this  gnqdi  is  die  scaled  range  of  the  detonation.  This  is  a  condnnation 
of  the  locatiai  of  die  deUxiaticxi  relative  to  the  frdlity  or  range  (R)  and  die  TNT 
equivalent  throw  weight  (W)  of  die  charge.  The  scaled  range  (SR)  is  defined  as: 


Figure  3.4  CcxiWep  Plot  of  Pressure-Tlnie  Ifistory  at  the  h^d-lfeight  of  die  Wall 
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SR  = 


1^1/3 


(3.1) 


Ihis  so  called  cube  root  scaling  is  attributed  to  Hopkinsan,191S,  and  is  well  documented 
and  aooqjted  in  the  literature  (Baker,1973;  Kingeiy,1966). 

The  weqxxi  used  in  this  work  is  assumed  to  be  a  general  purpose  bomb  detonating 
on  impact  at  a  near  perpendicular  angle  to  the  surftoe  of  die  ground.  Eiqilicit 
casing^fiagment  loading  effects  be  disregarded  but  will  be  accounted  for  by 
increasing  the  TNT  throw  wei^  of  die  weqxm  by  die  casing  weight  The  exduskn  of 
the  casing  effects  is  not  meant  to  infer  its  unimportance  but  is  an  initial  sim|difying 
assunpticxi  deemed  necessary  to  set  die  method.  As  discussed  in  Section  2.1,  there 
is  a  great  deal  of  intoest  in  the  community  for  furdier  study  of  case  fiagtnent  loadings 


Figure  3.5  QniWep  Plot  of  I¥essure-Time  rfistoiy  at  the  Top  of  the  Wall 


Figure  3.6  Positive  Phase  Shodc  Wave  Parameters  for  Hemi^iherical  INT  Bursts  at 

Sea  Level  (1^5-855,1986) 
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andthesynergistkeffeaofltoaadfit^^  llie  toinokigy  to  predKt  suc^ 
is  not  available  at  present  and  in  fiict  the  new  tn-serv4oe  oonventknal  weaixins  efifects 
design  manual,  under  development,  uses  the  above  assumption.  The  underlying  premise 
in  describing  the  distribution  of  fiagments  from  a  cased  ordnance  is  based  on  a  probability 
law  and  should  lend  itself  to  inclusion  in  sid)sequent  work. 

Given  the  attade  scenario  above,  die  burst  location  idative  to  die  RCX,  Mhicfa  will 
be  furdier  defined  in  Section  3.4,  m^  be  ddined  Ity  the  range  (R)  and  angie  of  inddenoe 
(a),  uhidi  are  functions  of  die  lateral  and  incident  standoffs,  x  and  z  reqiectivdy  as 
shown  in  Figure  3.7. 


Segments  of  die  program  Con  were  modified  aid  incoqxxated  into  the  load 
models  presented  here.  ConWep  takes  the  weqxxi  TNT  equivalent  dsow  wei^  and 


35 

detonatkii  location  and  calculates  die  saded  mge.  Widi  die  scaled  range  ConWep 
utilizes  equations  derived  fiom  the  curves  of  Figure  3.6  to  generate  the  parameters  vdikh 
may  feed  the  eni|»rically  based  mxfified  Friedlander  equation  (Baker, 1973)  to  produce 
the  ideal  blast  pressure-tmie  histcxy  in  air.  Ihe  modified  Frieiflander  equation 


*•  f* 

measures  time  fiom  the  time  of  arrival  of  die  blast  wave  to  the  point  in  space  of  interest 
Various  forms  of  this  equation  are  found  in  die  literature  and  range  fiom  triangular 
representation  with  two  variaUes,  to  this  exponential  etpiadon  widi  dvee  variaUes  to 
mult^le  eiqxinential  equations  with  five  variables  (Baker,1973).  Ihe  cfamoe  of 
empiiicalrqiesentation  is  cornmensurate  with  the  accuraty  desired  but  diat"probaW^  die 
best  ccMnpmnise  is  the  modified  Friedlander  equation,  since  it.allows  at^ustment  to 
conform  to  the  most  important  blast  wave  properties  and  yet  is  not  too  oomplex," 
(Baker,1973). 

When  die  blast  wave  hits  a  rigid  surfece  die  particle  velocities  ate  arrested  and  the 
pressure,  den^  and  tenqxrature  are  increased  tdxive  the  values  in  the  fiee  field  incident 
wave  (Baker,1973).  Ihe  magnification  resulting  fiom  die  reflection  of  die  blast  wave  is 
a  function  of  the  incident  wave  pressure  and  die  an^e  of  incidence  a.  As  defined,  the 
lower  the  angle  of  incidence  die  hi^ier  die  magnification  fector.  Magnification  fiictor 
curves,  as  seen  in  Figure  3.8,  have  been  developed  and  are  ctxnmon  throu^iout  die 
literature.  Ihese  curves,  which  were  generated  fiom  eiqperiniental  data  (Qrawford,1971), 


Figure  3.8  Reflected  Pressure  Coefficient  Versus  An^  of  Incidence 

Cnvi5-855,1986) 
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are  coded  in  ConWq>  and  were  also  inoorpontol  into  die  toad  modds  Two 

load  model  rqjresentadons  were  generated  and  investigated  in  this  work.  Hie  goal  of 
each  modd  was  to  best  rqxesent  the  nonuniformity  of  die  (fynamic  loading  on  a  stnicture 
from  a  near  fidd  conventional  weapon  detonation  and  lend  itsdf  to  interfradng  widi  die 
stracture  modd  i^iich  was  devdoped. 

3JL1  -  Tmiporal  ^^hve  Fonn  Modd 

In  Cotdiaq)  et  al, 1 985,  a  scaled  test  stnicture  was  subjected  to  the  near  field  (5  to 
15  feet)  detonadon  of  a  ^ical  conventional  IMT  throw  wd^  (250  to  1000  lbs)  wetqxxL 
Hie  test  structure  was  instrumented  with  pressure  gniges  as  shown  in  Figure  3.9.  Hie 
plot  of  the  average  pressures  at  various  times  along  the  vertical  gage  line  is  shown  in 
Figure  3.10.  Hiese  series  of  plots  represent  a  typical  set  of  blast  profiles  which  would 


Figure  3.9  Loc^oi  of  Pressure  Gauges  for  Test  Series  1  of  Cohharp,  et  al,  1985 


38 

be  produced  by  a  conventional  wn^xn  in  die  INT  Ifarow  wei^  and  standofif  ranges  of 
interest  in  this  stucfy.  Ibe  rqxoduction  of  diese  profiles,  while  varying  die  tnaximuin 
pressure  muhiplia,  represents  a  Q^pical  spatial  and  temporal  loading  serpience  needed  for 
ii^  into  die  continuous  beam  model  devdoped  in  section  32.  l*ierqxoductionofdiese 
profiles  for  the  purpose  stated  will  be  called  the  Temporal  Wave  Fonn  (TWF)  load  model. 

In  development  of  the  TWF  load  model,  the  plots  of  Figure  3.10  were  reproduced 
by  mult^e  orderpofynomialsastiiowninFigure3.il.  Basedonananafytisofpresaire 
time-histories  for  specific  wespons  enplr^ments  (Coldiatp  et  al,198S;  Wiri^  et  al,1988; 
Ifyde,1989)  and  running  numerous  simulations  with  ConWqi,  the  number  of  profiles 
required  to  represent  the  tenporal  and  spatial  load  history  associated  with  the  scaled 
ranges  ( range  (R)  ova  the  cube  root  of  die  INT  throw  wei^)  of  interest,  oe  tiiown 
Table  3.1.  These  allocations  are  based  on  the  &ct  that  die  fiatha  finm  the  structure  the 
we^xHi  detcHiates  the  moe  the  load  profile  resohbles  a  plane  wave. 


BHHHIHHiiHBIHHH 
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lable  3.1  Scaled  Kan|^  Versus  l^nies  Kequired 


The  prime  indicate  fa  allocating  profiles  and  assigning  (xessure  multipliers  is  a  dbeda 
of  the  total  inpulse  inparted  to  the  structure  fa  a  given  loading  scenario.  Ihe  inpulse 
inparted  to  die  beam  elonent  is  conputed  ty  summing  the  impulse  ppUed  by  eadi 
profile.  The  individual  profile  inpulse  is  determined  by  integrating  each  profile 
ployn(»nial  ova  its  assigned  duration.  The  generated  total  inpulse  values  are  conpared 


Figure  3.10  Pressure  Distribution  Along  Vertical  Gauge  Line  at  Various  Times 

(Cohhatp,  et  al,1985) 
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to  iiniMlse  calculations  i^tumed  by  GonWep  divides  a  structiffal 

element  into  sections  and  collects  the  in^xilse  values  shown  in  Figure  3.6  associated 
widi  the  scaled  range  to  each  section.  Table  32  shows  the  results  of  exhaustive 
ccHipdsons  of  die  type  depicted  above.  Ibis  table  shows  die  observed  relation  between 
the  number  of  profiles  required  and  the  percentage  of  the  maximum  reflected  pressure  and 
positive  jdiase  duration  allocated  to  each  profile. 

Given  a  scaled  i»ige,  the  TWF  modd  determines  die  required  number  of  profiles 
fitxn  Table  3.1.  The  TWF  model  then  takes  die  calculated  maxinnim  reflected  pre^ure 


Figure  3.1 1  Goieric  Loading  Profiles  Generated  Using  the  Temporal  Wave  Form 

Model 


ML 
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and  positive  phase  duration  and  ^lies  them  to  the  ^ipropriate  set  of  profiles  per  TaUe 
32.  The  RSG  program  calls  the  TWF  model  at  eadi  time  stq)  and  ^lies  die 
qjpixipnate  load  to  die  stnx:ture  model.  Figure  3.11  is  a  plot  ofa  case  where  five  profiles 
were  required  to  identify  the  loading  based  oi  die  scaled  range. 
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labie  3.2  Frotiie  Allocations 


32.2  -  Spatid/Temporal  Modified  FMedander  Equation  Model 


Hie  ability  to  vary  the  load  spatially  is  added  to  the  modified  Friedlander  equadcxi 
by  making  die  primary  iiqiut  parametos  funcdois  of  the  vertical  locadrxi  of  the  point  of 
interest  relative  to  die  burst  locadoL  The  time  of  arrival  and  positive  fdiaseduratiai 
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(tg)  become  functions  of  tiie  three  dimenskinal  scaled  range  (IDSR)  defined  in  equation 
3.3,  A^ch  includes  now  die  distance  up  the  R£X  (y)  from  the  ground  sur&ce.  The 
incident  pressure  (P^)  also  becomes  a  function  of  die  TDSR  .  The  resulted  incident 
pressure  is  magnified  die  tqyvoprtare  reflection  fiKtor  firan  Figure  3.8  using  as  iiput 
the  angle  of  incidence  calculated  using  die  x,y  and  z  coordinates  of  the  point  of  interest 
The  irqnit  needed  to  generate  die  required  parameters  at  ai^  point  along  die  ROI  are  the 
detonation  location,  x  and  z,  die  distance  vp  the  ROI,  y,  and  the  iwapon  throw  wd^ 
W.  The  TDSR  to  any  point  along  die  RCH  is  dierefore: 


IDSR  =  (33) 

WV3 

With  this  scaled  range  one  can  go  to  the  curves  of  Figure  3.7  and  read  o£f  die 
parameters  identified  above.  The  ability  for  ^latial  variation  is  coupled  with  the  inherent 
tenporal  variability  of  the  equatic»i  resulting  in  the  ^latial/tarporal  modified  Friedlando* 
equation: 


f  f-fjy)  1 


(3.4) 


The  value  P^)  is  the  reflected  pressure  seoi  at  the  point  y  ip  the  wall,  at  time  t  The 
Heaviside  function  H[t  -  t,(y)]  ensures  the  loading  at  a  ^lecific  point  is  not  tpplied  prior 
to  its  time  of  arrival. 


As  part  Hat  ncndunensionalizadoii,  tibe  hd^at  of  the  ROI  is  nonnalued  and  it 
is  given  a  unit  width.  Parameter  are  calculated  at  70  locations  the  wall  with  50 
Ming  on  the  lower  half  of  tiie  wall  >Ahere  the  loading  dianges  most  r^dly.  At  each 
location  an  equation  is  generated  and  die  parameters  stored  For  a  typical  12  foot  high 
wall,  loaded  by  a  1000  pound  bomb  at  a  15  foot  perpendicular  rtandofi^  Figures  33-3.5 
show  die  incident  pressure  time  histories  at  die  base  of  die  wall,  die  mid  point  and  the 
top.  These  figures  were  generated  in  the  program  ConWqi.  The  reflected  pressure  may 
be  calculated  by  multiplying  die  incident  pressures  on  diese  curves  by  the  angle  of 
incidence  multiplier  taken  fixxn  Figure  3.8.  At  each  integration  time  step  the  model  pulls 


Figure  3.12  Load  Profiles  Created  widi  the  Spatial  and  Tonporal  Friedlander 

Equatioi  Model 
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tiie  airent  knd  cn  tibat  segment  of  the  from  the  apixxifaiale  cuve. 

Hie  Friedlander  equation  model  can  also  be  used  to  generate  profiles  similar  to 
die  ones  recorded  in  Cohfaaip  et  al,1985.  Figure  3.12  shows  profiles  generated  fixxn  die 
Friedlander  model  using  die  test  parameters  fixxn  Coldiaip  et  al,198S  as  ii^xjt 
33  Load  Model  Selection 

For  die  purpose  of  diis  study  die  Spatial/Temporal  Mcxfified  Friedlander  Equation 
leant  itself  weU  to  alteration  into  die  structure  modd  and  solution  method  diosen.  The 
TWF  model,  though  deemed  ipplicable,  did  not  lend  itself  to  die  wei^ited  residual 
method  used  to  solve  die  equations  of  state. 

3.4  Stnictuie  Modd 

Reinfoii^d  ccxicrete  waU  structures  subjected  to  the  conventional  weapons  effects 
of  interest  in  this  woiic  have  been  shown  to  re^xmd  essentially  as  one  way  slabs, 
therefcHie  the  use  of  a  beam  model  to  rqiresent  die  re^xxise  cfaatacteristics  is  justified. 
It  is  desired  to  model  the  beam  so  that  die  re^xxise  can  be  evaluated  along  its  aitire 
length  in  such  a  way  as  to  include  the  higher  modes  of  vibration.  It  is  also  desired  diat 
the  model  not  be  too  detailed  as  to  cause  unwanted  ccHiiplexity  requiring  long  and 
eiqiensive  conqnita  runs,  as  in  a  ntxilinear  finite  element  analysis.  The  nonlinear 
response  of  the  structure  is  therefore  modeled  as  a  continuous  hysteretic  beam.  Thebeam 
model  is  takoi  to  rqiresait  die  vertically  centered  section  of  the  wall  referred  to  as  the 
ROI  in  Figure  3.7. 

The  support  c(Hiditi(xis  diat  best  rqxesent  the  re^xHise  of  similar  structures  has 
been  characterized  as  being  somei^here  between  fixed  and  sinply  supplied  (Ross  and 
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KiHwii]kler,1985).  To  gain  insig^  into  Ifae  method  proposed  the  exam|^  dev^jped  in 
dtis  work  assumes  sin^)ly  sifjpoited  conditions.  This  assumption  however  is  not  unique 
as  can  be  seen  in  the  literature  (Colthatp  et  al,198S;  Biggs,l%2).  Other  sippoit 
conditions  could  be  incorporated  into  the  model  at  a  later  date. 

If  the  system  is  modeled  as  a  general  beam  die  ecpudon  of  motion  of  the  system 
can  be  written  as 

+  pAv^  *  (3-5) 

where  the  ccanma  and  subscript  letters  rqxesent  die  derivative  of  that  term  widi  reflect 
to  the  subscripted  letter.  The  moment  curvature  relationshqi  can  be  defined  in  terms  of 
linear  and  nonlinear  pcxticais.  The  mcmient  term  of  equation  3.5  can  be  written  as: 

mOcJ)  .  *  (1-iVM  \A,\sl  (3.6) 

The  torn  M  rqxesents  the  nonlinear,  hysteretic  portion  of  the  total  moment  Ihe  term 
Aq  can  be  intopreted  as  the  po^-yield  to  pre-yield  nxxnent-curvature  ratio  or  for 
reinfcxced  conoete,  the  post-yield  to  inital  tangent  ratio  and  it  controls  the  degree  of 
nonlinearity  that  the  system  will  exhibit  (see  Figure  3.13).  For  exan^ile,  if  Aq  equals  1 
the  system  is  effectively  linear  vriiile  if  Aq  equals  =  0,  the  system  is  fully  nonlinear 
hystoetic. .  Substituting  equation  3.6  into  equation  3.5  returns: 

(l-i4o)M,„  +  +  pi4v,„  =  qJU)  (3-7) 

The  term  qj(x,t)  rqxesents  the  load  model  iipit  to  the  system.  The  hysteretic  moment 
curvature  relationship  may  be  eapessed  in  various  w^.  One  particularly  convenient 
form  is  the  rate  type  anooth  hysteretic  tyston  attributed  to  Bouc(1967), 
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I  M  m-i  M  _  A  I  M  in 

AfJ  M.  ^4),  'mJ 


(3.8) 


Here  is  tiie  ultiiiiate  moment,  howeva,  4;  is  the  curvature  viiere  the  initial  tangent 


curve  intercepts  tiie  horizontal  ultimate  moment  curve.  A  hj'steretic  curve  is 
shown  in  Figure  3.13.  The  terms  A|,  As  and  n  of  equation  3.8  contrcd  die  of  the 

hysteredc  curve  imrder  to  rq)licate  the  reqxxise  histoo'  of  a  real  world  system  (see  Baber 
and  Wen,1979  for  studies  on  hysteresis  diape  oontnd).  This  and  similar  modds  were 
developed  to  represent  re^xmse  of  structures  to  seismic  type  (tynamic  loadings  which 
would  normally  take  it  dirou^  several  cydes.  The  model  may  also  be  modified  to 
include  system  degradation.  At  this  time  diere  is  not  suffident  data  to  generate  an 
appropriate  degradation  iiiodel  for  the  ty^  of  loadings  uiideroondderationherdtL  The 
motion  of  die  system  is  now  totalty  defined  by  equations  3.7  and  3.8.  Closed  form 
solution  of  diese  equations  is  difiScult,  so  a  reasonable  ipproximate  solution  is  desired. 

PriOT  to  develqiing  a  solution  technique  the  equations  are  nondimensionalized. 
This  will  allow  the  analysis  of  ^'sterns  widi  various  reinforcing  and  concrete 
configuraticHis  by  manipulation  of  cxily  diree  input  variables.  The  eolations  become,  after 
rearranging  to  solve  for  the  derivative  of  the  di^lacement  with  reflect  to  time: 

~  4(5*^)  ~  (1  ~  ~  (3.9) 


The  details  of  die  nondimensionalization  are  in  die  ^ipendix  A 


(3.10) 
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Figure  3.13  Typical  Moiiiert-<Xirvature  Relationship 


Figure  3.14  Typical  Smootfi  Hysteredc  System 
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33  Snh^inM  TVImA|m> 

The  approach  chosen  to  solve  these  equaticns  is  tiie  method  of  wei^ited 
residuals(Cunninghain,1958).  Following  this  methodology,  a  set  of  complete  functions 
tfiat  satisfy  die  physical  boundaiy  conditioos  are  chosen  to  approximate  the  response  of 
die  system.  The  diylaocment  of  the  system  is  approodmaled  by  a  finite  sum  tenns  in 
die  form 

v(5,t)  =  £  ej(T)Tj(0  (3.11) 

*•1 

where  the  complete  set  of  finctions  4^^  are  chosen  to  be  the  set  of  dgen-functions  dtf 
satisfy  die  assumed  end  conditions.  In  a  like  manner  die  moment  solution  ms^  be 
approximated  as  a  series  taking  die  form: 

*|i 

=  E  ■  P  '2) 

!•! 

Substituting  the  assumed  solutions  of  the  associated  linear  ^sterns  into  the  differential 
equations  results  in  a  residual  difference  e.  The  residuals  of  equations  3.9  and  3.10 
beccxne; 


(3.13) 


and 
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ft-i  1-1  (-1 

-'*»(E  »«.«*.)  IE  •/?/!■ 

*-l  1-1 

Introdudiig 


“lit  *  Yi  **1A*  VI, 


(3.14) 


(3.15) 


ten  substituting  tiiis  variable  into  equatkxis  3.13  and  3.14,  results  in  three  first  order, 
partial  differential  equations.  Followii^  te  weigjhfted  residual  metiiod,  te  error  in 
equations  3.13  and  3.14  are  oilfaogcxializied  with  reflect  to  a  set  of  weighting  fiaictions 
witii  te  equations  taking  te  form 


0 

1 

fea(M*i((V(-0  (3.17) 

0 

Where  e^,  and  e,^  denote  te  residual  equations  3.13  and  3.14  above.  The  modd  solves 
equations  3.15,  3.16  and  3.17  simuhaneously  for  te  ooefiSdents  oc,  p  and  y  at  specific 
times  t.  Ihe  total  nunter  of  equations  to  be  solved  is  equal  to  twice  te  number  of 
terms  used  in  equation  3.11  phis  te  number  of  terms  used  in  equation  3.12.  Hie 
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oomplete  set  of  functions  used  as  the  appraximate  solution  also  represent  the  mode 
shq)es  of  a  lineody  reqpondiiig,  singly  si:;]p(xted  beam.  It  is  interesting  dust  even  with 
the  nonlinearity  of  this  pix)blem  significant  "modal"  re^ionse  occurs  when  the  detonation 
location  is  close  to  the  wall.  Once  the  coefficients  are  determined,  die  displacements, 
moments,  shears  and  curvatures  at  any  point  along  the  beam  and  at  any  time  diving  die 
reqxmse  histoiy,  can  be  evaluated. 

Input  variables  associated  widi  die  modd  ars  presented  in  table  33.  Results  of 
die  model  under  a  distributed  unit,  dynamic  load  were  compared  to  a  single  d^rcc  of 
fieedom  (SD(^  and  a  finite  elenxnt  representations  of  the  systenL  Ihe  SDCf  analysis 
was  tun  using  a  program  endded  Biggs  (Baker,19S9).  The  finite  element  model  consisted 
of  10  two-dimensional  beam  dements  using  die  program  ANSYS,  (Swanson  Analysis 
Systems,  1989)  which  represented  half  the  beam  recognizing  symmetry.  Deflection 
results  at  various  points  along  the  beam  were  compared.  Figures  3.15-3.17  show  a 
(xxi^iarison  of  maximum  deflections  at  duee  points.  Maximum  values  were  within  2%  of 
each  other  for  these  locatiixis. 

The  NONUN  model  required  seven  alpha  and  beta  terms  to  ensure  convergence. 
The  influence  of  the  higher  order  terms  is  most  evident  in  die  re^xxise  away  fixim  the 
midspan.  . 
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Model  Variables 

Variable  Name 

Deagjptiflp 

n,  AO,  Al,  A2 

Hysteresis  shape  &ctcxs 

N^andNy 

bbmber  of  inoment  Ot)  and 
di^daoement  (v)  tenns  in  the  sdudcn 

TS/Iterarions 

Hme  stq>  and  heratioos 

StmctimVsriribia 

^andl^ 

Steel  yield  and  concrete 
compressive  strengdi 

RHOS(qJ 

Pereentage  of  steel  leinfixoement 

XL,d 

Wall  height  and  duckness 

GAMC(Ye) 

Density  of  concrete 

Ddiveiy  and  Load  Variables 
ZVAL,XVAL 

ltiin0B  anH  rfiafawyiffg 

oenteriine  (d*  the  RCX 

WVAL 

TNT  equivalent  tfaroiw  of  weqxn 

Table  33  NONUN  Variables 


3.7  Effect  <^VaiyiiCtfieNiiiiier  of  Tmn  in  the  S(rfalk» 

Hie  effect  of  tiie  number  of  terms  N^and  H,,  used  in  the  solution  was 
investigated  for  a  typical  loading  scenario.  Assuming  die  validity  of  the  qiproximation 
incre^es  widi  die  number  of  terms  used  (Cunnin^iam,19S6),  ai  attend  was  made  to 
find  the  point  where  the  change  in  reqxmse  became  minimal  with  an  increase  of  terms. 
Using  the  desi^  of  experiments  (DOE)  methodology  and  die  progt:am  RS/Discova 
(BI^1992AX  an  experiment  was  designed  to  show  the  effect  of  varying  the  number  of 


AO  =  0.05 

o 

II 

< 

A2  =  0.5 

n  =  5 

TS  =  0.01 

Iterations  =  100 

fc  =  5000  psi 

&  =  60,000  psi 

ps  =  0.0(^5 

70  =  l^pcf 

XL=120  inches 

d  =  30  inches 

WVAL  =  20001bs 

ZVAL=10feet 

XVAL  =  0 

Table  3.4  NvNu  Ejqxrimoit  Variable  Settings 
and  K,  terms  on  the  response  statistics  generated  the  NONUN  model.  The 

paranKters  in  Table  3.4  were  ised  in  the  analysis.  Ihe  DOE  respcxise  surfece 
medKidoIogy  was  chosen,  with  a  Coitral  Conposite  Face-Centered  (CCF)  design  The 
resulting  relaticMiship  was  modeled  with  a  quadratic  equatiem.  The  CCF  design  required 
13  runs  to  characterize  the  response.  The  results  showed  that  die  deflection  and  curvature 
predictions  wore  not  nearly  as  sensitive  to  variatiois  in  the  number  of  terns 
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as  the  shear  aid  momaitpredicdons.  The  deflectim  and  shear  resuhs  are  ^lown  in  TaUe 
3.5.  Figure  3.18  is  a  plot  of  die  shear  model  developed  RS/Disoover  fixxn  the  data 
of  Table  3.5.  The  plot  is  of  shear  predictions  for  Nv  etpial  to  1  and  10  versus  a  range  of 
Na  The  95%  simultaneous  confidence  bounds  idiown  have  die  proper^  diat,  at  most,  5% 
of  the  models  predictions  will  have  a  single  confidence  interval  diat  does  not  contain  the 
true  estimate  for  the  prediction  (BBN,1992B). 

3.7  Oniiinrisonof  Response  Modd  Iti  Ibst  Resulli 

The  combined  oxiventional  wesson  load  model  and  sonihaidened  structure  model 
was  evaluated  ccxi^iarison  of  deflection  results  from  the  qiplication  of  actual  weqxxi 
loadings  (xi  structures  for  test  results  are  available  in  die  literature. 

The  tests  series  used  fix- (xxrpnrison  are  reported  in  CoUharp  et  al,1985,  Wri^  et  al,1987 
and  lfyde,1989.  Table  3.6  gives  the  results  of  the  model  simulations.  The  standoff 
referenced  are  incident  values.  Test  di^laconaits  were  recoded  i^ing  active  gages.  The 
diarge  weight  and  standoff  locations  fo*  the  Hyde  and  Coldiaip  tests  are  not  presented 
due  to  security  classificatioL  AltlKxigh  not  fully  validated,  diis  comparison  diows  that 
the  NONUN  code  returns  reasoiable  {xedicticxis  fo*  the  iiqxit  variable  ranges  of  intoest. 
The  re^xxise  statistics  goioator  (NCMJN)  returns  a  cnide  estimate  of  the  total  shear 
response  that  could  be  enhanced  by  inclusioi  of  a  Timoshenko  beam  element  in  the 
model.  In  a  like  manno',  the  inclusicxi  of  fiagmoit  effects  may  be  included  at  the  oude 
level  of  the  current  state  of  tedmology  in  this  area  These  two  additiois  should  add  to 
the  procedure's  ability  to  predict  shear  kills  explicitly.  As  stated  in  the  background 
informati(xi  there  are  numoxHis  researchers  working  (xi  die  problem  of  the  synergistic 
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efifects  of  blast  and  fiagments  on  die  re^jonse  of  die  wall  and  die  icsuhant  synergic 
shear  and  bending  re^Kxise.  The  hysteredcnmterial  model  may  also  be  modified  to  more 


Charge 

Weight 

lb 

TNT 

Stand 

off 

ft 

Cone. 

strengtn 

psi 

Wail 

in 

Wall 

Depdi 

in 

Test 

hfid-^ian 

Deflec. 

in 

Calculated 

D^ect 

in 

Testing 

Ref 

64 

20 

4700 

94 

12.6 

.1 

35 

WHght 

etal,87 

64 

5 

4700 

94 

12.6 

.64 

.76 

Wright 

etal,87 

194 

20 

4700 

94 

12.6 

.9 

.85 

Wri^ 

etal,87 

■ 

- 

5000 

158 

24.75 

2.46 

239 

Ifyde,89 

- 

- 

4700 

65 

12 

.94 

.94 

Coldiatp 

etal,85 

“ 

- 

4700 

65 

12 

.61 

.  .37 

Coldiaip 

etal,85 

- 

- 

4700 

65 

12 

1.94 

1.68 

Hyde,89 

la 

ble  3.6  K 

est  Data 

closely  repnesoit  the  response  charactmstics  of  rdnforced  concrete  A  majcH*  difQculty 
that  must  be  resolved  is  goieralizing  the  iQ'steresis  (or  failure  model)  to  pomit  interaction 
between  shear  and  flexure  failure  surfaces.  This  is  not  a  trivial  matter. 
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CSAFIIR4 


Stochastk  Andysis  A  Hadened  FIki^ 

A1  UnxhKtioii 

Tliiscfafyto  takes  the  scenario  discussed  in  Chqjter  3  and  presents  a  mdfacxtology 
to  predict  die  probcAiility  of  kill  of  a  structural  target  following  a  conventional  weqxxis 
attack.  The  load  and  structure  models  developed  in  Oif^iter  3  are  used,  i^iiere 
{^licaUe,  to  predict  tfae  probability  of  kill  of  a  section  of  the  structure  or  r^kn  of 
interest  (ROI)  given  an  attack.  The  attack  scenario  of  interest  is  an  attack  with  tactical, 
n(xi-precision  wes^pcms,  delivered  in  a  "Stick"  pattern  of  four  or  bonds  per  stick. 
The  target  structure  is  considered  killed  ifany  portion  offoe  outer  shell,  wall  and/or  roof 
is  conpxxnised  to  die  point  that  it  allows  the  infiltration  of  airborne  chemical  or 
biological  agents.  The  iiqHtt  statistics  required  to  assess  die  probalnlity  of  kill  are  given 
in  Table  4.2.  Limit  state  functions  are  ^eloped  fw  both  wes^xsi  delivery  accuracy  as 
well  as  structural  response  using  diese  iiqxit  variables.  Ihe  total  probabili^  of  kill  of 
the  structure  is  the  unicm  of  the  probability  of  kill  fw  each  weE^xxL  The  probability  of 
kill  for  eadi  weqxm  is  determined  by  evaluating  eadi  limit  state  function  dted  in  Table 
4.1.  The  probability  of  kill  for  each  limit  stafo  i,  can  be  written  as: 


* 


(4.1) 


Hoe  the  probability  of  kill,  ,  is  the  jHxxiuct  of  two  probabilities,  die  probability  of 
structural  kill  given  the  wesqxHi  hits  within  diat  kill  regicai  and  die  probability  that 
the  wesqxxi  hits  within  diat  regirm. 
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The  alnictoil  Idfl  Haak  sMe  fiadiocs  qtaHirfied  for  His  work  are  listed  and 


dwcmaed  briefly  in  Table  4.1,  developed  in  detail  in  Sections  42  and  4.3  and  presented 
esq^bdtfy  in  Section  43.1.  Figure  4.1  shows  the  response  regions  dted  in  T£d>ie  4.1. 
A  generic  atnidural  kill  limit  atsle  finctian,  of  a  set  of  vatifl^  x,  may  be 
further  defined  to  hatre  lomf  (fSFi,j^)md  *Ve8istanoe” 

diese  definitions,  die  liniit  atale  finction  is  in  equilibrium  if  LSFlfj^  =  LSF^j(x)  - 
If  fltBt  fits  syt4em  has  been  killed.  The /’mm  is  therefore 

die  probability  that  LSP(x)  ^  0  given  a  wBi|xm  detonation  in  that  kill  regiotL  The 
reqxxise  and  resistanoe  function  portions  of  the  limit  state  finctions  of  Tfole  4.1  are 
developed  in  Sections  42  and  43  le^iectiv^. 

Section  4.4  presents  die  probcftnlity  basis  and  mediods  of  determining  vKapoa 
delivery  accuracy  as  it  has  been  devdoped  by  die  Joint  Technical  Coordinating  Group/ 
Mmitions  Effectiveness  (JTCG/ME)  (JMEH1990).  The  probability  of  a  werqxxi  hitting 


Figure  4.1  Structural  Kill  Regions 
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wilhm  a  ^jedJBc  r^kxi  (Pm)  is  based  cn  these  methods. 


Unit  Stale  FVmctMNis  (LSF) 


LSF-1  Structural  hit,  cratering  and 
perforation  of  die  roof 


LSF-2  Near-field  deUxiation, 
breaching  and  severe  damage  of  the 
near  wall. 


LSF-3  Nfid-field  detCMiation, 
excessive  deflection  of  the  near  wall. 


LSF-4  K^field  det(»iati(xi,  direct 
shear  failure  of  the  near  wall. 


LSF-5  Mid-field  detcmaticxi, 
diagonal  shear  failure  of  the  near 
wall. 


Descriptioo 


Ihe  weapon  hits  die  structiae, 
detonates  and  creates  a  crater  that 
ma^  or  may  not  perforate  die 
stnicture. 


The  models  developed  in  Chapter  3 
will  not  handle  die  severe 
environment  and  structural  reqxxise 
modes  present  after  the  detontfion  of 
a  conventional  weqxm  close  to  the 
structure.  The  empirical 
relationsh^  of  eqpiation  4.4  will  be 
used  to  pr^ct  breading  and  severe 
damage. 


LSF-3  through  LSF-5 

Reqxxise  surfaces  are  developed  for 
the  diree  response  modes  using  the 
riiodels  presented  in  Chapter  3.  The 
reqxxise  sur&ce  equations  are  linked 
to  maximum  allow^le  re^xmse 
equations,  which  are  functions  of  die 
same  input  variables. 


Table  4.1  Limit  State  Function  Definiticxis 
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IVvo  Monte  Carlo  siniilflIioQ  prooedwes  an  developed  in  section  4.5.1  to  assess 
the  probability  of  kill  of  a  hardened  &cility  taking  into  account  weapon  ddiveo' 
uncertainties  as  well  as  uncertainties  in  the  pbyacal  characteristics  of  die  strucbire.  The 
procedures  Kilize  die  developed  limit  state  functions  to  ascertain  die  probability  of  kill. 
The  first  procedure  uses  response  sur&oes  devdoped  to  rcfdace  the  nonlinear  dynanuc 
response  code  (NONLIN)  ofChapler  3,  to  predict  the  response  for  three  levds  of  weqxxi 
INT  equivalent  throw  weig^  500, 100  and  2000  pounds.  The  second  procedure  calls 
the  NOMJN  code  directly  and  is  not  limited  by  weapon  throw  weight  A  discussicm  of 
die  Monte  Carlo  simulation  results  for  several  scenarios  are  presented  in  section  4.6  as 
well  as  a  (xmparison  between  direct  calculation  and  response  sur&ce  sdistitution  reaihs. 

A  procedure  is  discussed  in  section  4.52  to  accomidish  a  probability  of  kill 
analysis  using  reliability  mediods  to  evaluate  die  first  and  secotid  moinents  of  the  random 
variables  (C(xnell,1969;  Ang  and  Comell,1974;  Ang  and  TangU984)  associated  widi 
equaticm  4.1. 

42  Response  Equation  Devetopnient 

Single  re^xxise  equations  are  re  :  for  regions  Rb>l  and  Rh-2  in  Figure  4.1 
while  three  re^x)nse  equations  are  required  for  Rh-3.  These  relationships  are  developed 
in  the  following  three  sections. 

4.2.1  Direct  Ht  (Rb-l) 

)Vhen  a  weqxxi  lands  within  the  footyrint  of  the  target,  diat  structure  is  damaged 


Model  Varirite 

VarMhleNwie 

n 

1-12 

Hystnesis  Shqie  Factor 

AO 

0.(M).1 

N 

A1A2 

^1.0 

m 

N^andN, 

1-10  each 

Nunfaer  of  moment  Oi)  md 
di^iiaoement  (v)  tenns  in  sohrikn. 

TS/tteiBtions 

100-200  etc  . 

TimeStcfn  and  iterations 

StractiR  VaiiaUcs 

Fy 

S0,000>80.000  psi 

Steel  yield  strength 

Fc 

3,000-7,000  psi 

Concrete  compressive  strmgid) 

RHOS 

.0025-.1 

Percent  of  steel  leinforoement 

XL 

96-156  indies 

Wall  height 

Thidaiess 

12-48  indies 

Wall  thickness 

GAMC 

100-200  pcf 

Density  concrete 

Deliveiy  and  Load  Variables 

Zvai 

5-100 

Range  fiom  centnline  of  RCX 

Xval 

0-50  ■ 

Deflection  fixxn  centerline  of  the 

Rd 

Wval 

500-2000 

INT  throw  weight  of  wB^ion 

Table  42  Model,  Structure,  Ddhrery  Conditions  and  Load  VarM)les 
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based  <11  tiikdaiess  of  Ifae  roof  and  the  ttnw  of  the  wea^  Givai,  fir 
exunf^  the  target  as  shown  in  Figure  4.1,  with  the  dimensicns  of  a  typkal  NATO 
scpiadhxn  operadcns  fedlity,  60  feet  on  a  side.  Hie  aim  p(^  will  be  at  the  (xnter, 
feerefire  weapons  landing  widiin  range  and  desQection  offiets  of  30  feet  have  hit  the 
structure.  Ihe  US  Army  Protective  Construction  Design  Manual  (1M5-8SS,1988) 
provides  this  lelationshq)  for  crater  depth: 

Crater  Depth  =  +  (0.33X»^’^  (^-2) 

This  ecpiation  provides  prediction  of  the  crater  depdi  produced  fiom  the  detonation  of  a 
oonventionai  weqxxi  impacting  a  finite  durdmess  of  concrete  oonstnicted  idxwe  an  m 
void,  such  as  a  odlii^  or  roof  slab,  versus  a  slab  constructed  on  grade.  The  TNT  dsow 
weight.  FF,  and  die  depth  of  penetration  ,Dj, ,  prior  to  detonation,  and  the  overall  dqith 
(d)  or  thirdoiess  of  die  concrete  li^,  are  the  key  parameters  in  dus  emperical 
relationship.  Should  this  crater  <kpth  exceed  1/3  of  the  depth  (d/3)  of  the  concrete  h^er, 
that  layer  is  perforated  by  die  confeined  action  of  cratering  of  the  fiunt  face  of  the 
(xmciete  layer  fixxn  in^iact  and  detonation  and  ^»lling  of  the  back  fece  of  the  concrete 
layer  due  to  stress  and  crack  propagation  through  die  k^er.  In  this  stucfy,  only  weqxxis 
that  debxiate  on  in^iact  (noivenetrating  weqxxis  widi  contact  fuzes)  are  considered, 
thoefore  the  term  Dp  of  equaticxi  4.2  will  be  equal  to  zero.  These  fects  lead  to  the 
following  definition  of  the  limit  state  function  for  region  one: 

LSF-1  =  I  -  (0.33)  ♦ 


(4.3) 
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422  Ndii^ficldlktoiMlk»(Rb-2) 

The  models  developed  in  Qvapter  3  provide  an  e?q^cit  oonyiutationai  medxxl  to 
predict  die  response  of  die  structure  in  Figure  4.1.  However,  diete  are  limitations  in  die 
loading  and  re^xxise  portions  of  die  model  ^ch  preclude  its  use  for  very  close-in,  near- 
field  detonations.  A  inoddi^ikiiiiioacporates  die  foUowiiigfeaures  is  required  to  handle 
die  dramatic  loading  and  re^xxise  environment  comes  with  near-field  detonations: 

(1)  esqplidt  diear  deformation,  such  as  a  Tnnoshenko  beam  modd,  (2)  provisions  for 
explicit  fiagment  loading,  (3)  blast  wave  propagation  dirou^  concrete,  (4)  Mure 
surfoce  interaction  between  shear  and  bending  models.  In  lieu  of  such  a  modd,  vdiidi 
is  difficult  to  construct  at  diis  time,  the  National  Defense  Research  Committee,  NE^C, 
(White,1946),  equation,  sipplemented  by  more  recent  test  data  (McVay,1988),  which  was 
uitroduoed  in  Chapter  2,  will  be  used  for  this  r^kxL  The  equation 


w 


V3 


W 


Vi' 


w*C 


(4.4) 


provides  an  enpirical  damage  relationship  between  die  standoff  radius  r,  as  shown  in 
Figure  3.3,  and  a  wall  of  thidmess  t  given  a  cased  wepon  widi  TNT  equivalent  dirow 
weight  w,  and  case  weight  C  The  parameters  a  b  end  c  vary  according  to  dqgree  of 
damage  table  2.3)  to  be  predicted.  This  equation  provides  a  reasonable,  albeit 


enpirical  measure  of  the  response  in  diis  region.  Taking  the  wespon  c^  term, 
out  of  equation  4.4  provides  a  targeting  conservative  estimate  of  die  wall 
repexise  fa-  diese  cases.  In  Bessete,1988,  the  equation 
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r  =  a  ■  a*b*w^  (4-5) 

was  used  to  predict  the  ranger,  at  which  breaching  occurs.  This  equation  is  derived  fixxn 
ecpdon  4.4  after  reinovai  of  the  weafxxi  case  term,  w/(wMI).  Ihe  term  r  is  d^ned  as 
the  breadi  range,  r^  when  a  =  0.18  sodb  =  2.10.  When  these  paraneta  values  are 
substituted  into  equation  4.S,  die  limit  state  function  for  r^on  two  is  produced: 

LSF-2  =  r^-(0.18)(-^5^)  -  (0.378)wV3  i  0  (46) 


4b23  hfd-fidd  Donation  (RH-3) 

An  eiqierimental  design  or  desigp  of  esqieriment  approach  was  used  to 
develq)  respcxise  surfaces  to  replace  the  e?q>licit  use  of  the  con^xiter  model  NONUN 
develqied  in  Onaf^  3.  NONUN  was  used  to  predict  die  deflection  and  shear  re^xMise 
of  a  wall  section  or  region  of  interest  (RC8)  as  defined  in  section  3.4  and  shown  in  Figure 
3.4.  The  DOE  ^iqnoach  aids  in  the  selection  of  qiecific  values  and  pairings  of  the  ii^ 
variables  to  the  axiputer  model.  A  modified  central  ccxiqxsite  design  ((XD)  was  used 
to  represent  eadi  random  input  variable.  The  basic  (XD  for  k  variables  consist  of  a  2^ 
factorial  ctesign,  with  each  facUM*  at  tw«  levels,  the  maximum  and  minimum,  or  -1  and  +1 
(Peters(Mi,1985).  Table  4.2  lists  the  significant  variables  associated  with  die  NONUN 
respcmse  model  and  the  related  loachng  model. 

The  model  parameters  listed  in  Table  42  will  not  be  varied  in  die  stochastic 


analysis.  The  values  for  the  iQ'steretfo  variaUes  were  provided  from  Baixr  and 
Wen,1981,  and  are  intended  to  represeitt  a  typical  generic  l^steredc  material,  such  as 
reinforoed  concrete.  Work  by  Sues  et  al,l^,  extended  the  generic  faysteretic  model  of 
Baber  to  represent  reinforced  concrete  e>q>licidy  and  could  be  incorporated  to  refine  the 
model.  A  basic  modeling  uncertainty  niilti[^  may  be  ifJi^ied  to  the  modd 
representation  as  was  suggested  in  Twisdale,  et  al,198S,  but  will  not  be  included  in  this 
work.  The  effect  of  the  number  of  terms,  N,  and  used  in  die  sdution  was 
investigated  in  Section  3.3  for  a  typical  loading  scenario.  As  disaissed  in  Section  3.3  the 
number  ofterms  used  effects  die  tiiear  prediction  more  than  die  deflection  predictioa  For 
the  purpose  of  reqionsesurfooedeveloixnent  five  Hi  tenns  are  used.  Thisnunfoer 

of  terms  maximizes  the  deflection  predictions  and  returns  shear  predictions  within  95% 
of  dieir  maximunL  Table  4.3  shows  die  variables  sdected  to  mq>  as  random  predictors 
in  the  CCD  alcxig  with  their  associated  data  units  and  settings  This  table  is  taken  fixan 
the  DOE  software  package  RS/Discover  (BBW,1993),  used  in  this  analysis.  RS/Discover 


Name 

Ablxev 

Data  Units 

Settings 

INPUT 

Z_Range 

Z 

Feet 

5  to  50 

X  Deflection 

X 

Feet 

0to50 

Wall  Thickness 

W  T 

Indies 

24  to  36 

Wall  Stroisth 

W  S 

PSI 

3000  to  7000 

Wall_Hei^t 

W  H 

IiKhes 

108  to  168 

OUTPUT 

Max  Deflectirai 

M  D 

Indies 

IHrect  Shear 

S  DIR 

Kips 

Diagonal_Shear 

S  DIAG 

Kips 

Table  4.3  Reqxxise  Surface  Rand(»n  Predictors 
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recommended  28  runs  for  the  (XD  inchiding  ngd^wints.  Response  sur&xs  were 
devdoped  for  three  levels  of  weqxxi  TNT  weigjht,  500  and  1000  aid  2000  pounds.  The 
500  and  1000  pound  figures  oorre^xxid  to  the  INT  ecpiivalent  e^losives  found  in  two 
common  US  general  purpose  bombs  (MK83  and  MKM)  whereas  the  2000  pounds  is  a 
h^podiedcal,  yet  feasible,  weapon  weight  These  levds  were  chosen  purdy  to  devdop 
the  mediodology.  The  28  runs  were  accomplished  at  each  level  using  the  NCMJN 
prooedureL  RS/Discover  performed  a  least  squares  fit  to  the  data. 

The  procedure  used,  in  conjunction  with  RS/Discover,  for  fitting  the  respcmse 
surfiices  was  as  follows: 

(1)  Generate  data  and  feed  iitto  die  RS/Discover*s  multiple  regression  model 
generato-,  MULREG. 

(2)  Generate  a  model  for  the  five  input  variables  and  thdr  re^xxises.  Initially 
cxily  quadratic  terms  woe  included,  however  it  became  qparem  early  that  a  few  terms 
as  high  as  fourth  order  would  be  required.  The  initial  models  therefore  included  the  main 
effects,  the  main  effect  interactions  and  squares,  cubics  and  quarters  of  the  main  effects. 

(3)  The  MULREG  routine  walks  cme  through  a  process  of  refining  the  models. 
A  test  of  the  regression  models  goodness  of  fit  is  based  on  an  analysis  of  die  sums  of  die 
squared  residuals.  The  residual  sum  of  squares  fw  the'respraise  Y,  using  a  least  squares 
regression  technique,  is  given  as: 

-  TiiYt  -  (4.7) 

Here  the  y,  is  the  real  data  point  and  Y,  is  the  model  predictioi 

The  regressi(xi  sum  of  squares  fen-  ai^  model  is  defined  as: 
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»«»  *  ('*•*) 
Ikre  b  is  tile  vector  of  estimated  regression  coefiSdents,  without  die  constant  teem.  Hie 
tennA!'  is  the  /i  xp  matrix  of  predictor  variables,  iwheren  is  tiie  number  of  re^xxises  and 
/7  is  the  number  of  terms  in  the  model.  The  matrix  PCX)'^*  is  PCX)'^  witii  the  row  aid 
column  omitted  for  tiie  constant  term.  The  total  sum  of  squares  is  defined  as: 

sw  -  «  SS«ai. 

The  prqxxtion  of  tiie  variance  or  lack  of  fit  of  tiie  model  vhidi  can  be  attributed  to  tiie 
regressirxi  is  called  R  squared  or  and  is  defined  as: 

'  (“mm  “*««>)  /  “mm 

The  amount  of  variance  eiqilained  or  ac^usted  for  tiie  number  of  d^rees  of  fieedom  is: 

^^ADJ  ~  (^^TOIAL  ~  ^^RKsn^  /  ^^TOIAL 

=  1  -  (1  -  i?2)(n  -  c)/(n  -  p) 

where  (4.11) 

^^TOIAL  ~ 

^^KESID  ~  ^RRXirJ  (W  ~P) 

This  metric  is  affected  the  number  of  terms  in  the  model  that  (km't  contribute 
significantly  to  tte  fit  of  the  model.  When  nm-axitributing  terms  are  eliminated  fixmi 
the  model  this  term  inoeases  towards  ti^  value  of  R^.  Here  the  c  term  is  equal  to  1  if 
there  is  a  ccxistant  in  the  model  and  zero  if  not. 


’■M.  »««• 
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MULREG  also  recommends  iti^xinse  tmisfonnation  to  increase  the  modd  fit 
Reqxxise  variables  may  be  transformed  by  taking  either  their  natural  log,  square  root 
inverse,  or  inverse  squats  root  For  a  model  to  be  acoqsted  on  die  basis  of  goodness  of 
fit  had  to  be  greater  than  0.85. 

(4)  A  final  check  of  die  residuals  was  required  for  full  aooqjtance  of  the  model. 
Ideally  all  residuals  would  fall  within  3  standard  deviations  of  die  fitted  data  If  diis  was 
not  true  die  residuals  were  checked  to  be  normally  distribided  by  applying  either  a  chn 
square  or  >^lke-Slupiro  test  of  noimality.  If  the  residuals  had  a  normal  distribution  die 
model  was  accepted. 

(5)  If  the  models  still  did  not  pass  die  test  die  residuals  were  analyzed  and 
additi(»ial  runs  acccHiqilished  until  the  criteria  was  met  (Cuny,1994). 

Table  4.4  summarizes  the  nine  models  generated  using  RS/Discova.  The  models 


Response 

Regression  Technique 

Numiierirf’ Terns- 

Max  Deflection 

500# 

Least  Squares 

14 

1000# 

14 

2000# 

15 

Direct  Shear 

500# 

Least  Squares 

13 

1000# 

14 

2000# 

10 

Diagonal  Shear 

500# 

Least  Squares 

15 

1000# 

14 

2000# 

11 

Table  4.4  Re^xxise  Suri^  Summary 
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are  grouped  by  response  with  one  model  for  eadi  weapon  load  level  (ie.  500, 1000  and 
2000pounds).  Each  model  is  elalxsrated  on  in  Appendix  C  through  a  series  of  tables  and 
grq^  fixan  RS/Discovo-.  Presented  for  each  model  will  be  a  summaiy  of  the  moctel 
diaracteristics,  a  listing  of  the  model  terms  and  coefficients,  a  listing  of  the  analysis  of 
variation  (ANOVA)  table  and  a  gr^  of  the  model  residuals.  As  an  example.  Table  4.5 
presents  infcxmatim  on  die  500  pound  deflecticai  model  as  it  was  provided  by 
RS/Discover,  with  only  sli^  format  modifications.  Figure  42  shows  a  RS/Discover  plot 
of  the  residuals  fcH*  the  same  model. 

4J  Resistance  Equation  Development 

Deflection  and  shear  design  limit  states  are  defined  per  die  US  Air  Force’s 
Pit>tectiveC(Mistructi(»i  Design  Maiiual(Drake,etal, 1989).  Ftx  die  purpose  of  this  stucfy, 
reaching  the  ipper  design  limits  in  this  manual  will  be  used  to  define  the  limit  or  Mlure 
state  in  these  two  response  modes.  From  a  lethality  stand^int,  the  use  of  the  maximum 
allowable  design  equations  to  define  the  limit  state  still  results  in  a  lediality  ot  targeting 
nonconservative  estimate. 

43.1  Flexural  Resistance  Lindt 

Flexure  limits  are  adopted  based  on  testing  of  shallow  buried,  flat  roofed  structures 
(Kiger  and  Albitton,1980)  subjected  to  conventicMial  explosives.  Fa*  beams  with 
slenderness  ratios  (length  over  thickness,  L/T)  of  five  cr  more,  die  allowable  deflection 
over  length  (v/L)  is  v/L  =  0.10.  For  beams  of  L/T<  5,  die  allowable  v/L  =  0.06.  These 
figures  represoit  what  was  termed  sevoe  damage,  whidi  for  d^  purpose  of  diis  wcxlc 
translate  to  a  conpromise  of  the  wall  as  stated  in  Secdon  4.1. 


SOO  Poaad  DcflMdan  Model  Sumnary 


1  Model  Name: 

MR5S3 

2  Response  Transformation: 

Untransformed 

3  Method: 

Least  Squares 

4  Wei^ts: 

None 

5  Total  Number  of  Cases: 

29 

6  Number  of  Predictors: 

5 

7  Number  of  Unexpanded  Terms 

14 

8  Number  of  Excluded  Cases 

0 

9  Error  Degrees  of  Freedom 

IS 

10  Standard  Error  of  Respon 

6.013722 

1 1  ReUtive  PRESS: 

0.407118 

12  Root  Mean  Squared  Error 

2.984413 

Least  Squares  Coeffkieats 


Term 

Caeff. 

-Std.  £irer 

T-value 

Signif 

1  1 

-4.062295 

17.936065 

-0.23 

0.8239 

2  Z 

0.652205 

0.345298 

1.89 

0.0784 

3  X 

-0.42621 1 

0.088692 

-4.81 

0.0002 

4  W_T 

-2.505735 

0.583184 

-4.30 

0.0006 

5  W_S 

0.007857 

0.002558 

3.07 

0.0078 

6  W_H 

0.181769 

0.106889 

1.70 

0.1097 

7  Z*W_T 

-0.000047 

0.000021 

-2.22 

0.0421 

8  Z*W_H 

•0.002520 

0.001385 

•1.82 

0.0889 

9Z*W 

0.005176 

0.003065 

1.69 

0.1120 

10  X*W 

0.013339 

0.003119 

4.28 

0.0007 

1 1  W_T*W_H 

-0.000041 

0.000015 

•2.68 

0.0170 

12  W_H*W_T 

0.004851 

0.002283 

2.13 

0.0506 

13  Z**2 

-0.003385 

0.002456 

-1.38 

0.1883 

14  W_T**2 

0.022333 

0.009945 

2.25 

0.0402 

Least  Squares  Summary  ANOVA 

Source 

df  Sum  Su 

Mean  Sq. 

F-Ratio 

Sipiif. 

TotaKCoir.) 

28  1012.616 

Regression 

200  891.071 

44.554 

2.93 

0.0609 

Linear 

5  535.077 

107.0*5 

7.04 

0.0083 

Non-linear 

15  355.995 

23.733 

1.56 

0.2667 

Residual  ' 

8  121.544 

15.193 

Lack  of  fit 

7  121.544 

17.363 

Pure  error 

1  0.000 

0.000 

Model  obeys  hierarchy.  The  sum  of  squares  for  linear  terms 
is  computed  assuming  nonlinear  terms  are  first  removed. 


Table  4.5  Characteristics  of  tte  500  Pound  Re^xmse  Surface  Model 
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Figure  42  RS/Discova  Plot  of  tiie  Residuals  of  die  SOO  Pound  Deflection  Modd 
432  Shear  Resistance  Limits 

Shear  limit  states  are  defined  far  two  categnies  of  Mure,  direct  shear  and 
suppcnt  shear.  Two  conpcmaits  add  to  the  shear  cq)acity  of  a  sectitxi,  the  shear  cq>acity 
of  the  cmicrete  and  the  shear  stroigth  [xovided  the  presence  of  shear  reinforconent 
Span  to  depth  ratios  of  <  5  are  considoed  deq?  beams  and  will  be  handled  sqiarately. 

43.2.1  Support/Diagonal  Shear  Resistance  Limits 
433.1.1  Nonnal  Depth  Beams 


The  critical  section  ftH*  determining  shear  failure  is  at  a  distance  equal  to  the 
effective  depth  of  the  member  from  the  siqrprat  or  area  of  ccHiceiitrated  load  (Drake,et 
al,1989).  Fot  members  with  ^)an  to  dqjth  ratios  greater  than  5  die  shear  c^iacity  of 
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die  concrete  is  given  as: 


V,  .  (1.9^'  »  23OOP.0)  M  =s  3  (‘*  •2) 

This  is  die  ixaninal  shear  capadaty  per  the  American  Concrete  Institute  (AO)  Code 
Requirements  for  Reinforced  Concrete  (AO  31&-89).  Ihe  definition  of  is 

p,  =  A  (4.13) 

bd 

v4iere  A,  is  the  area  of  flexural  steel  and  is  the  beam  web  widflL  V  and  M  ate  the 
shear  and  mcanent  at  the  critical  section  and  d  is  the  effective  depth  of  tiie  critical  section. 
The  shear  strengdi  provided  1^  shem*  rdnforcement  is  given  as: 

V  =  ^  (sina  +  cosa)  ^  fiJ^'bd  (4.14) 

sb  - 

The  term  Ay  refers  to  the  cross  secticaial  area  of  the  shear  reinfotcement  witiiin  a  distance 
s.  Combining  the  limits,  the  diagcxial  sl^  limit  state  for  normal  depth  beams  becomes: 

The  total  nominal  ^lear  strength  (V^-  should  be  ccxiservative  as  it  assumes  shear 
reinforcement  vshich  may  not  be  present. 

432.12  Deq[)  Beams 

The  shear  capacity  of  mohbers  with  clear  span-dq}di  (L/d)  ratios  of  less  than  S 
is  finther  txdcen  down  into  L/d  <  2  and  2  ^  L/d  <  S.  The  normal  shear  strength  values 
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are  iniltqdicd  by  an  increase  fictor  to  obtam  tiie  fo  limit  state  equations: 


V,  .  tfjbj 

(4.16) 

> 

H 

+ 

o 

for  2  i  -  <  5 

D 

(4.17) 

A322  Direct  Shear  Resistance  Limits 

Direct  shear  Mure  usually  occurs  near  construction  johts  or  supports  and  are  a 
result  of  hi^  shear  forces  interacting  widi  preexisting  cracks.  IMer  the  dynamic 
loading  resulting  frcxn  conventional  wet^xxrs  detonations,  craddng  of  sections  is  not 
uncommon.  The  direct  ^lear  Mure  occurs  along  cracks  in  the  direction  ofthe  shear  load 
versus  inclined  crack  formations  and  Mure  planes.  Static  tests  of  monollthically  cast 
construction  joints,  for  typical  protective  constriiUion  sections,  has  diown  an  increase  in 
shear  o^ity  ova:  the  starxiard  AQ  equadai  (Karagmian  and  Case,1973)  and  has 
resulted  in  the  ]X‘qx)sed  limit  state: 

=  Q.WjbJi  +  1.4(P  +  i4^/p  i  QJSlfJbji  (418) 

4.4  Weapon  Deliveiy  Accunvty 

Figure  4.3  sho\vs  a  typical  inpact  pattern  fts*  a  stidc  of  four  weqxms  dropped 
from  a  single  aircraft.  The  dipersal  of  the  we^xxis  relative  to  die  center  of  die  stidc, 
\^ch  for  this  figure  is  also  die  aim  point,  is  a  function  of  the  release  conditions  Miiidi 
will  not  be  considered  in  this  wcxk.  Two  errors  are  associated  with  this  stick  inpact 
pattern.  They  are  the  aiming  error  and  the  weqxxi  ballistic  error.  Bodi  are  measured  in 
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tmm  of  mils,  <r  tfic  deviation  m  feet  per  1(X)0  feet  dait  range  fitim  weafxn 
to  die  aim  point  For  exan^le,  the  ballistic  error  may  be  provided  as  a  circular  error 
probable  (CEP)  of  5  mils  and  the  aim  point  error  may  be  given  as  a  CEP  of  20  mils.  The 
slant  range  would  also  be  given,  for  exan^le,  as  6000  feet  Fmn  Ch^iter  2  recall  the 
definition  of  die  CEP  as  being  the  radius  of  a  circle  in  which  50%  of  the  Ixxhbs  would 
land.  The  aim  point  error  standard  deviation  in  terms  of  die  given  CEP  is  defined  as: 
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_  CEPifi)  (4.19) 

d^/Uanm  fw,f  j  1774 

As  an  exan^Ie,  the  aim  point  of  Figiae  4.3  is  shown  as  (0,0).  If  the  slant  range  fixxn  the 


MFC 

mean  oonciete  strength 

SIGPC 

stamtard  <kvuti(»  the  concrete  strength 

MXL 

mean  wall  hei^  or  beam  element  length 

SIGXL 

standafd  deviation  of  the  wall  height 

MD 

mean  depth  wall  or  thickness 

SIGD 

standard  deviatkn  of  the  wall  thickness 

MXBl 

mean  X  coofdintie  bomb  one 

MYBl 

mean  Y  coordinate  of  txxsib  one 

MXB2 

mean  X  coordinate  of  bomb  two 

MYB2 

mean  Y  coordinate  of  bomb  two 

MXB3 

mean  X  coordinate  of  bond)  duee 

MYB3 

mean  Y  coordinate  of  bomb  three 

MXB4 

mean  X  coordinate  of  bomb  four 

MYB4 

mean  Y  coordinate  of  bomb  four 

SIGBOMB 

standard  deviation  for  each  bomb 
based  on  standard  ballistic  error  for 
that  type  of  bomb 

SIGRANGE 

standard  deviation  of  the  aim  point  in  die  range,  or  X 
direction  given  die  aim  point  is  at  0,0 

SIGDEFLECT 

standard  deviation  of  die  aim  point  in  the  deflection,  or  Y 
direction  given  the  aim  point  is  at  0,0 

Table  4.6  SIMTAC  Iipjt  Variables 
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wcapcn  rdease  p(^  to  tfie  aim  pc^  vias  SOOO  feet  and  ttie  CEP  was  20  imls,  fee  aim 


point  standmd  deviation  would  be  ec^  to: 


o— ■  =  o.-.  -  _ =  =  84.93  ft  (420) 

The  aim  point  error  ntay  also  be  provided  in  tenns  of  range  and  deflection  error  probable 
(REP  and  DEP)  terms,  in  feet 

Using  weqxn  1  in^jact  point  (WPN  1)  on  Figure  4.3  as  an  exaraple,  ifs 
detcxiadcxi  location  is  normally  distributed  in  fee  X  and  Y  directions  wife  a  mean  point 
ofin^nct  of(50,-40).  Hie  range  to  feat  induct  point  is  SOOO  feet  and  ballistic  accuracy 
(CEP)  is  5  mils.  The  ballistic  error  standard  deviation  for  werqxxi  one  is  feerefore: 


^wm-i 


(5*5) 

1.1774 


*  2123  ft 


(4.21) 


Fcm-  the  purpose  of  this  study,  means  and  standard  deviations  for  fee  stick  cento' 
or  aim  point  and  the  ballistic  error  of  eadi  weqxxi  will  be  provided  as  straight  iipit 
Current  automated  JMEM  tools  provide  the  stidc  pattern  as  well  as  the  probability 
statistics  given  the  we^xxi,  aircraft,  aira:aft  load  ccxifiguration,  and  the  weqxxi  drop 
sequoKe.  An  interface  between  these  existing  BASIC  programs  and  feis  work  would  be 
easily  acccxi^lished  This  was  not  dcxie  fcx*  fee  actual  weapon  delivery  accuracy  because 
portions  of  the  actual  data  is  classified  SECRET  or  CONFIDENTIAL 
4^  Probabilistic  Analysis  of  Kill 

The  limit  state  functions  defined  in  sections  4.1  through  4.3  are  evaluated  using 
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a  Monte  Cario  anadation  sdieme  to  predict  die  probdnlity  of  kill  of  a  structure  given 
mean  and  standard  deviations  of  the  tqxjt  vaudiles  in  Table  4.3.  An  alternative 
approadi,  not  in^lemented  herein,  is  to  perform  a  fira  order  -  second  moment  (FOSNd^ 
analysis  of  the  limit  state  functions.  A  discussion  of  diis  qiproadi  is  provided  in  section 
4.52 

45.1  Moide  Grio  Simubtion 

A  FORIKAN  program  was  written  to  perform  die  Monte  Carlo  simuladon  of  an 
attack  £^ainst  a  hardened  fecility.  The  source  code  for  the  program  SIMTAC  is  provided 
in  ^ipendix  D.  Two  versions  of  SIMTAC  were  developed.  The  first  version 
(SIMTACl)  uses  the  response  surfiioes  developed  in  this  diqiter  as  replacements  for  the 
ncxilinear  dynamic  analysis  routine  NCMJN,  developed  in  Oiqiter  3.  The  second 
vasi(m(SIMrAC2)interfkxs  directly  widi  the  NCMJN  program.  The  input  parameters 
of  the  program  are  given  in  Table  4.6.  The  following  is  an  overview  of  die  procedure 
SIMTAC  follows. 

Each  run  evaluates  the  kill  probability  fen*  the  four  weqxxis  in  the  sdde.  Sanpling 
for  the  concrete  stroigth  (FC),  wail  thickness  (d)  and  wall  hei^  PCL)  are  acconplished 
initially  and  are  valid  for  an  entire  stidc  san^ling.  That  is,  FC,  d  and  XL  are  held 
omstant  f<x  the  analysis  of  four  different  wetpms  and  their  locations.  Given  the  means 
and  standard  deviations  of  these  values  a  unifram  variate,  randcan  number  generates  is 
used  with  a  Log-Trig  transfomaticHi  process  to  produce  a  sanple  of  a  normally 
distributed  random  variable  (xi  eadi  call.  All  parameters  are  assumed  to  be  normally 
distributed  for  diis  study,  however  die  mediodology  is  not  restricted  to  diis  assumption 
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and  in  reality  mai^  of  diese  parametm  ^  not  be  nocmally  distri^^  The  indhidual 
weapon  locations  are  determined  tty  first  sampling  on  the  distribution  of  each  weq)on 
b^ed  (XI  its  ballistic  enor  properties,  A^ch  are  ii^  variables.  This  retisns  a  locaticxi 
(RANGE1,E£FLECT1)  relative  to  die  stick  center  for  diat  bond).  Subsequently,  die  aim 
point  error  is  applied  with  the  aim  point  standard  deviations  (sigrange  and  sigdeflect)  and 
RANGEIJEFLECTI  as  the  mean  location  coordinates.  The  RANGE^I^FLECT 
cocxdinate  of  the  current  weapon  is  dien  establidied.  This  stqi  is  accomplidied  four 
times  within  (xie  run. 

At  this  point  die  screening  of  the  detonation  locations  b^jns  to  determine  vhicdi 
limit  state  functicxis  are  required  be  evaluated  for  diat  weapon.  The 
RANGEJDEFLECT  cocxxiinate  is  first  checked  to  see  if  it  Ms  outside  die  structural 
damage  zone  identified  in  Figure  4.1.  If  it  does  a  "no  kill"  is  recorded  and  the  next 
we^xxi  evaluation  begins.  If  it  doesn't  the  analysis,  proceeds  with  the 
RANGE,DEFLECT  coordinates  being  scxeoied  to  see  vdiedier  the  weqxxi  hits  die 
structure  or  falls  within  the  breadi  zone.  As  discussed  in  section  4.2.1,  if  die  structure 
is  hit,  LSF-1 


LSF-1  =  I  '  (0.33)(»^^)  s  0  (4.22) 

is  evaluated.  If  poforation  is  achieved  a  "kill"  is  reccxded  and  die  next  weapon 
evaluati(xi  begins. 

If  the  we^n  falls  within  the  Ixeach  range  as  shown  in  Figure  4.1  and  calculated 
by  LSF-2, 


j 
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LSF-2  *  R  ..  -  Rg  iO 

n/V3 

=  Ji^-(0.1SX—)  -  (0378)W^  i  0 


(423) 


a  "kiU"  is  recorded 

If  a  weqxm  does  not  hit  die  stnicdjre  or  M  iwithin  the  breach  zone,  die  analysis 
continues  by  coiverting  the  RANGEJXFLECT  coordinates  into  local  coordinates, 
ZVAL,XVAL  TheZVAL»XVAL  coordinates  are  recpiiivd  as  iipd  to  ^I0^04  in  order 
to  perform  deflecdoi  and  shear  analysis.  Hiis  procedure  determines  the  closest  point  cm 
die  structure  ^di  would  see  the  maximum  loading  fiom  diat  weapoiL  The  comers  are 
(xxisidered  unsuscqjtible  to  deflection  or  diear  Mures  as  diey  are  calculated  Siould  a 
we^xm  fall  outside  a  distance,  in  feet,  (30  -  2.S'''(0  fixxn  the  centerline  of  die  structure, 
the  ZVAL^XVAL  coordinates  are  based  on  the  distance  fiom  the  detonation  to  the  nearest 
face  at  a  point  2.5  times  the  wall  thickness  or  dqith  (d)  fixxn  diejcomer.  This  estimate 
is  ccMiservative  and  is  based  on  the  &ct  diat  the  comers  are  normally  hi^y  reinfcxced 
with  steel  extoidng  a  minimum  of  a  wall  thickness  in  to  die  adjacent  wall. 

Deflection,  diagcsial  diear,  and  direct  shear  analysis  are  performed  using  the 
ZVAUXVAL  detcxiati(xi  locatiors  determined  above.  Should  any  limit  states  be 
^ceeded  the  procedure  junps  ahead  and  sets  die  kill  flag  and  goes  (XI  to  die  next  we^xxi 
or  begins  a  new  run  alter  reccxding  the  kill. 

There  are  two  limit  state  luncticxis  for  the  deflection  based  on  the  span  (XL)  to 
depth  (d)  ratio  of  the  wall.  If  the  deflection  re^xxise  surMe  equaticxis  developed  in 
section  42.3  return  predicticxis  that  exceed  the  limits  stated  in  section  4.3.1  a  "kiH"  is 


recofded.  These  limit  state  fiincdons  are  given  as: 


LSF-3  =  0.10  -  s  0 

XL 


XLId  ^  S 


(424) 


LSF-3  =  0.06  -  i  0  XI4d  <  5  (425) 

JCL 

There  are  dsee  limit  state  functions  for  diagonal  idiear  ^cfa  are  also  based  on  the 
span  to  dqjth  ratio.  Ifthe  diagonal  shear  re^xxisesur&ce  returns  (I^agSiear)  predictions 
Mhicfa  exceed  the  limits  stated  in  section  4.32,  a  "kill''  is  recorded.  These  diagonal  ^lear 
limit  state  functions  are  givoi  as: 


LSF-4  =  11.5  -  DiagSkear  s  0  XLfd  >  5  (4-26) 


LSF-A  =  -  DiagShear  i  0 


(4.27) 


2.  ^<5 

d 


LSF-4  =  S^'bjd  -  DiagShear  i:  0  ^  <2  (4.28) 

^  d 

The  following  single  limit  state  functicm  defines  the  critaia  fa-  a  direct  shear 


(DirShear)  kill. 
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LSF-S  =  0.51  fjbd  (4^) 

After  all  runs  are  ocxv9)leted  the  total  numba*  of  kills  is  divided  die  total 
number  of  we^xxis  delivered  to  give  the  percent  probability  of  kill  for  a  single  attack 
with  the  given  input  characteristics. 

45.2  First  Older- Second  Mnnnd  Method  Discusskm 

This  section  is  broken  down  into  two  parts.  The  first  part  assesses  die  prdiability 
of  kill  for  entire  stick  of  wetpons  and  die  sectxid  part  looks  at  die  individual  weapons 
probability  of  kill. 

45uLl  nobriiility  of  MU  For  a  Stick  of  \¥eap(His 

This  section  discusses  a  procedure  for  assessment  of  die  probability  of  loll  of  a 
structural  target  and  a  given  attadc  scenario  using  a  first  order  second  moment  (FOSNO 
differential  method.  The  total  j^obability  of  kill  firrmi  a  stidc  of  we^xxis  is  the  union  of 
the  probability  of  kill  of  eadi  individual  weqxxi,  as  shown  hoe. 

P[^:i  U  U  /Tj  U  U  =  PiK^)  +  P{K^)  +  +  PiK^ 

-  PiKfiK^)  -  PiKfaCy)  -  Pf^jfi^y 

+  +  pf^TjfWjriAy  (4.30) 

+  PC^inATjlXl  PiKf\Kf\K^l 
-  P(KpKf\Kf\K^ 

A  weqxxi  kill ,  Kf ,  of  a  stick  of  weqxms  is  not  an  indqiendent  event  due  to  the 
tie  of  e^h  weqxin  to  the  aim-point  error  of  the  stick.  The  individual  kill  events  are  also 
not  mutually  exclusive.  These  facts  lead  to  equaticm  4.30  as  a  definiticn  of  {xobability 
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of  kill  for  a  stick  of  4  weffxxis.  EXie  to  the  correlation  ofthe  4  kill  events  tiieprobidnlity 
of  kill  fix)m  ons  weqxxi  is  not  indqiendent  of  another  wet^xxi,  but  is  given  by: 

P(K,  n  K^)  =  PiKJK^)  PiK^  (4.31) 

The  co-variance  matrix  vriiidi  describes  the  correlation  between  die  individual  bomb  kills 
is  unknown  at  this  time. 

4.5^  bidividual  Ihobalality  of  KlI 

Ihe  probability  of  kill  of  an  individual  weqxm  is  the  union  of  die  probability  kill 
by  way  of  each  of  the  S  limit  state  flincticxis  defined  in  Table  4.1.  A  sin^e  functim 
whidi  transiti(Mis  over  the  three  regicxis  is  required  for  evaluation  of  the  probability  of 
kill.  The  limit  state  functions  developed  need  to  be  linked  lay  ig^propriate  transition 
functi(His.  Develqxnent  of  a  model  whidi  acccxnmodates  response  of  the  structure 
eiqilicitly  beginning  at  the  structure  face  would  eliminate  the  transiticxi  problon.  The  fact 
that  determinaticm  of  the  breadi  range  is  a  functicxi  of  the  random  ii^  Variables 
(equation  4.6)  adds  to  the  conpl^ty.  The  correlaticHis  between  the  deflecti(m  and  ^lear 
kill  events  are  also  not  defined  at  this  point  but  are  required  for  complete  evaluaticm.  Fot 
these  reasons  acconplishment  of  a  FOSM  analysis  for  eadi  weqxxi  is  inqpropriate  at 
this  time. 

4.6  Siimdation  Results 

Initially,  the  Mcxite  Carlo,  re^xinse  surface  rqilaconent  program,  SIMTAC,  is 
conpared  to  runs  made  with  the  Monte  Carlo  program,  SIMrAC2,  which  calls  the 
NONUN  p-ogram  directly.  Delivoy  ccMiditions  and  distributions  woe  held  ccxistant  for 


1  TNT 

1  Weight 

1  (pounds) 

Respmise 
Sufiaee  800 
Ribs 

lleqptHBe 

S«&ce4800 

Rms 

Respome 

SiirfKeSOOO 

Runs 

NONUN  1 
Code 

800  Runs  1 

7.6% 

63% 

53% 

6.4%  1 

23% 

Z4% 

23% 

23%  1 

Table  4.7  Reqxmse  Sur&oe  Resuhs  versus  Actual  NCMJN  Code  Results 
all  simulaticxis  and  rqxesent  typical  accuracy  for  inventory  weqx)ns.  Table  4.8  shows 

the  delivery  accuracy  used.  A  conpdson  for  500  and  1000  lb  TNT  weapons  is  given 
in  Table  4.7.  Ihis  conqjarisoi  reveals  diat  as  die  number  of  simulatiocis  increases  die 
ability  of  the  re^xHise  surface  to  rqilicate  NQNUN  results  inpxive.  iliese  results  ^ipear 
reasonable  in  lisjbt  of  die  amplified  central  compo^  desigp  (CCD)  method  used  to 
generate  the  response  surfaces.  It  should  be  noted  that  the  8000  re^xinse  sur&ce  model 
runs  took  qpuximately  two  and  one  half  minutes,  wiiereas  the  800  NO^QJN  runs  took 
^q^iroximately  six  hours  to  conplete.  Ihis  is  a  significant  reducdcxi  in  time  with  an 


Aim  Punt 

Bomb 

2 

(coonO 

Bomb 

3 

(coonQ 

Bomb 

4 

(cooid) 

III! 

Mem 

Sigma 

Range 

(feet) 

Sigma 

Deflect 

(feet) 

I 

70 

70 

15,40 

-2030 

15,-60 

Table  4.8  SimulaticHi  Delivery  Oxididcxis 
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aocxfjtable  reduction  in  accuracy.  A  more  refined  and  robust  response  sur&ce  may  be 
developed  using  a  hig^  fiactional  fiictoriai  experimental  design  metiKxl,  however  such 
a  model  is  not  warranted  for  demonstration  of  die  procedure  put  foidi  in  diis  work.  It 
should  be  noted  diat  no  shear  kills  were  recorded  for  simulations  run  to  die  2000  lb. 
TNT  weight  level.  Several  fiictors  m2Qr  have  coritributed  to  diis  outoxne.  The  following 
are  a  few  possibilities: 

&  A  shear  loll  occurs  primarify  for  detonatioris  at  dose  range  to  the  target  The 
fact  that  the  procedure  established  uses  an  enprically  based  breaching  equation  (section 
42.1)  to  predict  kills  at  most  of  these  same  ranges,  precludes  die  calculation  of  a  diea* 
kill. 

b.  The  NONUN  model  calculates  on  bending  shear  re^xxise  and  as  such  may 
under  calculate  die  total  shear  response. 

c.  F(n-  the  kill  criteria  inclosed,  the  beam  sections  under  investigation,  loaded  by 
the  weqxxis  of  interest,  m£^  not  exceed  the  limit  states. 

In  OTder  to  give  insist  into  the  relative  effect  of  varying  the  characteristics  of  the 
random  structural  variables,  a  series  of  runs  woe  acconqilished  as  shown  in  Table  4.9 
al(mg  with  their  results.  A  base  case  is  shown  in  the  first  row  of  Table  4.9.  Eachrandcxn 
structural  variables  mean  was  set  at  (xie  standard  deviation  frcxn  die  maximum  and 
minimum  values  while  the  otho-  variables  were  maintained  at  dieir  base  case  levels.  A 
preliminary  analysis  of  these  results  reveals  that  die  probability  of  kill  is  most  sensitive 
to  the  thidoiess  of  the  wall  and  roof  and  least  sensitive  to  die  concrete  strengdi  This 
type  of  analysis  is  required  in  order  to  extract  die  characteristics  which  positivdy  effect 
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the  Pk  so  that  greata  efifot  is  spent  atten^idng  to  characterize 
intelligence  resources. 

This  mediodology  will  provide  die  military  decision  makers  die  piobabili^  of  kill 
of  a  single  sortie  or  aircraft  dropping  a  stidc  of  four  weqxins.  The  objective  of  the 
attadc  on  a  specific  target  will  then  be  wd^ied  against  die  nundier  of  aircraft  diat  will 
fly  a  sntie  against  die  facility.  For  exan^le,  if  die  objective  is  to  kill  a  fticility  using 
the  criteria  established  in  this  study,  oiie  sortie  provides  approxiinately  a  15%  probability 
of  kill.  If  this  is  not  acceptable  a  commander  must  make  die  decision  to  send  more  pilots 
and  planes  against  that  target  to  adiieve  die  probability  of  success  he  desires. 


Concrete  Stnmgdi 
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Thickness  ' 
(indies) 
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12 

24 

6 
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— 
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12 

18 

6 
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12 
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6 
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12 

24 

6 

4.0 

&0 

1^^311 
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12 

24 

6 

53 

103 

KM 
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12 

24 

6 

5.0 

9.6  1 
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KM 
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12 

24 

6 

43 

-d 

Table  4.9  Simulation  Results 
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CHAFURS 

CbodiBKHis  and  Recomendilioiis 

5.1  OmchBioiis 

A  procedure  was  developed  diat  returtis  a  probability  of  kill  of  a  hardeiied&cility 
taking  into  account  two  types  of  uncertainties:  we^xxi  delivery  accutacty  and  structural 
characteristics  (x  intelligence  uncertainties.  Ihe  kill  criteria  ate  based  on  die  structural 
re^xxise  of  the  dicility  exceeding  predetemuned  limits  which  rqvesent  the  achievement 
of  the  attadc  objective.  Perfect  knowledge  is  rarely  known  about  the  structural 
diaracteristics  of  a  target  once  a  conflict  is  initiated.  Analysts  tasked  to  perform  pre¬ 
attack  wetqxxis  analysis  and  post-attadc  we^xxis  effectiveness  must  be  able  to  rqxxt  to 
their  siqieriors  realistic  probabilities  of  achieving  die  objective  of  an  airbome  strike  on  a 
target  Curroit  mediodologies  do  incorporate  weipon  delivety  accuratty,  however  duty 
overlook  uncotainties  in  target  struAural  diaracteristics  which.can  make  a  dramatic 
differoice  in  a  ]:xobability  of  kill  prediction. 

The  following  ^^ecific  items  woe  acctxiplished  under  this  effort 

a.  A  nonlinear,  nondimaisional,  oxitinuous  hysteredc  beam  model  was  developed 
to  represent  a  secticxi  of  a  hardened  fadlity  subjected  to  conventional  weqxms  effects. 
The  mocfel  returns  re^xmse  calculaticsis  across  the  height  of  die  section  as  required  to 
fx'ovide  infOTmation  for  determining  the  kill  state  of  a  hardened  target  The 
ncmdimensionalizaticm  allows  for  ease  of  parameter  irput  and  serves  the  stochastic 
analysis  well  where  structural  charactmstics  are  (xxitinuously  changed. 

b.  Existing  etipirical  models  wtudi  generate  conventional  wespon  blast  pressure 
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tune  histories  as  a  function  of  the  INT  eqpuvalent  throw  wng^  were  modi^ 
a  function  of  ^»ce  as  well  as  time.  A  new  model  was  generated  that  i^tumstiie  pressure 
at  any  point  a  wall  section  as  a  function  of  time,  ^reoe  and  an^e  of  incidence.  Ihis 
type  of  load  rqxesentatirxi  was  required  to  feed  toe  continuous  beam  model  referenced 
above.  Ihecoinbinationofthe  beam  model  and  toe  load  model  was  termed  tiieNCMJN 
code. 

c.  Robust  statistical  modds  or  response  surtoces  (RS)  were  derived  from 
NONUN  code  output  calculated  fix»n  typical  ccxnbinations  of  werqxxi  throw  weig^  and 
range  from  the  target,  target  wall  hei^  depth  and  concrete  compressive  strength.  The 
use  of  a  design  of  experiments  (DOE)  or  erqperimental  design  qjproadi  to  toe  RS 
develqjment  ensured  the  RS  closely  reheated  the  ii^  data  across  tiie  parameter  ^race 
ofintoest  The  resuh  was  a  niultidimensional  RS  returned  toe  structural  re^xxise 
given  a  set  of  the  five  parametos  stated  above.  The  RS  rqilaced  the  full  NONUN  code 
within  the  Mcxite  Carlo  simulaticxi  {xogram  whidi  calculated  tire  target  probability  of  kill. 
Use  of  the  RS  replacement  models  allowed  a  simulati(»i  to  be  run  in  less  tiian  1%  of  the 
time  required  to  run  the  simulaticm  with  NONUN.  This  time  saving  is  crucial  to  the  use 
of  this  tool  in  a  dynamic  wartime  envirmmaiL 

d.  The  two  simulaticm  programs  develc^red  (SIMTACl  and  SIMrAC2)  are  the 
first  models  uhich  take  into  account  structural  characteristics  as  random  iiqrut  variables 
in  additicxi  to  the  traditicmal  wetqxm  ^livery  accuracy.  SIMTACl  uses  the  response 
surfaces  which  are  good  cxily  fex*  the  range  of  the  iipit  parameters  frexn  which  toey  were 
derived.  Answers  after  8400  siniulaticxis  are  returned  in  tpproximately  two  and  one  half 
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mimites  on  a  PC  with  an  80486  microprocessor  running  at  33  m^ahertz.  Should  real 
world  circumstances  take  one  out  of  the  valid  iiqxit  ranges  for  the  re^x)nse  sur&ces, 
SIMrAC2,  wiiidi  calls  the  NONUN  code  directly,  may  be  used  without  ii^  restrictions. 
The  time  to  run  ^  simulations  with  SIMrAC2  takes  qjproximately  six  hours.  It  should 
be  noted  that  response  surfaces  may  be  developed  for  any  input  parameters  and  their 
ranges  using  an  ejqjerimoital  design  tpproach.  Use  of  other  models  and  medxxis,  such 
as  die  finite  element  method,  may  also  be  used. 

It  has  been  shown  that  uncertainties  in  die  structural  characteristics  of  a  target  nuy 
significandy  effect  its  response  to  conventional  weqxms  and  die  determination  of  the 
resulting  probability  of  kill  given  an  attadc.  The  use  of  robust  response  sur&ces  to 
replace  conplex  analytic  procedures  allows  for  timely  calculation  of  probabilities  of 
damage  in  spite  of  using  a  simulation  technique  sudi  as  Monte  Carlo.  In  conjunction 
with  a  ^ial-tenp(xal  load  model,  as  presented  herein,  a  total  shear  re^xmse  procedure 
is  required  to  [Hodict  shear  feilures  at  locaticns  of  hi^  load  c(xicaitrati(xis  aw^  fixxn  die 
supports.  A  ccxitinuous  beam  moderrqxesoitation  will  allow  re^xxise  calculation 
against  a  myriad  of  potential  kill  oitoia  whidi  need  mne  dian  a  single  nodal  re^xxise. 
An  exan^le,  which  would  easily  be  acccxnmodated,  would  be  a  kill  criteria  based  on  a 
maximum  length  of  the  wall  secti(xi  exceeding  a  limiting  deflection.  The  methodology 
presented  will  accommodate  studies  to  single  out  the  most  critical  tandcxn  structural 
variables  and  their  critical  ranges  to  allow  the  prq)er  emphasis  to  be  placed  (xi  variable 
significance  in  a  targeting  analysis  and  data  gadiering  exercise. 
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5L2  RBomnDMndirioiis 

The  remaining  sections  of  ^  chapter  discuss  recommendations  for 
improvements,  modifications  and/or  additions  to  die  models  and  procedures  presented. 

5u2.1  StnictmlMMld 

TheNGMJNpcogtam  is  sufiBdently  general  in  its  nondimensioiiBl  form  to  handle 
a  variety  of  cross-sections  and  material  properties  which  are  found  in  die  world-wide 
hardened  fiidlity  cornmunity.  As  stated  in  Chipter  3,  NCMJN  returns  a  crude  estimate 
of  the  total  shear  response  that  could  be  enhanced  by  inclusion  of  a  Timoshoiko  beam 
element  in  die  model.  This  fiict  is  edianced  Ity  the  &ct  diat  the  inodd  never  returned  a 
shear  kill  as  defined  in  Qis^iter  4.  In  a  like  manner,  the  includon  of  fi:agment  effects  may 
be  included  at  the  crude  level  of  die  current  state  of  technology  in  diis  area.  These  two 
additions  should  add  to  the  procedures  ability  to  predict  shear  kiUs  esqilidtly.  As  stated 
in  the  background  infcxmatioi  there  are  numerous  researchers  working  on  the  problon 
of  the  synergistic  effects  of  blast  and  fingmoits  on  die  response  of  the  wall  and  die 
resultant  synergistic  shear  and  balding  respoise.  The  hysteretic  material  model  m^  be 
modified  to  riKxe  closely  represent  the  respcx^  characteristics  of  reinforced  ccsicrete  as 
cited  in  secticxi  4.2.3.  A  majo-  difficulty  diat  must  be  resolved  is  generalizing  the 
hystoesis  (ot  failure  model)  to  pomit  interactim  between  shear  and  flexure  fidhne 
surfaces.  This  is  not  a  trivial  matto’.  These  areas  warrant  the  furtha  investigatitxi  that 
will  add  to  the  robustness  of  the  cunent  predictive  techniques,  however  these  &cts  do  not 
detract  fixnn  the  significance  of  the  mediodology  presented  herein. 

As  discussed  in  section  4.6,  die  procedure  presented  may  be  used  to  hi^g^tt  die 
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critioJ  stmctual  imoertaiiides  most  dramadcally  affect  tfie  probelMlity  of  kill  of  a 
particular  class  of  targets.  Iii  this  ligiit,  a  screening  pcxx^durem^  be  used  to  reccxnmend 
random  variables  vhich  may  be  modeled  fay  a  more  robust  version  of  NONUN  cm*  a 
higher  fidelity  predictive  tool  sudi  as  a  nonlinear  finite  element  code.  Generation  of  die 
respcxise  sur&ces  using  die  revised  randcxn  variables  and  a  more  robust  predictive  code, 
would  greatly  reduce  the  uncertainties  inherent  in  the  process  and  reduce  the  overall 
mcertainties  of  the  analysis  provided  to  air  canq^gi  decision  makers. 

5^  ReqxHise  SurfSace  Replaorant 

One  of  the  objectives  of  this  work  was  to  provide  a  procedure  diat  returned 
probabilities  of  kill  in  as  near  real  time  as  possible.  Hiis  requirement  ccxnes  fixxn  the 
need  to  be  able  to  assess  kill  probabilities  chiring  a  conflict  such  as  I^ERT  STORM, 
in  which  there  is  not  the  time  or  the  inclination  to  wait  for  such  an  assessment  In  the 
same  breath,  the  commander  vdio  wants  an  answor  "now”  also  wants  a  high  degree  of 
ccMifidence  that  the  answer  provided  is  accurate.  The  re^xxise  sur&ce  replacement 
technic{ue  pxivides  this  flexibility  and  ccMifidoice.  Better  reqxxise  sur&ces  may  ahv^ 
be  develcped  as  discussed  in  Gi^ito*  4,  however  the  ^licability  of  this  method  to  the 
analysis  perfcxmed  herein  was  clearly  shown. 

5^  Pmbabilistic  Analysis  Technic|iie 

The  use  of  Monte  Carlo  Simulation,  when  used  in  ccxijunction  with  a  reqxxise 
surface  rqilacemait  models  ran  fairly  c]uickly  (xi  a  DOS  based  personal  ccmputer  with 
a  80486  microprocessor  running  at  33  megahotz.  Eighty-four  hundred  simulaticxis  would 
run  in  qiproximately  two  and  one  half  minutes.  In  die  absence  of  covariance  matrices 
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for  ddivery  axuracy  variables  as  well  as  die  correlation  between  the  diear  and  deflection 
response,  an  esqilicit  integration,  first  (xder,  second  nxxnent  (FOSM)  method,  or  other 


derivative  mediod,  can  not  be  acconplished. 
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Equations  3.9  and  3.10  of  section  3.4  entitled,  Stnictute  Model,  are 
nondimensionalized  with  KspeOt  to  the  ultimate  bending  moment,  M,,  and  the 
curvature,  <|>i ,  vdiere  the  initial  tangent  curve  intercqjts  tiie  horizontal  ultimate  mnnent 
curve.  As  shown  in  Figure  3.12,  is  on^  tiie  ultimate  moment  ^len  Ao  ==  0.  For 
tills  model  Aq  =  .05.  The  ultimate  curvature  and  nondimensional  moment  are  defined 
as: 


The  ultimate  moment  is  defined  as: 


K 


*.'(#-  iiHr/4  -  f )  ♦ 


The  nondimensicxial  length  and  di^lacement  are: 


£ 

L 


¥  » 


The  nondimensional  time  is  defined  as: 


1  ■ 


>\here 


The  ncndinienaGnal  load  and  load  per  init  lengdi  are: 


/• 


tst 

K 


Incorporating  diese  conversions  into  equations  3.7  aid  3.8  results  in  equations  3.9  and 


3.10  of  the  same  section. 
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AiyendhB 


Noniiiinr  Load  and  Stradnri  RespcMiBe  Modd  Cbde 


PROCiRAM  RESPCONG 

<XCXXXXXX(XCXXXX(XCXXXX(XCXXXXCCCCCXXXXXXXXXXD^^ 

C  This  program  will  goierate  die  one  way  slab  or  beam  response  cf  a 
C  given  reinfixced  concrete  protective  constructed  fittility  subjected  to 
C  blast  loadiiig  from  a  delink  conventional  weapon  at  an  infad  locatka 
C  The  loading  oiviroranent  is  generated  from  a  modified,  empirical  code 
C  based  on  the  code  OOMWEP  provided  by  die  Aniiy  Waterways  Eaqxriment 
C  Station  (WES).  The  kiading  emnionmeit  b  fed  inlo  a  fayalcntic, 

C  condinjous  bonn  model  whi^  sdves  fcr  deflections,  moments,  shears 
C  and  curvatures  using  a  weighted  readual  scdution  method.  Response 
C  statistics  mt^  be  generated  at  a  uso*  inputted  matiber  of  prants  along 
C  the  beam  element 

C  The  main  loading  subroutine  b  entitled  LQADCDC  while  the  main  oon- 
C  tinuous  beam  model  subroutine  is  entided  NCNUN.  NCMJN  calb  an 
C  IMSL  roiaine  DGEAR  which  solves  a  set  of  three  first  order  differential 
C  equations  at  each  time  step. 

C 

cccxxcxxxccccoxxtxccxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx: 

c 

C  Variables 
C 

C  NV  •  The  number  of  displacement  terms  to  use  in  the  solution 
C  NU  -  The  number  of  moment  terms  to  use  in  die  solution- 
C  ITER  -  The  number  of  iterations  that  the  solution  will  be  carried  to. 

C  XVAL  •  The  coordinate  left  or  right  of  the  centerline  of  the  region  of 
C  interest,  which  is  "O"  at  the  centerline  where  the  detontation 

C  is  located. 

C  YVAL  -  The  vertical  coordinate  measured  up  the  region  of  interest, 

C  which  is  "0"  at  the  ground  elevation. 

C  ZVAL-  The  distance  perpendicular  to  the  region  of  interest  where  die 
C  detonation  is  located. 

C  RJ  -  The  coordinate  up  the  region  of  intoest  where  the  response  will 
C  be  calculated. 

C  LOAD  -  The  nondimensionalization  fector  which  multiplies  the  load 
C  terms. 

C  TIME  -  The  nondimensionalization  factor  which  multiplies  the  time  toms. 

C  AO-  The  post-yield  to  pre-yield  moment-curvature  ratio  which  controls 
C  the  degree  of  nonlinearity  that  the  hysteretic  model  will  exhibit 
C  AI,  A2,  P  -  Control  the  shape  of  the  hysteretic  curve. 

C  NUM  -  Run  reference  nurtiber  which  is  inserted  into  the  name  of  the  output 
C  file  (ie.  NLCnum"). 

C  WVAL  -  The  TNT  throw  weight  of  the  chai^ge. 

C  XL-  The  height  of  the  wall  or  lengdi  of  the  beam  element 
C  TS  -  The  time  step  at  which  returns  solutions. 

C 


o«  o  nn  oo  o  n  on  on 
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CHARACTER  NUM*3 
INTCGERNVJJUJTER 
REAL  XVAL,YVAUWVAlJUajOAD,TTME 
C 

COMMON  m  XVA1^AUWVAL4TER,TSJU JJlMnja3J=C^ 
COMMON  /CONV/  LOAD,11ME 

COMMON  /PARAMS/  PRV(70XAV(70XTAV(70XTOV(70)^YDATA(70) 
COMMON  /NONL/AOA1>^24‘J^VJ^IUJ1JXXNV,SCONV 
EXTERNAL  IVPAG 
C 

kk>=100 
call  inpul(kk) 

OPEN(UNrr=17JlL&>#IljOT//ta^ 

OI>B4(UNn^l6JlLEF'SUMM//ia^ 

do  lOl^kJI 
CAll  INPUItl) 


CALL  LOADCDC(XVAL^AL,WVAL) 

CALLNONLIN(XL) 

CONTINUE 

END 


CXXCCCCCOCCCOCCOCCOCCCCXXXXOOCOOOOOOOCOCOOOOOOOOOOOOCOCOCOOOCCC 
SUBROUTINE  NONlJN(xl) 


INTEGER  P,N>!ETHINlDEXJ>MC(35XIER^J4VJ^JTER,TOTE»(iIRP 

/  RES 

CHARACTER  NUM*3JTJIMl*44=LNM2M 


REAL  X  WY(35XWK(596XX,TOUXENDAA1,A2^0,TS4VXSPJ40M, 

/  CURV,SHTlP,YPRIME(3SXYM(10XRJ,tinaxs3iiax,dtnax,linKS,lina}^ 
/dsx^ 

COMMON  /PARAMS/  PRV(70XAV(70XTAV(70XTT)V(70XXUYDATA(70) 
COMMON  /NONL/AOwAl  A24»JW3^41,DCONV,SCONV 
COMMON  /IN/  XVAl^AUWVALJTER,TSJU3^UMjn4U)TOT 
EXTERNAL  FCNKNJ JX5EAR 
oPEN(UNIT=12,nLE='CHECK!//CHARaW8y/CHAR(KK+48)) 

PI  =  3.14159 
5  CONTINUE 
C 

A0=.05 

Al=5 

A2=.5 

P=5 


c 

c  WRnEdV)^-';^ 
c  WRnE(12,*)’NU  =  ’JJU 
c  WRnE(12,*)  TIERATIONS; HER 

c  WRnE(12,*)  TIME  STEP: TS 

c  WRnE(12,*)  '  A0UV1>K2: 

c  WRnE(12,*)  -N.-JP 

c  WRnE(12,*)  •XVAUZVAUWVAL’,  XVAUZVAUWVAL 

c  WRTIECV) 

C 

IT  =  nER-l 
TaiE  =  2*NV  +  NU 
Q  =  2*NV 
C 

DO  10 1»1,TOIE 
Y(I)  =  0 

10  CXDNTTNUE 

C 

C  Y(1-NV)  =  GAMMA(K) 

C  Y(NV-2NV)  =  ALPHA(K) 

C  Y(2NV-2NV-fMJ)  =  BETA(L) 

C 

C  DO  15  1=1JW^ 

C  YaH2/(a*PD*(X2-Xl)))*(-COSa*PI*X2)  +  CXDSa*PI*Xl))*Xl 
C15  CXDNllNUE 
C 
C 

N=1UIE 
METH=1 
MITER  =  0 
X=0 

TOL=001 
H  =00001 
INDEX  =  1 
RP  =  TS*nER 
C 

C  Diagonal  shear  calucalte  at  the  efiective  depth  (d-3)  of  the 
C  concrete  section  (nondimensionalized)  away  from  the  support. 
PRINT*,  D^ 

DSX=(D-3y(XL*12) 

Z=TS 

smax=0 

dsniax=0 

dmax=0 

tniaxs=0 

liiiaxs=0 

tmaxcN) 

ImaxcH) 

DO20KKK=  14TER 
0*0 

XEND=Z 


CALL  IXjEAR(N.FCNJO«XHY.XENIXTOL>lElHJ^4r^^ 

\  JtER) 

110  FORMAT(F12.5) 

DO30KKKK»l.ll 

DISIH) 

MOI^ 

SHm 

CURV=0 

diagshm) 

dii^sh2»0 

C 

DO  40 1=1;^ 

DISM)ISPfSINa*H*RJ)*Y(NV-rt) 

CmV<mV-(a*PI)**2)*SINa*H*RJ)^(NV+I) 

40  OGI^nNUE 
C 

DO  50 1=1  J4U 

SH=StHI*H)*OOSa*P!*RJ)*Y((^ 

50  CONIINUE 

if  (disp  .gt  dmax)  then 
dniax=disp 
tmaxd»7 
linaxd=g 
ENCXF 

IF  (ABS(M0M)  .GT.  .9999)  WEH 
D055I=1J4U 
print*,  mom 

Y(0<-I)=(.9999/ABS(M0M))*Y(QH) 

55  CONTINUE 
MOM=0 
SH=0 
DO  60  1=1 

MOM=MOIVffSINa*PI*RJ)*Y(OH) 

SH=SH+(l*PI)*OOSa*PI*RJ)*Y(OfI) 

60  CONTINUE 

PRINT*,'MOMENT  AND  SHEAR  HAVE  BEEN  OORRECTEDI'jacMKxn 
ENDIF 
C 

if  (abs(sh).gt.  abs(smax))  then 
smax=A^sh) 
tmaxs^ 
lmaxs=^ 
endif 
c 

rj=rj+.l 
30  continue 

do  65  i=l,NU 

DiagshI=diagdil-Ki*P>)*co9(I*pi*dsx)*y(qH) 

diagsh2=dia^t2+(i*^)*cos(I*^*(l-da))S(Q^0 
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coiituiue 

if  (riwCduigjdil)  .gt.  absCdsniBx))  then 


endif 

if  (abs(diaG^)  .gL  abs(dsniBx))  then 

dsniMF  ahfl(dMgsh2) 

endif 

ZfZ+TS 

CX}NTTNUE 

DMAX_ODMAX*DOONV 

SMAX_O(SMAX*SC0NVyi000 

dsn)ax_cKdsm8x*sconvyiOOO 

Fdnnt(D^SJ^2xJSJ,2XJF83^116^12.6^123^ 

/  F72;2XF7:j^l2.6^XJFlld) 

write(16330)  KKK^AL^AUWVAI434’C,XU3MAX,SMAXJDMAX_QSM^^ 
FORMAT(8F82) 
xxl=od*12 

WRnE( 17,400)  zval,xvtd,d,fc,xxU^dAXjC,^4AX_c43Siiiax  c 
RETURN 


(XCXXCXXCCCCCCtXCCCXXXXXXCXXXXXXCCXXXXXXXCXXXXXXXXCCXXX^ 
SUBROUTINE  Fa'J(HX,Y.YI»RIME) 

OC  The  fdlowing  is  an  exoeipt  fitm  ^  IX£AR  convnents  oonoeming  the  sub' 

CX^  routmeFCN. 

COXXXXXCXXCCXXCCCXXXXOCOCCXXXXXXCCXXXXXXCXXXXXXXXXXXOOOC^^ 

C  FCN  -NAME  OF  SUBROUTINE  FOR  EVALUATING 

FUNCTIONS.  DGEA0180 

C  (INPUT)  DGEA0190 

C  THE  SUBROUTINE  ITSELF  MUST  ALSO  BE  PROVIDED  DGEA0200 

C  BY  THE  USER  AND  IT  SHOULD  BE  OF  THE  DGEA02I0 

C  FOLLOWING  FORM  DC£AQ220 

C  SUBROUTINE  FCN  (N,X,Y,YPRIME)  DGEA0230 

C  REAL  X,Y(N),YPRIME(N)  DGEA0240 

C  DGEAQ2S0 

C  DGEAQ260 

C  DGEAQ270 

C  FCN  SHOULD  EVALUATE  YPRIME(IX..,YPRIME(N)DGEA0280 

C  GIVEN  N,X,  AND  Y(1X...,Y(N).  YPRIMEa)  DGEA0290 

C  IS  THE  FIRST  DERIVATIVE  OF  YG)  WITH  DGEA0300 

C  RESPECT  TO  X  DGEA0310 

C  FCN  MUST  APPEAR  IN  AN  EXTERNAL  STATEMENT  IN  DGEA0320 

C  THE  CALLING  PROGRAM  AND  N,X,Y(IX..,Y(N)DGEA0330 

C  MUST  NOT  BE  ALTERED  BY  FCN.  DGEA0340 

CC 

CCCCCCCXXXXXCCXXXXXOCCCXXOOCXXXCOCCCCCCOCXXXXCOCtXXXXXXXCCCCCOCC 
INTEGER  HPJSIV WW,XT4>T,XTG,CiXTB4^ 

C 

REAL  FUNl  JUN2JUN3,A0A1.A2.Y(N),YFRIME(NXH,CSG(I00XCSB(100) 
/,CS(500XTOTALl(21),TOTAL2(21XPARTl(2IXPART2(2IXCJ>I,TOn(2IX 


/I0r2(21),XJjOAIXlO)J,TOrrAU69XSU^ 

/^jOADUUME 
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C 

COMMON  /NONU  A0wM,A2J»J4VJJUJ1 
CXDMMON  /CONV/  LOADUUME 

CXDMMON  /  PARAMS/  PRV(70XAV(70XTAV(70XTDV(70XXL,YDATA(70) 

C 

T1  =  'IlME'*(TAV(e9)  +  T0V(69)) 

Q  =  2*NV 

(XCXXXXXXXXXO(XCXXXXXXCXXXXCXXCCCCOCCCCXXXXXXXXXXXXXXXXXXXXXXXXX: 
C  The  sidnoutiiie  UMTIS  takes  die  cuntnt  fonn  oi the  finctioas  FUNl,  FIJN2 
C  and  FUN3  and  finds  dieir  points  of  inflectkn  across  the  beam  and 
CALL  IJMrTS(Y.CSG.CSB,XmXTG) 

CALL  SORlXCSG.CSB.Xnj,XIB,CS^ 

C 

IX)200K=1JW 
IF  (X  .GE  Tl)  THEN 
LOAD(K)  =  0 
ELSE 

C  WRITE(12,101)  FUN4(K4XF(I,XXKLX 
002101  =  1,69 

TOTALa)  =  FUN4(KJ)*F(I,X) 

210  CONTINUE 
SUMl  =0 
SUM2  =  0 
D0220W,50 

IF  ((J  .EQ.  1)  .OR.  (J  .BQ.  50))  THEN 
SUMl  =  SUMl  +  .5*T0TAL(J) 

ELSE 

SUMl  =  SUMl  +  T0TAL(J) 

ENDIF 

220  CONTINUE 

CCC  MULTIPLY  BY  THE  INCREMENT  "DELTA  X"  =  .0102043 
SUMl  =  SUMl  *.0102043 
C 

DO  230  J  =  51,69 

IF  ((J  .EQ.  51)  .OR.  (J  .EQ.  69))  THEN 
SUM2  =  SUM2  +  .5*TOTAL(J) 

ELSE 

SUM2  =  SUM2  +  TDTAL(J) 

ENDIF 

230  CONTINUE 

SUM2  =  SUM2*.0263158 
C 

LOAD(K)  =  SUMl  +  SUM2 
C  WRI're(12,*)  'KLOAIXK)  = ',  KLOAD(K) 

ENDIF 

200  CONTINUE 
C 
C 
C 


c 

IF(NV.GT.NU)THEN 
DO  10  K-1 J4U 

YPRIME(K)=-A()^(K*PI)**4)*Y(NV+K>Hl-A0)*((K*PI)**2)*Y(O^) 
/  +2*LOAD(K) 

10  CONIINUE 

DO  15  K=WflJW 

YPRIME1(K>=.A0*((K*PI)**4)*Y(NV4K)  +  2*IjOAD(K) 

15  CX3NIINUE 
ELSE 

DO20K=lJiV 

YPRIME(K)=^A0*((K*PI)**4)*Y(NV+K>Kl-A0)*((K*PI)**2)*Y(OfK) 
/  +  2*LOAD(K) 

20  CXKnNDE 
ENOF 
C 
C 

DO30K=l^ 

D“NV+K 

YPRIME(L)  =  Y(K) 

30  CONTINUE 

C 

C 

DO40K=l>IU 

DO50I=OT.l 

FIM) 

D060  0CS(I).CSa+lX.05 
PT*PT+1 

lF(P.GE2)-mEN 

TOTl(PT>=FUN  1(K,C)*ABS(FUN2(C,Y))*(ABS(FUN3(CY))**(P-1))* 
/  FUN3(C,Y) 

ELSE 

T0T1(PT>=FUN1(1CQ*ABS(RJN2(C,Y))*FUN3(C,Y) 

ENDIF 

TXDT2(PT>FUNI(ICQ*FUN2(C,Y)*(ABS(FUN3(C,Y))**P) 

60  CONTINUE 
PARTIOH) 

PART2a)=0 
DO  70  K1=1,PT-1 

IF  ((KI.EQ.1)  .OR.  (KI.EQ.(PT-1)))  TH^ 
PARTia)=PARTIflh.5*TCITl(KI) 

PART2(l)=PART2(I)+.5*TOT2(KI) 

ELSE 

PARTia)=PARTia^TCITl(KI) 

PART2(I>=PART2(I>fTOT2(KI) 

ENDIF 

70  CONTINUE 
H  =  .05 

PARTia)  =  H*PARTl(I) 

PART2(D  =  H*PART2a) 

50  CONITNUE 
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T0TALl(K)-0 

T0TAL2(K)«0 

DO80I=OT-l 

T0TAL1(K)=TCITAL1(K)+PART1(I) 

T0TA12(K>-T0TAL2(K)+PART2(I) 

80  CONITNUE 
40  C(»mNUE 
C 
C 

IF  (NU  .GT.  NV)  THEN 
DO90K-1J4V 
J=(2*NV>fK 

YPRIME(J>=K(K*PI)**2)*Y(K)  -2*AI*TCITAL1(K) 

/-2*A2*T0TAL2(K) 

90  CONIINUE 

DO  100  K=^IV+1^ 

J=(2*NV>fK 

YPRIME(J)=-2*A1  *TOTAL1(K>2*A2*T0TAL2(K) 

100  CONTINUE 
FI<CT? 

DO110K=lJ4U 

J=<2•NV>^-K 

YPRIME(J)=K(K*PI)**2)*Y(K)  -  2*A1*TCITAL1(K) 

/  -  2*A2*TOTAL2(K) 

110  CONDNUE 
B4IXF 
C 

RETURN 

END 

CCCCCCCCCCCCC(X(XC(XCCCCXXXCCOCXXXCCCCXXXXCCOOCCCCCCCCOOOOCCCCCC 
REAL  FUNCTION  Fa,T) 

C 

C  This  function  takes  the  parameters  computed  in  the  subroutine 
C  LOAIXDC  and  returns  the  reflected  blast  pressiBC  at  a  given  point 
C  for  agiven  time  after  blasrt  arrival  at  the  \^I.  This  function 
C  represents  a  spatial  Friedlander  equation  \^hich  returns  die  reflec- 
C  ted  pressure  at  point  "I"  and  time  "T"  using  die  equation; 

C 

C  Pa,T)  =  PRVO)  •  H[T  -  TAV(I)]  *  [1  -  ((t  -  TAV(I)yTOV(I))] 

C  *  exp  {  -AVO)  *  ((T  -  TAV(1^V(I)) 

C 

C  Where  the  term  H[T  -  TAV(1)]  represents  a  Heaviside  function  which 
C  prevents  the  load  from  being  applied  prior  to  its  arrival  time. 

C 

C  The  terms  LOAD  and  TIME  are  nondimensionalization  fectors. 

C 

CCCCCCCOCOCCOCOOOOCCOCOOOOCCXXXXXXCCCCtXCCCCCCOCCCCOCOCCCCOOOC^ 

rrrxrrrrrrrrr 

REAL  Y.T4>1  J>2J>3J>4,WGHr4jOAD,TlME 
INIIUNSICEXP 
C 


CMMMON  /PARAMS/  PRVCTOXAVCTOXTAVCTOXTOVCTOXXUYDATACTO) 
(COMMON  /CONV/  LOAIXTIME 
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C 

PI  =  PRV(I)*1jOAD 
P2  =  TAV(I)*T1ME 
P3  =  T0V(I)*T1ME 
P4  =  AV(I) 

C 

P5  =  P2  +  P3 

IF  (T  XT.  P2  .OR.  T  .GT.  P5)  THEN 
F  =  0 

F  =  PI  *  (l-((T-P2)/P3))*EXP(-P4*((T-P2yP3)) 

ENESF 

C 

C 

RETURN 

END 

<X(XCCCCXXXXXOCmX(XCXXXXXXXXXCXXXXXXCCOCCCXXXXCCC(XCCC^ 
SUBROUTINE  IJMITS(Y,CSG,CSB,X1B,XTC) 

C 

C  This  subroutine  takes  the  fuctions  FUN2  and  FUN3  and  detnniines  their 
C  inflection  points  with  the  cunent  set  of  coeficients  T”.  It  takes  the 
C  cunent  set  of  "Y"  coeficients,  calculates  the  inflections  points  and 
C  stores  them  in  the  array  CSG  for  FUN2  and  CSB  for  FUN3.  The  number  of 
C  inflection  points  is  recorded  as  int^en  XTQ  and  XIB  for  FUN2  and  FUN3 
C  respectively.  Ihese  inflection  points  are  used  to  define  the  limits 
C  for  the  numerical  integration  of  the  equations  which  contain  FUN2  and 
C  FUN3  to  preclude  integration  over  a  range  which  crosses  die  x  axis. 

C 

C  Variables 
C 

C  XTG  -  Nunfoer  of  inflection  points  in  the  function  FUN2. 

C  XTB  -  Number  of  inflection  points  in  the  function  FUN3. 

C  CSG  -  Array  of  FUN2  inflection  point  locations. 

C  CSB  -  Array  of  FUN3  inflection  point  locations. 

C  L  -  The  length  of  the  beam  element 

C  Z  -  Temporary  array  of  terms  in  calculating  current  value  of  either 
C  FUN2  or  FUN3. 

C  TSZ-  Term  which  determines  the  existance  of  an  inflection  point  by 
C  dividing  the  current  function  value  ty  the  previous. 

C 

cc<xcccccccccccccccccx:ccccccc(x(xc(x(x<x(x(x0cxxaxx(xc(x(x(xccx: 

ccccccccccc 

INTEGER  XIB,XTGJS[U3W,(iPX 
REAL  Y(35XA0^1>V2,CSB(100XCSG(100),Z(ll),ZCJ»I 
C 

COMMON  /NONI7  A0,A1,A2J>J^J4UJ»I 

Q=2*NV 

XIG=0 

XTB=0 
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DOS  1=1,11 

ZOH) 

5  CXKIINUE 
L=0 
RJ=0 

DO10J=l,ll 

L=L+1 

DO20I=14W 

(X(X(XCCCCCCCCCOOOC(X(XCCXXXXXOCCaXXXXXXXXXXXXXXXXXXXXX^^ 

CC  ThisisFUlSC. 

caxcxxxxxxxxxxxxxcaxccxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxy^ 

Z(L)=Z(LHa*PI)**2)*SIN(I*PI*IU)*Ya) 

20  CXKnNUE 
RNRH-.l 
10  CX3NT1NUE 

c 

J=0 

DO  so  RJ=1,10 
J=J+1 

IF  (ABS(Z(J))  .LT.  .00001)  GO  TO  SO 
Z1C=Z(J)/Z(J+1) 

IF  (ZC  .LT.  0)  THEN 
XTOXTGfl 
CSG(XTC)  =  (RJ/10)-.0S 
ENEHF 

SO  (DOI^IINUE 
DOSS  1=1,11 

2oy^ 

55  CONTINUE 
L=0 
RK=0 

DO  30  K=  1,11 
L=L+1 

DO  40  1=1, NU 

CCCCCCC(XCCCCCCCCCCCC(XCXXXXXXXXCCCXXC(XCXXXX:(XCCC(XCCCCCCCCC 
CC  ThisisFUN3. 

CCCCCCCCCCCCCCCCCCCCCCO(XOCOCOOCCCCCXXCCCCCCCCCCOCCOOOCCCCCCC 

Z(L)=Z(L>4-SINa*PI*RK)*Y(OfI) 

40  CONTINUE 
RK-RK+.l 
30  CONTINUE 
C 

K=0 

DO60RK=l,10 

K=K+1 

IF  (ABS(Z(K))  .LT.  .00001)  GO  TO  60 
ZC=Z(K)/Z(K+1) 

IF  (ZC  .LT.  0)'raEN 
X1B=XTB+1 
CSB(XIB)  =(RK/10)-.05 
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ENIXF 

60  CX3NT1NUE 
C 

REIURN 

END 

CXXXCCCXXXXXXCXXXCCCCXXCCCXXXXXXXXXXXXXXCXXXXXXXXXCCCCCCC^ 
SUBROUTINE  SORT)[CSG,CSB,XrG,XTB.CS^ 

C 

C  S(^T  takes  the  sets  of  inflection  points  generated  in  LIMITS  md 
C  int^rates  diem  into  a  single,  sorted  array  CS  of  size  XT. 

C 

<XCXXXXXXXXXXCXXC(XC(X(X(XCCXXXXXXXC(XCXXCXXC(XCCCXXCXX^^ 

cccxxxxcccc 

INTEGER  XrGOT.WaWJ4U,XA 
REALCSG(100XCSB(100)MCS(200XA0,A1,A2 
C 

XT=XrG+XTB+2 
CS(1)=0 
CS(XT)=1 
DO  10  I=1,XTG 
CSa+1)  =  CSG(I) 

10  CONTINUE 
DO20I=lOT 
CS(XTrHl+I>=CSBa) 

20  CONTINUE 
DO70l»lXr-l 
DO60J=I+l,Xr 

IF  ((CS(I>CS(J))  .LE  0)  GO  TO  60 
M=CSa) 

csa)=cs(j) 

CS(J)=M 

60  CONTINUE 
70  CONTINUE 
C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCOOCCCCCCCCCCCOCOCOCC 
FUNCTION  FUN1(K,C) 

C 

INTEGER  K,! 

C 

REALPLC 

C 

COMMON  /NONIV  A0,A1,A24*,NV,NU4>I 
C 

FUNl  =SIN(K*PI*C) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCOCCCCCCCCCCCCCCC 
FUNCTION  FUN2(C,Y) 

C 
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INTEGER  NV 
C 

REALCJ>1,Y(35) 

C 

COMMON  /NONL/  A0A1,A2J»J^^J>I 
C 

FIJN2-0 
DO  10 1-I JW 

FUN2  -  FUN2  -  (a*P0**2)*SINa*PI*Q*Y(I) 

10  CONI1NUE 
C 

RETURN 

END 

CXXXXXCXXXXXOOCCOCXXOOCXXOCCCOO(XCXXXXXXXOOCCCCOCCCXXXXXXXCCCCOC 
FUNCTION  FUN3(C;Y) 

C 

INTEGER  NUJW,Q 
C 

REALPI,Y(35),C 

C 

COMMON  /NONU  A0^1>^2J>JSIVJ^UJ»I 
C 

FUN3  =  0 
Q  =  2*NV 
DO  10 1=1  J4U 

FUN3  *  FUNS  +  SIN(I*PI*C)*Y(QH) 

10  CONITNUE 
C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCOCCOOCOOCCCCCCCOCOOCCIXCOCCCCCCCCCCOOCC 
FUNCTION  FUN4(K,1) 

C 

INTEGER  K,I 
REALP1,C 

COMMON  /NONU  A0^1^J>4WJ^,PI 

COMMON  /PARAMS/  PRV(70)AV(70XTAV(70),TOV(70)^YDATA(70) 

C 

C  =  YDATA(1) 

FUN4  =  SIN(K*P1*C) 

C 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SUBROUTINE  FCNJ(N^Y,PD) 

C 

INTEGER  N 
REAL  Y(N),PD(N,N)^ 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCCCOCCCCCCCCCCOCCCOOCCOCC 
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SUBROITIINE  lX)ADCIXXXZ,1Nn 
COXXXXXXXXXCXOCCXXXXXX:  4  JANUARY  91 

cxxxxxxxxxoccxxxxxxxxxxxxxxxxocxx: 

c 

C  This  subroudne  is  a  derivhive  of  source  code  taken  from  the  program 
C  CXDNWEP,  whidi  was  developed  by  die  Anny  Waterways  Experiment  Station 
C  (WES)  as  described  in  the  comments  below  taken  from  the  CX3NWEP  source 
C  code.  The  routines  and  fimctioos  marked  are  also  used  in  this  prograia 
C  The  body  of  this  routine  comes  from  the  CX]NWEP  subroutine  BLAST  which  is 
C  described  below,  after  CONWEP,  in  the  original  comments.  KLAST  was 
C  modified  to  meet  die  specific  needs  of  this  program.  It  provides  die 
C  am^  of  data  required  to  specify  die  pressure  time  history  at  a  point 
C  using  the  Modifieid  Friedlander  Equation.  It  provides  these  data  amys 
C  to  tte  function  F(I,T)  which  returns  the  pressure  at  any  point  on  the 
C  r^kn  of  interest  and  any  time. 

C  The  subroutines  and  functions  asteridced  are  used  in  diis  program 
C  as  well. 

C 

ccccxxccxxcxxxaxxxxxxctxccxxxxxxxxxxxxxxxxxxxcxxxxxcxxxxxxxxccc 

cccccxxxxxxc 

c 

C  PROGRAM  ConWep 

C 

C  SOURCE  USAEWES  /  SS  (Structures  Division) 

C  D.  W.  Hyde 

C  (601)634-2758 

C 

C  PURPOSE  Calculate  a  variety  of  conventional  weapons  effects 
C  according  to  TMS-8SS-1,  "Fundamentals  of  Protective 

C  Design  for  Conventional  Weapons” 

C 

C  LAST  UPDATED.  April  13,  1989 
C 

C  REQUIRED  SUBROUTINES: 

C  ABOTTS  -  Airblast  options  menu  for  PANTO,  PANDW,  &  RANDW 
C  *  BLAST  -  Find  impulse,  time  of  arrival,  duration,  etc. 

C  CASED  -  Brompts  user  for  weapon  parameters  for  cased  weqxms 
C  CRATER  -  Calculate  crater  dimensions. 

C  *  TOCAY  -  Find  decay  coefficient  for  Friedlander  equation. 

C  FRAG  -  Calculate  fragmentation  effects. 

C  INSITO  -  Calculate  quasi-static  pressure  from  an  internal 
C  explosion. 

C  JHOOK  -  Calculate  path  of  a  projectile  through  soil. 

C  LOGO  -  Write  intoductory  screen. 

C  MARRAY  -  Nfriltiply  an  array  of  values  by  a  constant.  Used  for 
C  poforming  unit  conversions  prior  to  plotting. 

C  NEWTNT  -  Prompts  user  for  explosive  parameters. 

C  NOSE  -  Dsplays  nose  shape  ^ors  for  common  projectile  shapes. 

C  PANTO  -  Given  peak  pressure  &  range,  find  chai]^  wei^ 

C  PANDW  -  Given  peak  pressure  &  charge  wei^  find  range. 

C  PAUSE  -  Pronqits  user  to  press  <CR>  before  continuing 
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PEN  -  Calculate  pene(nti<n  of  a  projectile 
various  taq^  matrrials 

RANDW  -  Given  and  charge  find  peak  pressure. 

*  RCCXF  >  Blodc  data  routine  cordaining  digitized  values  for 

reflected  pressure  coefficients  at  difibent  angles 
(tf  inckfcnoe. 

REFLBC  -  Finds  pressure  and  impulse  distributions  on  a  vwll. 
SHAre  -  Find  penetration  of  shaped  charges  into  various  matls. 
SHOCK  -  Calculate  pressure,  dispbcement,  acceleration  due 
to  a  buried  charge  detonation. 

SMARMS  -  hompts  user  for  small  arms  prc^ectile  parameters. 
SIDOUT  -  Standard  output  routine  for  all  routines  that  produce 
X,Y  amys 

SIUNTT  -  IVompt  user  for  unit  system  to  ua^  set  conversion 
foctors. 

SUBSUR  -  Cakailate  airUast  due  to  a  shaUow4wried  cstplosion. 
ITiT  -  Di^lays  a  list  of  common  ex{dosives  and  associated 
constants. 

TUNNEL  -  Calculate  aiiblast  attenuation  in  a  tunnel. 

WEAPCM  -  Select  a  weapon  for  use  in  other  sidxoutmes. 

WEAPS  •  Block  Data  routine  containing  parameters  for  weapons 
catalog. 

REQUIRED  FUNCTIONS  : 

*  nNC  •  Given  a  scaled  lar^  and  burst  config^oation,  retuns 

tile  pedr  inddeit  pressure  in  psi. 

PREF  •  Returns  the  normally  reflect  pressure  in  psi. 

PVEL  •  Returns  the  tiiodc  fiont  velocity  in  kfps. 

*  TARR  •  Returns  the  scaled  arrival  time  in  msec  /  lb**l/3 

*  TDUR  •  Returns  the  scaled  positive  phase  duration  in  ms/lb**l/3 

XIMPR  -  Returns  the  scaled  reflected  impulse  in  psi*ins/ib**l/3 

*  XIMPS  -  Returns  the  scaled  iiKident  impulse  in  psi*ms/lb**l/3 
ZFROMP  -  Given  an  incident  pressure  (psi)  and  burst  configuration, 

returns  the  scaled  range  in  ft/Ib**I/3. 


XCCCCCCCCCCCCCCCCCCCCCCXXXXXXXXXXCCOCOCCCOOCO 

:cccccccccccccccccc 

SUBROUTINE  BL\ST1(SURF) 


SOURCE:  USAEWES  /  SS 
D.  W.  Hyde 
(601)634-2758 


LAST  UPDATED.  27  April  1988 


PURPO!^  Use  the  equations  fiom  BRL  Technical  Report 
ARBRL-TR-Q2S5S  to  find  the  incident  and  reflected 
impulse,  reflected  pressure,  time  of  arrival,  duration, 
and  pressure-time  history  for  a  conventional 
exitiosion. 
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C  DESCRIFIIQN  OF  VARIABLES: 

C  A  -Decay  coefficient  of  incident  preasure 

C  vs.  time  hisloiy,  vdiere 

C  P(t)*PSO*(l-T/TO)*exp(-A*T/IO) 

C  B  -  Dec^  coefficient  of  fefkic^  |xessu^ 

C  vs.  time  histofy 

C  IMPO  -  Peak  incident  impulse,  psi-msec 
C  IMPR  -  Peak  reflected  impulse,  psi-msec 
C  PSO  -  Peak  incident  overpressure,  psi 
C  n(I,i)  -  Incident  pressure  at  tone  station  i 
C  n(24)  -  Incident  impulse  at  time  station  i 
C  PI^i,i)  -  Reacted  pressure  at  time  station  i 
C  m(2,i)  -  Reflected  impulse  at  time  station  i 
C  mo  -  Peak  reflected  pressure,  psi 
C  R  -Raoge,ft 

C  SURF  -  .TRUE  for  surfooe  burst,  JALSE  for  air  burst 
C  T(l,i)  -  Time  at  station  i 

C  Tl(2,i)  -  Identical  to  T(l,iX  kqjt  for  compatilMlity  W  plot 
C  routineonly 

C  TA  -  Arrival  time,  msec 
C  TO  -  Positive  phase  duration,  msec 
C  WTNT  -  Equivalent  weig^  (rf^TNT,  lb 
C  Z  -  Scaled  r^  ft/lb**l/3 
C  ZLXXj  -  Logarithm  of  scaled  rar^ 

C 

CXXXXXXCOOCXXIOOOOCXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXDCCOOOCCCC 

ccxxxxxxxxxxr 

CHARACTER  NUM*3 

REAL  MAXP14O,CONST,SSPOLY(10),STAT(10)4jOGRS 
&  ,YINC1,YINC2 
REAL*8  DECAY 
INTRINSIC  EXP 
C 

common  /in/  xval,zval,wval,iter,ts,ij4ium,iM(,d,fc4y 

COMMON  /PARAMS/  PRV(70),AV(70),TAV(70XTDV(70),XL,YDATA(70) 

COMMON  /REFCCVRa20,39),ATABLE(39)J>TABLE(20) 
open(unit=20,file='Ioad'//num) 

C 

C  BLOCK  DATA  REFLEC 

C  THIS  BLOCK  DATA  SUBPROGRAM  CONTAINS  REFLECTED  PRESSURE  OOEFFIdENTS 
C  FROM  RGURE  3-3  OF  TM5-855.  THESE  DATA  WERE  DIGITIZED  BY/FOR 
C  PERSONNEL  OF  THE  NAVAL  CIVIL  ENGINEERING  LABORATORY  FOR  USE  IN  THE 
C  MICRO  PROGRAM  SHOCK  THE  ORIGINAL  BLOCK  DATA  SUBPROGRAM  CONTAINED 
C  THE  SINES  OF  THE  INVERSE  OF  THE  COEFHCIENrS;  THIS  VERSION  HAS  THE 
C  ACTUAL  COEFndENTS. 

C  COMMON/REFCCyRC(2039),ATABLE(39XPTABLE(20) 

C  5000  PSI 

DATA  (RC(  U)4=U9)  /1225,12.00,1 1.65,1120,10.65,10.05, 9.40, 

&  8.85,  8.55,  825, 7.95,  7.90,  7.85,  7.80,  7.60,  7.50,  735,  7 35, 

&  7.50, 9.00,  9.05,  9.05,  8.35,  720,  4.70,  2.45,  220,  2.05,  1.45, 
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A  US,  UO.  120, 1.IS,  1.12, 1.Q2,  1.00. 1.00;  1.00;  1.0(V 
C  3000  PS 

DATA  (RCX  24)4-139)  /10.80,10.70,1035.10j0. 9.95, 9.40,  8.80, 
A  835, 8.00, 7.75, 730, 7.45, 7.40, 735, 730, 735, 7.40,  830, 

&  8.60,  8.65,  8.60,  8.50,  7.70, 6.40. 435, 145, 230, 105,  1.45, 

A  135, 130, 130,  1.15, 1.11  101 100, 1.00;  1.00, 1.00/ 
C2000PSI 

DATA  (R0(  34)4-139)  /10.00, 9.95, 9.80, 935. 930, 8.75,  830, 
&  7.75,  7.50,  735,  7.05,  7.00, 6.95, 6.90, 635, 6.85,  7.10,  8.00, 

A  830,  835, 8.15, 8.05,  735, 635, 430, 145, 230, 105, 1.45, 

A  135, 130, 130, 1.15, 1.11  101  100. 1.00, 1.00, 1.00/ 

C  1000  PS 

DATA  (RQ  44)4-139)  /  8.60. 8.45,  835, 735, 7.65,  730, 630. 
A  6.60;  6.40, 635, 6.15, 6.10, 6.05, 6.05, 6.00;  6.15, 7.10;  7.60; 

A  735, 7.45, 735, 730, 635, 5.75, 3.90, 145, 230, 105, 1.45, 

A  135, 130, 130, 1.15, 1.ll  1.01  100,  1.00, 1.00,  1.00/ 

C  500  PS 

DATA  (RC(  54)4-139)  /  7.80,  7.65,  7.45, 735, 6.95, 6.70, 635, 
&  6.05,  5.85,  5.75,  5.70,  5.75,  5.80,  5.80, 6.05, 6.50,  7.10,  7.00, 

A  6.85, 6.75, 6.65, 635, 6.00,  535, 3.75, 145, 120, 105,  1.45, 

&  135, 130, 130, 1.15, 1.11  101  100, 1.00, 1.00, 1.00/ 

C  400  PSI 

DATA  (RCX  64)4-139)  /  7.00, 6.95, 6.85, 6.75, 6.55,  6J0, 6.00, 
&  5.753.633333.63.73.8,6.05,635,6.6,63, 
&6.4,63,63,6.13.6.4.93.63.45333.0S,1.45, 

A  US,  130. 130, 1.IS,  1.11  \Sa,  1.00, 1.00;  1.00;  1.00/ 

C  300  PS 

DATA  (RCX  74)4=139)  /  6.65, 6.60, 630, 635, 630, 5.90, 5.65, 
&  5.40,  530,  535,  5.40,  5.50,  5.60,  5.80, 6.00, 6.05,  5.95,  5.85, 

&  5.75,  5.70,  5.60,  535,  5.15, 4.60,  3.50, 145, 230,  lOS,  1.45, 

A  US,  UO,  130,  1.15,  1.11  101  100,  1.00,  1.00, 1.00/ 

C  200  PSI 

DATA  (RQ  84X1=139)  /  6.00,  5.90,  5.85, 5.75,  5.65,  5.45,  530, 
&  4.95, 4.90, 4.95, 5.15,  530,  5.40,  530,  5.45, 5.40,  535,  5.10, 

&  5.05,  5.00, 4.95,  4.85,  4.55, 4.10,  3J0, 145, 235, 230, 1.90, 

&  1.85,  1.80,  1.60,  1.55,  1.50,  1.40,  135,  130,  1.15,  1.00/ 

C  150  PSI 

DATA  (RC(  94)4=139)  /  5.60,  530,  5.45, 535, 530,  5.05, 4.80, 
&  4.60, 4.60, 4.65, 4.90,  5.10,  5.15,  5.10,  5.05, 4.90, 4.75, 4.65, 

A  435, 4.45, 4.40, 435, 4.05, 3.70, 3.00, 230, 115, 2.10,  1.80, 

&  1.75,  1.70,  1.65,  1.60,  130,  135,  130,  1.15,  1.10,  1.00/ 

C  100  PS 

DATA  (RC(104)4=139)  /  5.00,  4.95, 4.85, 4.80, 4.65,  4.55,  4.40, 
&  430, 435, 430, 4.65,  4.75, 4.70,  4.60, 4.40, 430,  3.95,  3.75, 

A  3.70, 3.60,  3.55, 3.45, 335,  2.95, 2.40, 2.00,  1.90,  1.85,  1.70, 

&  1.65, 1.60, 130, 1.45, 1.40, 130, 135, 132, 1.10, 1.00/ 

C  TOPS 

DATA  (RQl  U)4=U9)  /  4.45, 4.40, 435, 435,  4.15,  4.10, 4.00, 
&  3.85, 3.85,  3.95, 430, 435, 435, 4.10, 3.85, 3.60, 335, 330, 

&  3.10, 3.05, 3.00, 2.90, 2.70, 345, 305,  1.80,  1.75,  1.70,  1.55, 

&  1.50,  1.45,  1.35,  130,  135,  133  130,  130,  1.10,  1.00/ 


C  SOPS 

DATA  (RC(12jy-U9)  /  4.00, 3.90. 3.85, 3.80, 3.70, 3.65, 3.55, 
A  3.50, 3.50, 3.60, 4.10, 4.10, 4.00, 3.80, 3i0, 3.25, 195, 180, 

&  170, 2.65, 160, 155, 130, 2.05,  1.85,  IJSS,  1  JO,  1.45, 135, 

&  131 130, 115,  m  110,  1.15,  1.11  1.10,  1.05,  1.00/ 

C  30PS 

DATA  (RC(134)4=139)  /  335, 330, 315, 310, 3.15, 3.10, 3.10, 
&  3.10, 3.15, 310, 3.50, 3.60, 3.55,  3.40, 310, 2.95, 175, 2.60, 

A  2.55, 150, 2.45, 2.40, 210, 2.00, 1.70,  135, 130, 1.45, 135, 

A  131  130, 115,  m  110, 1.15, 1.11  1.10,  1.05,  1.00/ 

C  20PS 

DATA  (RC(144)4=139)  /  3.00,  3.00,  3.00, 195, 190, 185, 185, 
&  190, 2.90, 195, 3.05, 3.10, 310, 315, 3.40, 315, 3.00, 190, 

A 185, 175, 170, 165, 145, 210, 1.90,  1.70, 135, 1.60;  130, 

A  1.45, 1.40, 130, 117, 115, 122, 110, 110;  1.10;  1.00/ 

C  lOPS 

DATA  (RC(154y=139)  / 150, 2.50, 230, 150, 230, 230, 230, 
&  230, 150, 2.50, 2.50, 2.55, 2.60, 2.65, 170, 175, 3.00, 310, 

A  315,  315, 310, 3.15, 2.90, 2.60, 210,  1.90,  1.85,  1.80, 1.65, 

A  1.60,  1.55,  1.50,  1.45,  1.40,  130, 115,  110;  1.15, 1.00/ 

C  5PS1 

DATA  (RC(16jy=139)  /  210, 210, 210, 210, 210, 210, 210, 
&  210,  210, 210, 210, 210, 210, 215, 230, 235, 235, 237, 

&  2.40, 2.45, 2.47, 2.52, 2.75, 315, 185, 145, 235, 230, 2.00, 

A  1.95,  1.90,  1.70, 1.65,  1.60,  1.40,  135,  130, 110,  1.00/ 

C  2PS 

DATA  (RCX17jy»139)  / 110, 110, 110, 110, 110, 110, 110, 

&  110, 110, 110, 110, 110, 110, 110, 110, 110, 110, 110, 

A 110, 110, 110, 110, 115, 210, 2.40, 3.10, 3.10, 3.05, 2.60, 

&  150, 2.40, 110,  2.00,  1.90,  1.60,  1.51  147,  110,  1.00/ 

C  1  PSl 

DATA  (RC(18,1),I=139)  /  2.05,  2.05,  2.05, 105,  2.05,  2.05,  2.05, 
&  2.05,  2.05, 2.05, 2.05,  2.05, 2.05, 2.05, 105, 2.05, 2.05, 2.05, 

A  2.05, 2.05, 2.05, 2.05,  2.05, 2.05, 110, 215, 230, 235, 3.10, 

&  3.15, 3.10, 2.65, 2.50, 135,  1.90,  1.75,  1.65,  130,  1.00/ 

C  0.5  PSl 

DATA  (RC(19,1),1=139)  /  2.05, 105, 105, 2.05, 2.05, 105, 105, 
&  2.05,  2.05, 2.05, 2.05, 2.05, 2.05, 105, 105, 2.05, 2.05, 2.05, 

A  2.05, 2.05, 2.05, 105,  2.05, 2.05, 2.05, 105, 2.05, 110, 215, 

&  230, 2.40, 3.05, 3.10,  3.05, 230, 2.15, 100,  1.45,  1.00/ 

C  01  PSl 

DATA  (RC(204XI=139)  /  2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 100, 

&  2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 100, 2.00, 2.00, 

A  2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 2.00, 100, 

&  2.00, 2.05, 2.10, 115, 210, 3.05, 3.05, 195,  1.70,  1.00/ 

DATA  ATABLE/  0.0,  5.0,10.0.15.010.015.030.034.036.038.0; 
&  40.0,40.5,41 .0,41 .5,42.5,43.5,45.0,46.0,46.5,47.0,47.5,48.0, 


&  79.0,80.0,81.0,85.0,90.0/ 

DATA  PTABLE/  5000.3000.1000., 1000.,  500.,  400.,  300.,  200., 
A 150.,  100.,  70.,  50.,  30.,  20,  10  .  5.,  2,  1., 
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ft  0.5,  0.2/ 

C 

C 

YINCl -XL/98 
YINC2 -XL/38 
YDATA(l)-0 
DO  5  I  -  2,69 

IF  a  IE  50)  THEN 
J-1-1 

YDATA(I)  -  YDATA(J)  +  YINCl 
ELSE 
J  =  I-1 

YDATA(I)  =  YDATA(J)  +  YINQ 
QCIF 

5  CONIINUE 
C 

C  CXDMPLTIE  CUBE  ROOT  OF  TNT 
WJ  =  TNT**(173.) 

C 

DO  10 1-1,70 

RANG  =  SQRTPC*X  +  Z*Z  +  YDATA(I)*YDATA(I)) 

RANGS-RANG/W3 
C  COMPUTE  ANGLE  OF  INODB^ 

CO  =  AMINl(l.,Z/RANG) 

RT0D-45yATAN(l.) 

ALPHA  -  AO0S(OO)^RT0D 

C  FIND  INCIDENT  PRESSURE  AT  THIS  POINr,  NEED  THE  LOG  OF  THE  SCALED  RANGE 
LOGRS  =  ALOG10(RANGS) 

MAXn  =  PINC(LOGRS) 

C 

C  FIND  PRESSURES  FROM  TABLE  THAT  BOUND  PSO 
rm  "Jon  kk  =  i 

1F(MAXP1  .GE  PTABLE(KK))  GOTO  201 

200  CONITNUE 

201  IPI  -  MAX0(IJCK-1) 

IP2  =  MIN0(20.KK) 

C  FIND  ANGLES  FROM  TABLE  THAT  BOUND  ALPHA 
DO205KK-  1,39 

1F(ALPHA  .LE  ATABL£(KK))  GOTO  206 

205  CONITNUE 

206  lAl  =  MAX0(]  JCK-I) 

1A2  =  M1N0(39,KK) 

1F(1A1  .EQ.1A2)THEN 
FACTA  =  0.0 
ELSE 

FACTA  =  (ALPHA  -  ATABLE(IA1))  /  (ATABLB(IA2>.ATABLE1(IA1)) 

ENDIF 

Cl  =  RCXIPIJAI)  +  (  RCaP14A2)  -  RC(1P14AI)  )  *  FACTA 
C2  -  RC(IP2M1) +  ( RC(IP2M2)  -  RC(IP2M1)  )•  facta 
IFOPl  .BQ.  1P2)  THEN 
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REF-Cl 

pi.q? 

Rff  »  Cl  +  (C2-C1)  •  (  MAXPI  -  PTABLEOPl) )  / 
ft  ( FTABLE(IF2>FTABLE(IPI) ) 

ENDIF 

C 

PRV(D  =  MAXPI*REF 
CCCCCOCCC 

COCCOOCCC  INSERT  FUNCTIONS  FROM  CONWEP 
OCCCCCOOC 

IO»XIMPS(1jOGRS)*W3 
TAV(I)  =  TARR(LOGRS)*W} 

TOVG)  =  TIXJR(LOGRS)*W3 
ID-'TOVG) 

AV(I)  -  DECAYCMAXPUQTD) 
write^O^)  i4)rv(tXtav(i),t0v(iXt(Uv(t) 

10  CONIINUE 
22  FORMATKIS^FllS) 

DO  15  J-1.69 

YDATA(J)  -  YDATA(jyXL 
15  CONTINUE 
C 

RETURN 

C 

END 

COOOOCOOOOCXXXXXXOCOOOCOCCCCCCOOOOCCOOOCXXXXXXXXXOOOOOOOOOOOOOOC 
REAL  FUNCTION  IWaZLOG) 

C 

C 

C  PURPOSE:  F^^ID  THE  INODENT  PRESSURE  DUE  TO  TIIE  DETONATION 
C  OF  A  1  LB  EQUIVALENT  TNT  CHARGE  EQUATIONS  ARE  FROM  BRL 
C  TECHNICAL  REPORT  ARBRL-TR-02555.  PRESSURE  IS  RETURNED  IN 
C  UNITS  OF  PSI. 

C 

C  DESCRIPTION  OF  VARIABLES: 

C  SURF  -  .TRUE  FOR  SURFACE  BURST,  .FALSE  FOR  AIR  BURST 
C  ZLOG  -LOGARITHM  OF  SCALED  RANGE 
C 

PARAMETER  (NS=1 1,  NF=8) 

REALCSURF(NS+1) 

DATA  CSURF  /  1.9422502013,  -1.6958988741, 
ft  -0.154159376846,  0.514060730593, 

ft  0.0988534365274,  -0293912623038, 

ft  -0.0268112345019,  0.109097496421, 

ft  0.00162846756311,-0.0214631030242, 

ft  0.0001456723382,  0.00167847752266/ 

C 

U  =-0.756579301809+  1.35034249993*ZLOG 
PINC»CSURF(NS+1) 

DO10I=NS,l,-l 
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10  P!NC-P!NC*U  +  CSURF(I) 

P!NC-10**PINC 

RETURN 

END 

ccxxxxxxxxxxxxxmxxxccxxxxxxxxxcxxxxxxxxxxcocccxxxxxxxcxxxxxxxx: 

DOUBLE  PREaSION  FUNCIION  DE1CAY(F04Q.1D) 

REALM  P0J0.1D 
REALM  A^'AJTAJTOI 
C 

C  FIND  RATE  OF  DECAY  FOR  PRESSURE  ASSUMING; 

C 

C  PCD  =  PO  •  [  1  -  (T-TAyiD]  *  EXP[-A*(T-TAyiDl 

C  (FRIEDLANDER!S  EQUATION) 

C 

C  WHEREAISADBCAYOOEFFIOQrr.  INIBGRAT1NG  THIS  EQUATION 
C  OVER  TIME  GIVES  THE  IMPULSE; 

C 

C  I0  =  P0*TD*[A  +  EXP(-A)- 1  ]/A**2 
C 

C  F1NDF(A)  =  A**2-P0M1>I0*[A  +  EXP(-A)- 1]»0 
C  OR: 

C  A**2/[A  +  EXP(-A).  l]-PO*TTyi0»0 

C 

C  FOR  LARGE  A  EXP(.A)  APPROACHES  0.  AND 
C  A**2/(A- 1) APPROACHES A+ I 
C 

C  INITIAL  GUESS  A  *  P0*TIVI0. 1 
C 

PTOI  =  P0*TI>10 
C  INITIAL  GUESS: 

A  =  PT0I-  I. 

1  FA  =  A*A  -  Froi*(A  +  EXP(-A)  -  I.) 

FPA  =  2*A  -  PTOI*(l.  -  EXP(-A)) 

A  =  A  -  FA/FPA 

IF(ABS(FA)  .GT.  l.I>6)  GO  TO  I 

DECAY  =  A 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCOCOOCCCCCCOOCCCCCCCCCCOCCCOOOCOCCCCOOCCC 
REAL  FUNCTION  TARR(ZLOG) 

C 

C  PURPOSE:  FIND  THE  SCALED  TIME  OF  ARRIVAL  FOR  THE  DETONATION  OF  A 
C  1  LB  EQUIVALENT  TNT  CHARGE  EQUADONS  ARE  FROM  BRL 
C  TECHNICAL  REPORT  ARBRL-TR-02555.  ARRIVAL  TIME  IS  RETURNED 

C  IN  MSBC/LB**(l/3). 

C 

C  DESCRIPTION  OF  VARIABLES: 

C  ZLOG  -LOGARITHM  OF  SCALED  RANGE 
C 

PARAMETER  (NS^) 

REALCSURF(NS4-1) 
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CATACSURFA0.173607601251.  1J57064962SS. 

A  0.0S249279«64S.  -0.196563954086^ 

ft  -0.0601770052288,  0.0696360270891, 

ft  00215297490092,  -0.0161658930785, 

ft  -0.00232531970294, 0.00147752067524  / 

U  =  -0.755684472698  +  1 J7784223635*ZIXX} 

TARR  =  CSURF(NSfl) 

DO10I  =  NS,l,-l 

10  TARR  =  TARR*U  +  CSURF(I) 

TARR*  10.**TARR 
REIURN 

END 

(XXXXXXXCCXXXXXCCXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXCXXXCXXXXXXXX: 

REAL  FUNCnON  1DUR(ZLOG) 

C 

C  PURPOSE  FIND  THE  SCALED  DURAHON  FOR  TOE  DETONAHON  OF  A  1  LB 
C  EQUIVALENT  INT  CHARGE  EQUATIONS  ARE  FROM  BRL  TECHNICAL 

C  REPORT  ARBRL-1R.-02555.  DURATION  IS  REIURNED  IN 

C  MSEC/LB**(l/3). 

C 

C  DESCRIPTION  OF  VARIABLES: 

C  ZLOG  -LOGARITHM  OF  SCALED  RANCX 
C 

PARAMETER  (NS1=5J4S2=8JSIS3=5) 

REAL  CSURF1(NS1+1X  CSURF2(NS2+1X  CSURF3(NS3+1) 

DATA  CSURFl  /  -0.728671776005,  0.130143717675, 
ft  0.134872511954,  0.0391574276906, 

ft  -0.00475933664702,-0.00428144598008  / 

DATACSURF2  /  0.20096507334,  -0.0297944268976, 
ft  0.030632954288,  0.0183405574086, 

ft  -0.0173964666211,-0.00106321963633, 

ft  0.00S62060030977, 0.0001618217499, 

DATACSURR/  0.572462469964,  0.0933035304009, 

ft  -0.0005849420883,  -0.00226884995013, 

ft  -0.00295908591505, 0.00148029868929  / 

1F(ZL0G  -LT.  0.4048337)  THEN 
U  = -0.1790217052  +  525099193925*ZLOG 
TTXJR  =  CSURF1(NS1+1) 

DOll  I=NS1,1,-1 

11  TDUR  =  TDUR*U  +  CSURFia) 

ELSEIF(ZLOG  .GE  0.4048337  .AND.  ZLOG  .LT.  0.845098)  THEN 
U  =  -5.85909812338  +  929%28861  l*ZLOG 
TDUR  =  CSURF2(NS2+1) 

D012I=NS2,1,-1 

12  TDUR  =  TDUR*U  +  CSURF2a) 

ELSE 

U  =  -4.92699491 141  +  3.46349745571 ‘ZLOG 
TDUR-CSURF3(NS3+1) 

D013I=NS3,1,-1 

13  TTXJR  =  TDUR*U  +  CSURF3a) 
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QOF 

TOUR- 10**T1XJR 

RETORN 

END 

(XXX(XCXXXXXCXXXXXXXCXXOCXXXXXXXXXXXXCXXXXXXXXXXXXXXXCXXXXOCCXXX 
REAL  FUNCTION  XIMPS(ZLOG) 

C 

C  PURPOSE:  FIND  TNE  SCALH)  INCTDENT  IMPULSE  FOR  THE  DETONATION  OF 
C  A  1  LB  EQUIVALEt^rr  TNT  CHARGE  EQUATIONS  ARE  FROM  BRL 
C  TECHNICAL  REPORT  ARBRL-TR-Q25SS.  IMPULSE  IS  RETURNED  IN 
C  UNITS  OF  PSI*MSBC/LB**(l/3). 

C 

C  DESCRIPTION  OF  VARIABLES: 

C  ZLOG  -LOGARITHM  OF  SCALED  RANGE 
C 

PARAMETER  (NS1-4NS2-7) 

REAL  CSURF1(NS1+1XCSURF2(NS2+1) 

DATA  CSURFl  /  1.S71S9240621.  -0JQ2992763686, 

&  0.17133S64S23S,  0.0450176963051, 

&  -0.0118964626402/ 

DATA  CSURF2  /  0.719852655584,  -0384519026965, 

&  -0.0280816706301,  0.00595798753822, 

&  0.014544526107,  -0.00663289334734, 

&  -0.00284189327204, 0.0013644816227  / 

IF(21XXi  LT.  0382017)  THEN 
U  -  0.832468843425  -t-  3.Q760329666*ZLOG 
XIMPS  =  CSURF1(NS1+1) 

DOll  I=NS1,1,-1 

1 1  XIMPS  -  XIMPS*U  +  CSURFIG) 

ELSE 

U  =  -2.91358616806  +  2.40697745406*ZLOG 
XIMPS  =  CSURF2(NS2+1) 

DO  12  1  =NS2,1,-1 

12  XIMPS  =  XIMPS*U  +  CSURF2(I) 

ENDIF 

XIMPS  =  10.**XIMPS 

RETURN 

END 

CC(X(XC(X(XC(X(X<X(XCXX(XCXXXXXXXXXXX(XC(X^ 

C  The  remainder  of  the  code  (xxitains  the  first  oder  siiiiultaneous 
C  equatitm  solver  DGEAR  and  its  ^sociated  subroutines. 
CCCCC(XCrCCC(XCCCCCCCC(XCC(XCCC<XCXXXXXXCXXX(X(XCXXX^ 


^G0'40Mni(^UJ(0MOU>0>'s]0\U1it!>U>t0l->OU>00>J(Mnil^U)t0H 
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Anxndix  C 

ReqpODM  Surf  ac«  Model  InfomiatloQ 


500  lb  RMqnoM  Surface  -  MOMLZN  DA3A 


1  ZVAL 

2  XVAL 

3  wjr 

4  W_S 

5  W_H 

6  M_D 

7  SDIR 

8  SDIAG 

15 

0 

12 

3000 

168 

19.20 

7.01 

5.81 

15 

0 

12 

3000 

108 

8.17 

10.14 

6.03 

15 

0 

12 

7000 

168 

15.74 

8.45 

7.00 

55 

0 

36 

7000 

108 

0.00 

0.00 

0.00 

55 

0 

36 

3000 

168 

0.22 

19.30 

6.94 

5 

55 

36 

3000 

168 

0.05 

10.23 

5.38 

55 

55 

36 

3000 

108 

0.00 

0.00 

0.00 

5 

0 

24 

7000 

108 

1.84 

44.38 

26.24 

55 

0 

12 

7000 

168 

3.40 

4.61 

3.96 

5 

55 

12 

7000 

168 

1.07 

2.72 

2.29 

55 

55 

12 

7000 

108 

0.00 

0.00 

0.00 

5 

0 

36 

7000 

168 

1.23 

103 . 17 

79.30 

5 

55 

36 

7000 

108 

0.00 

0.00 

0.00 

55 

55 

36 

7000 

168 

0.00 

0.00 

0.00 

5 

28 

24 

5000 

144 

0.40 

15.53 

6.17 

55 

28 

24 

5000 

144 

0.12 

10.02 

4.98 

30 

0 

24 

5000 

144 

1.73 

33.26 

4.07 

30 

55 

24 

5000 

144 

0.07 

8.07 

3.56 

30 

28 

12 

5000 

144 

3.68 

3.29 

2.58 

30 

28 

36 

5000 

144 

0.23 

29.86 

7.78 

30 

28 

24 

3000 

144 

0.89 

16.79 

6.41 

30 

28 

24 

7000 

144 

0.70' 

16.17 

7.02 

30 

28 

24 

5000 

108 

0.21 

18.29 

7.77 

30 

28 

24 

5000 

144 

0.78 

16.86 

6.34 

30 

28 

24 

5000 

144 

0.78 

16.86 

6.34 

5 

0 

12 

7000 

108 

23.00 

10.00 

8.00 

5 

0 

12 

3000 

168 

8.00 

7.00 

5.00 

5 

0 

36 

7000 

108 

0.79 

120.29 

59.17 

5 

0 

36 

3000 

168 

1.49 

65.75 

41.05 
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1000  lb  Raqpoos*  Surfae*  -  NQNLZN  TSOX 


0 

1  ZVMi 

2  XVAL 

3  W_T 

4  W_S 

5  W_H 

6  MJD 

7  SDIR 

8  SDIAG 

1 

5 

0 

12 

3000 

168 

45.00 

10.00 

8.00 

2 

55 

0 

36 

7000 

108 

0.00 

0.00 

0.00 

3 

55 

0 

36 

3000 

168 

1.29 

37.41 

15.61 

4 

5 

55 

36 

3000 

168 

0.24 

17.90 

5.89 

5 

55 

55 

36 

3000 

108 

0.00 

0.00 

0.00 

6 

5 

0 

24 

7000 

108 

6.07 

105.92 

90.44 

7 

55 

0 

12 

7000 

168 

12.18 

5.39 

4.31 

8 

5 

55 

12 

7000 

168 

2.90 

3.31 

2.61 

9 

55 

55 

12 

7000 

108 

0.47 

5.79 

3.39 

10 

5 

0 

36 

7000 

168 

3.99 

176.59 

164.29 

11 

5 

55 

36 

7000 

108 

0.00 

0.00 

0.00 

12 

55 

55 

36 

7000 

168 

0.00 

13.47 

4.21 

13 

5 

28 

24 

5000 

144 

0.95 

26.95 

5.41 

14 

55 

28 

24 

5000 

144 

1.18 

19.20 

4.05 

15 

30 

0 

24 

5000 

144 

5.44 

19.15 

9.30 

1€ 

30 

55 

24 

5000 

144 

0.47 

13.08 

4.50 

17 

30 

28 

12 

5000 

144 

13.70 

7.32 

5.16 

18 

30 

28 

36 

5000 

144 

1.22 

43.96 

24.79 

19 

30 

28 

24 

3000 

144 

3.60 

21.66 

6.76 

20 

30 

28 

24 

7000 

144 

2.91 

41.70 

8.03 

21 

30 

28 

24 

5000 

108 

1.43 

30.44 

16.00 

22 

30 

28 

24 

5000 

144 

3.18 

23.86 

7.52 

23 

30 

28 

24 

5000 

144 

3.18 

23.86 

7.52 

24 

55 

25 

12 

7000 

108 

2.96 

5.41 

3.28 

25 

15 

0 

12 

3000 

108 

13.04 

12.62 

7.59 

26 

25 

0 

12 

3000 

108 

13.77 

11.17 

6.70 

27 

25 

0 

12 

3000 

168 

32.81* 

7.31 

6.03 

28 

55 

0 

12 

3000 

108 

5.88 

5.58 

3.79 

29 

5 

55 

12 

3000 

108 

1.40 

4.29 

2.44 

30 

55 

55 

12 

3000 

168 

7.21 

3.44 

2.71 

31 

5 

0 

12 

7000 

108 

22.00 

7.00 

4.00 

aii^u>(OHou>oo^<riuittkU>toHO 
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2000  lb  SMponM  Surfae*  •  NaNLZN  DlklX 


1  ZVAL 

2  XVAL 

3  wjr 

4  W_S 

5  W_H 

6  M_D 

7  SDIR 

8  SDIAG 

5 

0 

12 

3000 

168 

35.50 

8.00 

6.00 

55 

0 

12 

3000 

108 

2.87 

7.10 

4.02 

5 

55 

12 

3000 

108 

0.78 

3.62 

2.70 

55 

55 

12 

3000 

168 

4.78 

3.09 

2.49 

5 

0 

36 

3000 

108 

1.73 

122.03 

71.97 

55 

55 

12 

3000 

108 

0.52 

4.58 

3.04 

5 

0 

12 

7000 

108 

12.00 

15.00 

9.00 

55 

0 

12 

7000 

168 

6.24 

6.54 

5.30 

5 

55 

12 

7000 

168 

1.71 

3.75 

3.19 

55 

55 

12 

7000 

108 

0.02 

4.54 

2.82 

5 

10 

36 

7000 

168 

0.35 

66.75 

27.87 

55 

0 

36 

7000 

108 

0.00 

0.00 

0.00 

5 

55 

36 

7000 

108 

0.00 

0.00 

0.00 

55 

55 

36 

7000 

168 

0.00 

0.00 

0.00 

5 

28 

24 

5000 

144 

0.66 

20.05 

6.02 

55 

28 

24 

5000 

144 

0.44 

13.00 

4.70 

30 

0 

24 

5000 

144 

3.25 

40.52 

7.94 

30 

55 

24 

5000 

144 

0.21 

9.70 

4.50 

30 

28 

12 

5000 

144 

6.85 

9.44 

6.82 

30 

28 

36 

5000 

144 

0.53 

40.42 

15.61 

30 

28 

24 

3000 

144 

1.75 

24.74 

4.13 

30 

28 

24 

7000 

144 

1.41 

20.50 

6.52 

30 

28 

24 

5000 

108 

0.59 

25.14 

8.34 

30 

28 

24 

5000 

144 

1.54 

25.31 

5.13 

30 

28 

24 

5000 

144 

1.54 

25.31 

5.13 
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500  lb  RMponM  Surfac*  llrxial  mforiiiatloa 


tM0D55D  12R  X  1C  15 -JAN- 94  8:56 

Page  l 

Acc^ted  model  for  reaponse  M_D 
0  1 


1  Model  Name: 

2  Renxxise  Transfoxmatlcti: 

3  Method: 

4  Wel^tS: 

5  Total  Number  of  Cases: 

6  Nunter  of  Predictors: 

7  NUnter  of  l]k]e}q)anded  Ter 

8  Ntnrber  of  Bxcliided  Cases 

9  Error  Degrees  of  Freedom 

10  Standard  Error  of  Respcn 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


hSiSSa 

Untransformed 
Least  Squares 
None 
29 
5 

14 
0 

15 

6.013722 

0.407118 

2.984413 


ICOEF55D  14R  x  4C  .  '  15- JAN- 94  8:56 

Page  1 


Least  Squcures  Coefficients 

0  Term  1  COeff.  2 

,  Response  M, 

Std.  Error  3 

Model  MR551_COPX 

T-value  4  Signif . 

1  1 

-4.062295 

17.936065 

-0.23 

0.8239 

2  Z 

0.652205 

0.345298 

1.89 

0.0784 

3  X 

-0.426211 

0.088692 

-4.81 

0.0002 

4  W- 

-2.505735 

0.583184 

-4.30 

0.0006 

5  W 

0.007857 

0.002558 

3.07 

0.0078 

6  W  H 

0.181769 

0.106889 

1.70 

0.1097 

7  Z*W 

-0.000047 

0.000021 

-2.22 

0.0421 

8  Z*W  H 

-0.002520 

0.001385 

-1.82 

0.0889 

9  Z*W 

0.005176 

0.003065 

1.69 

0.1120 

10  X*W 

0.013339 

0.003119 

4.28 

0.0007 

11  W  *W  H 

-0.000041 

0.000015 

-2.68 

0.0170 

12  W  H*W 

0.004851 

0.002283 

2.13 

0.0506 

13  Z**2 

-0.003385 

0.002456 

-1.38 

0.1883 

14  W**2 

0.022333 

0.009945 

2.25 

0.0402 

No.  cases 

=  29 

R-sq.  = 

0.8681 

RMS  Error  = 

2.984 

Resid.  df 

=  15 

R-sq-ad],  = 

0.7537 

Oond.  NO.  = 

175.6 
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tAN0V55D  7R  X  5C  15-JAN-94  8:51 

Page  i 


Least  Squares  Sumnary  PHKNk,  Response  M  Model  MR551 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  ToteQ(Oorr.) 

28 

1012.616 

2 

Regression 

20 

891.071 

44.554 

2.93 

0.0609 

3 

Linear 

5 

535.077 

107.015 

7.04 

0 . 0083 

4 

Non-linear 

15 

355.995 

23.733 

1.56 

0.2667 

5 

Residual 

8 

121.544 

15.193 

6 

Lack  of  fit 

7 

121.544 

17.363 

7 

Pure  error 

1 

0.000 

0.000 

R-sq.  s  0.8800 
R-sq-adf.  =  0.5799 

Model  obeys  hieraxohy.  The  sum  of  squares  for  linear  terms 
is  conputed  assuming  nonlinear  terms  are  first  removed. 


tMQD55S  12R  x  1C 
Page  1 
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15-JAN-94  8:56 


Accepted  model  for  response  SDIR 
0  1 


1  Model  Name: 

2  Response  Transformaticn: 

3  Method: 

4  Heights: 

5  Total  Number  of  Cases: 

6  Number  of  Predictors: 

7  Number  of  Uneiqianded  Ter 

8  Number  of  Excluded  Cases 

9  Error  Degrees  of  Freedom 

10  Standard  Error  of  Reapcn 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


MR5S4 

untransformed 

Least  Squares 

None 

29 

5 

13 

0 

16 

29.162779 

0.77632 

7.503846 


iaOEF55S  13R  X  4C  .  »  15-JJ«J-94  8:57 

Page  1 

Least  Squares  Coefficients,  Response  S,  Model  MR554 


0  Term 

1  Coeff. 

2  Std.  Error 

3  T-vsQue 

4  Signif. 

1  1 

-52.895570 

40.769655 

2  Z 

1.942077 

0.466188 

3  X 

-2.154419 

0.523765 

4  W 

6.267415 

0.977596 

5  W 

-0.005148 

0.085973 

6  W  H 

-0.076856 

0.263059 

7  Z*X 

0.025219 

0.003388 

7.44 

0.0001 

8  Z*W 

-0.000208 

0.000052 

-3.99 

0.0011 

9  Z*W 

-0.071017 

0.009085 

-7.82 

0.0001 

10  X*W  H 

0.012470 

0.003068 

4.06 

0.0009 

11  x*vr 

-0.035524 

0.006581 

-5.40 

0.0001 

12  W  *W  H 

0.000069 

0.000037 

1.86 

0.0819 

13 

-0.016012 

0.006058 

-2.64 

0.0177 

No.  Ccises 

=  29  R-S(J. 

»  0.9622 

RMS  Error 

=  7.504 

Resid.  df 

=  16  R-sq-ad] . 

=  0.9338 

Oond.  NO. 

=  156.7 

tANOVSSS  7R  X  5C 
Page  l 
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15-JAN-94  8:55 


Least  Squares  Sumnaxy  P^KUA,  Response  S  Model  MR554 


0 

Source 

1  df 

2  Sum  Sg. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  Total  (C3orr.) 

28 

23813.09 

2 

Regression 

12 

22912.17 

1909.35 

33.91 

0.0000 

3 

Linear 

5 

14107.21 

2821.44 

50.11 

0.0000 

4 

Non-linear 

7 

8804.96 

1257.85 

22.34 

0.0000 

5 

Residu2d 

16 

900.92 

56.31 

6 

Lack  of  fit 

15 

900.92 

60.06 

7 

Pure  error 

1 

0.00 

0.00 

R-sg.  «  0.9622 
R-sg-adf.  >  0.933B 

Model  obeys  hierazx±y.  The  sum  of  squares  for  linear  terms 
is  confuted  assuming  nonlinear  terms  are  first  removed. 
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tMQDSSSD  i:!R  X  1C  15-JMI-94  8:56 

Page  l 

Accepted  model  for  response  SDIAG 
0  1 


1  Model  Name: 

2  Response  Transfomiatian: 

3  Metkxl: 

4  Vfeic^ts: 

5  Total  NUnioer  of  Cases: 

6  Nuinber  of  Predictors: 

7  NUnober  of  Uhexpanded  Ter 

8  NUnber  of  Exclxided  Cases 

9  Error  Degrees  of  Freedom 

10  Standard  Error  of  Reqxxi 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


^Si555 

Untransformed 

Least  Squares 

None 

29 

5 

15 

0 

14 

18.270403 

0.791464 

3.622131 
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Least  Squares  Coefficients,  Response  SD,  Model  MR555 


0  Term 

1  Coeff. 

2  Std.  Error 

3  T-value 

4  Signif. 

1  1 

214.375478 

41.805261 

2  Z 

0 . 967757 

0.217279 

3  X 

0.780201 

0.194868 

4  W 

-0.386121 

0.628712 

5  W 

-0.009305 

0.002934 

6  W  H 

-3.129392 

0.549253 

7  Z^X 

0.011766 

0.091568 

7.50 

0.0001 

8  Z*W 

-0.000126 

0.000026 

-4.94 

0.0002 

9  Z*W 

-0.033258 

0.003835 

-8.67 

0.0001 

10  X*W 

-0.000109 

0.000023 

-4.78 

0.0003 

11  X*W“ 

-0.030166 

0.003404 

-8.86 

0.0001 

12  W  *W  H 

0.000056 

0.000018 

3.06 

0.0085 

13  W  *W 

0.000298 

0.000054 

5.56 

0.0001 

14  W  H*W 

0.006501 

0.002973 

2.19 

0.0462 

15  wrH**2 

0.010443 

0.001918 

5.45 

0.0001 

No.  caises 

=  29  R-SC[. 

=  14  R-sq-adf. 

=  0.9803 

RMS  Error  »  3.622 

Resid.  df 

=  0.9607 

Oond.  No. 

s  494.4 
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Least  Squares  Suninary  MKNA,  Response  SD  Mcxiel  MR555 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  Total  (Oorr.) 

28 

9346.614 

2 

Regression 

14 

9162.936 

654.495 

49.89 

0.0000 

3 

Linear 

5 

5307.318 

1061.464 

80.91 

0.0000 

4 

NOn-linesu: 

9 

3855.618 

428.402 

32.65 

0.0000 

5 

Residual 

14 

183.678 

13.120 

6 

Lack  of  fit 

13 

183.678 

14.129 

7 

Pure  error 

1 

0.000 

0.000 

R-sg.  «  0.9803 

R-sg*adf.  ■  0.9607 

al  obeys  hierardiy.  Ihe  sum  of  squares  for  linear  terms 
coRputed  assuming  nonlinear  terms  are  first  removed. 
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Accqpted  model  for  rMponae  M_D 
0  1 


1  Model  Name: 

2  Repnonse  Transfomatlon: 

3  Meunod: 

4  Hei^lts: 

5  Total  Nbnober  of  Cases: 

6  Ndnber  of  Predictors: 

7  NUnloer  of  Itoeaipanded  Ter 

8  NUiriaer  of  Excluded  Cases 

9  Error  Degrees  of  Freedom 

10  Standard  Error  of  Respcxi 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


PLQfI2002 

Dkxtransformed 

Least  Squares 

None 

25 

5 

14 

0 

11 

7.242622 

0.208364 

1.690644 
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Least  Squares  Coefficients,  Response  M,  Model  PlLC7r2002 


0  Term 

1  Coeff.  2 

Std.  Error 

3  T- value  4  Signif. 

1  1 

-5.969702 

7.072281 

2  Z 

-0.187083 

0.108461 

3  X 

-0.186465 

0.102248 

4  W 

-C  525068 

0.327016 

5  W 

-0.001926 

0.000538 

6  W  H 

0.437493 

0.045519 

7  Z^X 

0.006853 

0.000781 

8.78  0.0001 

8  Z*W 

0.011265 

0.001819 

6.19  0.0001 

9  Z*W  H 

-0.002346 

0.000720 

-3.26  0.0076 

10  X*W 

0.01190C 

0.001793 

6.64  0.0001 

11  X*W  H 

-0.002847 

0.000658 

-4.33  0.0012 

12  W*W 

0.000056 

0.000028 

2.02  0.0686 

13  W*W“H 

-0.011774 

0.001523 

-7.73  0.0001 

14  W**5 

0.019214 

0.005685 

3.38  0.0061 

No.  cases 

®  25  R-sq. 

=  0.9750 

RMS  Error  »  1.691 

Resid.  df 

=  11  R-sq-ad]. 

=  0.9455 

Cond.  No.  =  115.9 

tANCfUAZ  7R  X  5C 
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Least  Squares  Sunmaiy  PXKNA,  Respcnse  M  Model  PL0T2002 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  P-Ratio 

5  Signif. 

1  Total  (Oorr.) 

24 

1258.934 

2 

Regression 

13 

1227.493 

94.423 

33.03 

0.0000 

3 

Linear 

5 

656.755 

131.351 

45.95 

0.0000 

4 

Non-linear 

8 

570.738 

71.342 

24.96 

0.0000 

5 

Residual 

11 

31.441 

2.858 

6 

Lade  of  fit 

10 

31.441 

3.144 

7 

Pure  error 

1 

0.000 

0.000 

R-sq.  -  0.9750 
R-sq-adf .  «  0.9455 

Model  obeys  hierarchy.  The  sum  of  squares  for  linear  terms 
is  conputed  assuming  nonlinear  terms  are  first  removed. 
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Accepted  model  for  response  SDIR 
0  1 


1  Model  Name: 

2  Response  'n^nsformatlGn: 

3  Met&xl: 

4  Weights: 

5  Total  Number  of  Czises: 

6  NUnber  of  Predictors: 

7  Number  of  Uheaqpanded  Ter 

8  Number  of  Excluded  Cases 

9  Error  Degrees  of  Freedom 

10  Standard  Error  of  Reapon 

11  Relative  IttBSS: 

12  Root  Mean  Squared  Error: 


PLCrr2003 
Utitransformed 
Least  Squares 
None 
25 
5 
14 
0 

11 

26.471305 

0.880064 

3.393651 
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Least  Squares  Coefficients,  Response  S,  Model  FL0T2003 


0  Term 

1  COeff. 

2  Std.  Error 

3  T- value  4  Signif. 

1  1 

118.627133 

20.455120 

2  Z 

-0.640760 

0.255962 

3  X 

-0.560721 

0.120671 

4  W 

4.347709 

0.345749 

5  W 

-0.020032 

0.002916 

6  VTH 

-1.021689 

0.113097 

7  Z*X 

0.014785 

0.001509 

9.80 

0.0001 

8  Z*W 

-0.029421 

0.003347 

-8.79 

0.0001 

9  Z*W  H 

0 . 008676 

0.001782 

4.87 

0.0005 

10  X*W 

-0.043227 

0.004627 

-9.34 

0.0001 

11  X*W 

0.000106 

0.000025 

4.31 

0.0012 

12  w*vr 

-0.000246 

0.000064 

-3.82 

0.0028 

13  W  *R  H 

0.000151 

0.000017 

8.97 

0.0001 

14  Z^*2 

-0.010774 

0.002796 

-3.85 

0.0027 

No.  caises 

=  25  R-SCJ. 

=  0.9925 

RMS  Error 

=  3.394 

Resid.  df 

=  11  R-sq-ad^. 

=  0.9836 

Oond.  No. 

=  160.5 

tANCfUA3  7R  X  5C 
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Least  Squares  Sumnary  KtKNh,  Response  S  Ftodel  PL0T2003 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  Total  (Corr.) 

24 

16817.52 

2 

Recpression 

13 

16690.83 

1283.91 

111.50 

0.0000 

3 

Linear 

5 

9335.91 

1867.18 

162.10 

0.0000 

4 

Non-linear 

8 

7354.92 

919.37 

79.83 

0.0000 

5 

Residual 

11 

126.69 

11.52 

6 

Lack  of  fit 

10 

126.69 

12.67 

7 

Pure  error 

1 

0.00 

0.00 

R-sg.  «  0.9925 

R-sq-adJ.  «  0.9836 

Model  obeys  hierarchy.  The  sum  of  squares  for  linear  terms 
is  conputed  assuming  nonlinear  terms  are  first  removed. 
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Accepted  nodel  for  xespcnse  SDIAG 
0  1 


1  Model  Name:  PLQT2004 

2  Response  Transformation:  untransfoimed 

3  Method:  Least  Squares 

4  Weights :  Ncne 

5  Total  Number  of  Cases:  25 

6  Number  of  Predictors:  5 

7  Number  of  unexpended  Ter  14 

8  Number  of  Excluded  Cases  0 

9  Error  Degrees  of  Freedom  11 

10  Standard  Error  of  Respcn  14.366108 

11  Relative  PRESS:  0.481828 

12  Root  Mean  Squared  Error:  4.099204 
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Least  Squares  Coefficients,  Response  SD,  Model  PL0T2004 


0  Term  1  Coeff. 

2  Std.  Error 

3  T-value  4 

Signif . 

1  1 

81.481744 

22.482048 

2  Z 

-0.311996 

0.156940 

3  X 

-0.896116 

0.300395 

4  W 

0.981792 

0.626264 

5  W 

-0.009560 

0.003540 

6  W  H 

-0.557030 

0.130852 

7  Z*X 

0 . 007351 

0.001821 

4.04 

0.0020 

8  Z*W 

-0.025476 

0.006004 

-4.24 

0.0014 

9  Z*W 

0.000083 

0.000034 

2.46 

0.0319 

10  X*W 

-0.009849 

0.003876 

-2.54 

0.0274 

11  X*W  H 

0.005363 

0.002045 

2.62 

0.0237 

12  W*W 

-0.000220 

0.000076 

-2.92 

0.0140 

13  W  *W  H 

0.000075 

0.000020 

3.70 

0.0035 

14  W^*2 

0.035083 

0.014657 

2.39 

0.0356 

NO.  cases  =  25 

R-sq. 
R-sq-ad] . 

=  0.9627 

RMS  Error  = 

^  4.099 

Resid.  df  =  11 

=  0.9186 

Cond.  No.  = 

:  151.2 
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Least  Squares  Sunmary  ANOVA,  Response  SD  Model  PL0T2004 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  Total  (Ctorr.) 

24 

4953.242 

2 

Regression 

13 

4768.403 

366.800 

21.83 

0.0000 

3 

Linear 

5 

2193.701 

438.740 

26.11 

0.0000 

4 

NOn-linear 

8 

2574.703 

321.838 

19.15 

0.0000 

5 

Residual 

11 

184.838 

16.803 

6 

Lack  of  fit 

10 

184.838 

18.484 

7 

Pure  error 

1 

0.000 

0.000 

R-sq.  =  0.9627 

R-sq-ad] .  =  0 , 9186 

Model  obeys  hieraxx±y.  The  sum  of  squares  for  linear  terms 
is  cxxiputed  assuming  nonlinear  terms  are  first  removed. 
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Accepted  model  for  response  M_D 
0  1 


1  Model  Name: 

2  Response  Transformation: 

3  Method: 

4  weights: 

5  Total  Nunoiber  of  Cases: 

6  Number  of  Predictors: 

7  Number  of  Une>^>anded  Ter 

8  Number  of  Excluded  Ceuses 

9  Error  Degrees  of  Freedom 

10  Standard  Error  of  Respon 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


PIJ0T3502 

Untransfoxmed 

Least  Squares 

None 

31 

5 

15 
0 

16 

10.193068 

0.889395 

1.967447 
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Least  Squares  Coefficients,  Response  M,  Model  PLOT3502 


0  Term 

1  Coeff. 

2  Std.  Error 

3  T- value  4  Signif. 

1  1 

-24.193203 

7.342000 

2  Z 

-0.160725 

0.096374 

3  X 

-0.193163 

0.088468 

4  W 

-2.285052 

0.348288 

5  W 

0.006236 

0.001121 

6  W  H 

0.651004 

0.051447 

7  Z*X 

0.003661 

0.000760 

4.82 

0.0002 

8  Z*W 

0.011528 

0.001862 

6.19 

0.0001 

9  Z*W  H 

-0.002129 

0.000668 

-3.19 

0.0057 

10  X*W 

0.014530 

0.001632 

8.90 

0.0001 

11  X*W  H 

-0.002875 

0.000590 

-4.87 

0.0002 

12  W*W 

0.000091 

0.000023 

4.03 

0.0010 

13  W*W  H 

-0.003531 

0.001506 

-2.34 

0.0323 

14  W  H 

-0.000062 

0.000008 

-7.57 

0.0001 

15  W^*2 

0.022051 

0.005694 

3.87 

0.0013 

No.  cases 

=  31 

R-sg. 

=  0.9801 

RMS  Error 

=  1.967 

Resid.  df 

=  16 

R-sq-adJ. 

=  0.9627 

Cond.  No. 

=  116.9 

tiaiC)VA22  7R  X  5C 
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Least  Squares  Sumnary  ANOVA,  Response  M  Mcxiel  FL0T3502 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  Toted  (Corr.) 

30 

3116.959 

2 

Regression 

14 

3055.026 

218.216 

56.37 

0.0000 

3 

Linear 

5 

1926.404 

385.281 

99.53 

0.0000 

4 

Non-linear 

9 

1128.621 

125.402 

32.40 

0.0000 

5 

Residual 

16 

61.934 

3.871 

6 

Lack  of  fit 

15 

61.934 

4.129 

7 

Pure  error 

1 

0.000 

0.000 

R-sg.  s  0.9801 
R-sg-ad^.  =  0.9627 

Model  obeys  hierazohy.  The  sum  of  squares  for  linear  terms 
is  cottputed  assisning  nonlinear  terms  are  first  removed. 
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Acc^ted  model  for  response  SDIR 
0  1 


1  Model  Name: 

2  Rerocnse  'n^nsformation: 

3  Method: 

4  Weights: 

5  Total  Nunter  of  Cases: 

6  Number  of  Predictors: 

7  Nundser  of  Uhexpanded  Ter 

8  Nuiriber  of  Excluded  Cases 

9  Error  Degrees  of  Freedom 

10  standard  Error  of  Respcn 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


PIi0rr3503 

untransfomned 

Least  Squares 

None 

31 

5 

10 

0 

21 

35.056863 

0.694524 

14.629961 
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Least  Squares  Coefficients,  Response  S,  Model  PL0T3503 
0  Term  1  Coeff.  2  Std.  Error  3  T-value  4  Signif. 


1  1 

92.986651 

2  Z 

0.047889 

3  X 

-0.125998 

4  W 

5.550065 

5  W 

-0.029305 

6  W  H 

-0.993305 

7  Z*X 

0.029550 

8  Z*W 

-0.066073 

9  X*W 

-0.063864 

10  W_*W_H 

0.000221 

No.  cases  =  31 
Resid.  df  =  21 

R-sg 

R-sq-eUiJ 

42.109917 

0.331640 

0.302151 

0.631001 

0.007962 

0.304785 


0.005391 

5.48 

0.0001 

0.012802 

-5.16 

0.0001 

0.011397 

-5.60 

0.0001 

0.000056 

3.92 

0.0008 

0.8781 

rae  Error 

=  14.63 

0.8258 

Cond.  NO. 

=  87.33 

141 
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Least  Squares  Sumnary  ANDVA,  Response  S  Model  PLGrr3503 


0 

Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

1  Total  (Corr.) 

30 

36869.51 

2 

Regression 

9 

32374.76 

3597.20 

16.81 

0.0000 

3 

Linear 

5 

16442.47 

3288.49 

15.36 

0.0000 

4 

Non-linseu: 

4 

15932.29 

3983.07 

18.61 

0.0000 

5 

Residual 

21 

4494.75 

214.04 

6 

Lack  of  fit 

20 

4494.75 

224.74 

7 

Pure  error 

1 

0.00 

0.00 

R-sq.  B  0.8781 
R-sq-adf.  ■  0.8258 

Model  obeys  hierarchy,  sum  of  squares  for  linear  terms 
is  conputed  assuming  nonlinear  terms  are  first  removed. 


0 


Accepted  ncdel  for  response  SDIAG 


1 


1  ftodel  NSme: 

2  Response  Transfonnatlon: 

3  Method: 

4  Heights: 

5  Total  Nundber  of  Cases: 

6  Number  of  Predictors: 

7  Nuniber  of  UnescpaiKied  Ter 

8  Nuirber  of  Excl\ided  Cases 

9  Error  Degrees  of  Freedon 

10  Standard  Error  of  Respcn 

11  Relative  PRESS: 

12  Root  Mean  Squared  Error: 


PLOT3504 
Untransfoimed 
Least  Squares 
None 
31 
5 
11 
0 

20 

32.609024 

0.705662 

13.499653 
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Leeist  Squares  Coefficients,  Response  SD,  Model  PL0T3504 


0  Term 

1  Coeff. 

2  Std.  Error 

3  T-valxie 

4  Signif. 

1  1 

436.520515 

120,302831 

2  Z 

0.079125 

0.306166 

3  X 

-0.106779 

0.280290 

4  W 

5.149210 

0.587511 

5  W 

-0.028980 

0.007374 

6  W  H 

-6.187558 

1.745447 

7  Z*X 

0.028887 

0.004975 

5.81 

0.0001 

8  Z*W 

-0.065194 

0.011828 

-5.51 

0.0001 

9  X*W 

-0.060722 

0.010574 

-5.74 

0.0001 

10  W  *W  H 

0.000220 

0.000052 

4.21 

0.0004 

11  W_H**2 

0.018801 

0.006223 

3.02 

0.0067 

No.  Ceises 

=  31  R-SC[. 

=  20  R-sq-ad] . 

=  0.8857 

RMS  Error  »  13.5 

Resid.  df 

=  0.8286 

Oond.  No. 

x:  364.6 
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Least  Squares  Sumnary  ANDV;^,  Respcsise  SD  Model  PL0T3504 


Source 

1  df 

2  Sum  Sq. 

3  Mean  Sq. 

4  F-Ratio 

5  Signif. 

Toted  (Corr. ) 

30 

31900.45 

Regression 

10 

28255.64 

2825.56 

15.50 

0.0000 

Linear 

5 

12471.85 

2494.37 

13.69 

0.0000 

Ncn-linear 

5 

15783.79 

3156.76 

17.32 

0.0000 

Residual 

20 

3644.81 

182.24 

Lack  of  fit 

19 

3644.81 

191.83 

Pure  error 

1 

0.00 

0.00 

R-sq.  ■  0.8857 

R-sq-adJ.  «  0.8286 

Model  obeys  hierarchy,  nie  ston  of  squares  for  linear  terms 
is  cornputed  assiiming  nonlinear  terms  are  first  removed. 
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Monle  CMo  Sjmiiaiinn  F^Dfram  -  SIMIAC 

PROGRAM  MONTE2 

CC  lUS  PROGRAM  RUNS  A  MONTE  CARLO  SIMULAHON  OF  AN  ATTACK  OF  A 
CC  HARDEN  FAQUry  WITH  A  STICK  OF  CONVaniONAL  WEAPONS. 

CC 

CC  The  input  variables  for  each  attack  are  as  folbws: 

CC  MFC  -  mean  concrete  strength  fc' 

CC  9CFC-  standard  deviatkm  the  concrete  strength 
CC  MXL  -  mean  wall  height  or  beam  element  length 
CC  SIGXL  -  standard  deviation  of  the  wall  heig^ 

CC  MD  -  mean  depth  of  wall  or  thidmeas 

CC  SIGD  •  standard  deviation  of  the  wall  tfaickneas 

CC  MXBl  •  mean  X  coordinate  of  bomb  one 

CC  MYBl  •  mean  Y  coordinate  (rf*  bomb  one 

CC  MXB2- mean  X  coordinate  of  bomb  two 

OC  MYB2  -  mean  Y  coordinatB  of  bomb  two 

CC  MXB3  -  mean  X  coordinate  of  bomb  three 

CC  MYB3- mean  Y  coordinate  of  bomb  three 

CC  MXB4- mean  X  coordinate  of  bomb  four 

CC  MYB4  •  mean  Y  coordinate  of  bomb  four 

CC  SIGBC^IB  -  standard  deviation  for  each  bomb  based  on  standard 

ballistic  for 

CC  that  type  bomb 

CC  SIGRAb^  -  standard  deviation  of  foe  aimpoint  in  the  range,  or  X 
direction 

CC  given  the  aim  point  is  at  0,0 

CC  SIGDEFLECT  -  standard  deviation  of  foe  aimpoint  in  the  deflection, 
or  Y 

CC  direction  given  foe  aim  point  is  at  0,0 
CC  FC  -  cunent  value  of  concrete  strei^ 

CC  D  -  cunent  value  for  the  depth  or  foidoiess  of  foe  wall 

CC  XL  -  cunent  value  for  the  height  of  the  wall 

CC  RANGEl  •  cunent  value  of  the  rai^  for  the  hh  bond)  after  samplii^ 

for 

CC  foe  ballistic  error 

CC  RANGE  -  cunent  value  for  the  range  for  the  hh  bomb  after  sampling 
for 

CC  aim  point  error 

CC  KILL  -  flag  variable,  if  KILL  =  1  stmctioe  is  killed,  if  KILL  =  0 
structure 

CC  is  not  killed 

CC  ZVAL^XVAL  -  tranformed  coordinates  for  deflection  and  shear  analysis 
CCC 

REAL  MFC>!XUMD,SIGFC,SIGXL,SIGD,XB1,YB1,XB2,YB2,XB3,YB3, 

& 

SIGBJlANGEJ3EFlECTJlANGEIJ3EFLECT13ICa(ANGE,SIGDEFLECr,XVAL, 

&  CHECK1,SCHECK1,SCHECK2,SCHECK34>ENETRATEJXAGSHEAR, 

&  DEFLBCnON,WVAL3J),XB4,YB4,ZVAUDFT4XRECrSHEAR 


cc 

INIHjBI  RUNSLKHTTJCBREACHKimKTJCSIOIUC^^ 

OC 

CHARACTER  NUM*3 
KBREAOH) 

KHTD-O 

KDEFLECM) 

KSHDOH) 

KSHDIAO-O 

B»12 

C 

CC 

990  FORMATa5.17F72) 

OPEN  (UNIT-12JFILB-^fOKIHN) 

READ(12,'(A)r)NUM 

READ(12.*)  IRlJNS>4FC3IGFC,MXi:^GXU^^Ca3;XBl,YBl.XB^ 

&  ^,YB3^,YB4^(B,SIGRANC£^GDEFLECr,WVAL 
CC 

OPEN  (IJMT“134TL&^M0NT//NUM) 

OPEN  (lMT=144TLEF^SUMNf//NUM) 

RUN=0 

lOLL^O 

WRITE(13,*)  WVAL 
DO  100 1=1  JRUNS 
CC 

CC  Each  nnvnU  evaluate  the  idUprobUityfiv  Ok  fiwwofxns  in 
the  stick. 

CC  Sampling  for  tfie  concrete  stfo^  (PCX  wall  thickness  (D)  and  wall 
height 

CC  are  accomplished  initially  and  are  valid  lor  an  entire  slide 
sampling.  That  is, 

CC  FQD  and  XL  are  held  constant  for  the  analysis  of  four  different 

weapon 

CC  locations. 

CC 

CC  Sanity  checks  are  placed  on  the  values  of  FCJ3,XL  Also,  the 
response 

CC  surfeces  are  only  valid  over  a  specific  range. 

CC 

FC  =  RANNORM(MFC,SIGFQ 
1F(FC.LT.2500)FC=2500 
XL  =  RANN0RM(MXUS1GXL) 

IF  (XL  .GT.  168)  XL=168 
IF  (XL  .LT.  96)  XL=96 
D  =  RANN0RM(MD,S1GD) 

IF(D.GT.36)D=36 
IF(piT.  12)I>=12 

CC  The  individual  weapon  locations  are  detennined  by  first  samplii^  on  the 
CC  distribution  of  each  weapon  based  on  its  ballistic  error  prop^es. 

CC  This  returns  a  location  (RANC£L  DEFLECT! )  relative  to  the  stick  center  for 
CC  that  bomb.  Subsequently,  the  aim  point  error  is  applied  wMi  the 


OC  tim  poilt  steidani  devutiom  (sigrngB  and  a^deflect)  wilfa  RANGEIJSEFLSCTI 
C£  as  die  mean  kxatkn  ooonfinates.  Ite  RANCSyXFIJBCnON  coodiaale  of  the 
CX  cunent  weapon  is  then  establised.  This  8tq>  is  accomplished  four 
times 

(X  withinonenm. 

CC 

D025  J»1.4 

DTHEN 

RANGE1»RANN0RM(XB1^GB) 

RANCB-RANNORM(RANC£l^GRANGE) 

DEaECTl>RANN0RM(YB13ICB) 

DEHJECP4(ANN0RM(DEF^^ 

ELSEIF(J£Q.2)‘niEN 

RANGEl-itANNOi<M(XB23ICiB) 

RA^IGB-RAN^iORM(IL^NaEl.^^ 

DEFLECTl=«ANNORM(yB23IGB) 

DEF1£CT=41ANN0RM(DEF1H^ 

ELSEIF(J£Q.  3)1^4 

RANGE1=RANN0RM(XB331CB) 

R>^NG&«Ara40RM(RANGE13IC]RA^^ 

DEF1£CT1»RANN0RM(YB3^CB) 

DEFl£CT=dlANr«RM(D^^ 

ELSE 

RANGEl»RANNORM(XB4^GB) 

RANG&4UNNORM(RA^^13IGR^ 

DEFLBCT1*«ANN0RM(YB4.SIG8) 

DEIU1CT=RANN0RM(DE^ 

ENDIF 

CC 

5  WRITE(I3  *)  RUN,KIli4UhK3EjmjBCTJCJ)^ 

CC 

CC  The  RANGE,DEFLfiCT  is  initially  dieck  to  see  if  it  is  outside  the 
CC  effective  range  of  the  weaponjiere  assumed  to  be  80  feet 
CC  The  analysis  proceeds  with  RANGSJ[£FLECT  coonhaitBS  beii^ 
screened  to 

CC  see  whether  they  hit  the  stnicture  or  fell  within  the  breadi  zone. 

CC 

IF  ((ABS(RANGE)  .GT.  80)  .OR.  (ABS(DEFLECI)  .GT.  80))  GO  TO  15 
CC 

IF  ((ABS(RANGE)  .LT.  30)  J^ND. 

&  (ABS(DEFLECT)  .LT.  30))  THEN 
PENETRATE  =  (0.33)*(WVAL**(333))  -  (D*(0333)) 

IF  (PENETRATE  .GE  0)  THEN 
KHIT=KHIT+1 
GOTO  10 
ENDIF 
GO  TO  15 
ENDIF 
CC 

DFT-D/12 

XA=0.18 


AB^)378 

RB-XA*((WVAL*%667)yDFI)  -  AB*(WVAL**(  333)) 

CHECKRB»RB+30 

WRni(I3,*)  RANGEJ3EF1£CTJ<B,CHECKRB 

IF  ((ABS(DEFLBCT)  .GE  30)  .AND.  (ABS(DEFLBCD  IE.  affiCKRB).AND. 

&  (ABS(RANGE)XE.CHBCKRB))'niEN 

KBREAOMCBREACHfl 

CG  fOlO 

ENIXF 

IF  ((ABS(RANGE)  .GE  30)  j\ND.  (ABS(RANGE)  XE  CHBCKRB)  .AND. 

&  (ABS(DEFl£CT)i£.  CHBCKRB))  11104 
KBREACH^^CBREACHfl 
GOTO  10 
ENIXF 
CC 
OC 

CC  If  a  weapon  does  not  hit  the  structiae  or  fidi  within  the  breach  zone,  its 
CC  RANGEDEFLBCT  coonlinates  are  transfonned  into  local  coordinates  ZVAEXVAL 
CC  in  order  to  perfonn  deflectkn  and  diear  analysis.  This  procedure  deienranes 
CC  the  closest  point  on  the  stnictue  wfaicfa  would  see  the  maxinun  loading  fiom 
OC  thatweapoa  Ihe  comets  are  conadered  unsucceptiMe  to  deflection  ordiear 
CC  fiuluies  as  diey  are  caloilated.  As  diown  in  Ogure  4.1,  should  a  weapon  ftll 
CC  outside  a  distance  30  •  2.S*D  fiom  die  coderiine  of  the  stnictuie,  the 
CC  ZVAEXVAL  coordinates  are  based  on  the  distance  fiom  the  detonation  to  the 
CC  nearest  fiice  at  a  ^int  23  times  the  wall  dikdcness  (D).  This  estimate  is 
CC  ocnserra^andisfaaaedindiBftctdHtdiBOcnienarencxmBlly 
CC  hi^ily  reinforced  with  steel  extending  a  minimum  of  a  wall  thidmess  in  to  the 
CC  adjacent  wall.  The  subroutine  LOCATE  performs  the  determination  of  the 
CC  ZVAUXVAL  coordinates  for  each  weaport 
CC 

CALL  LOCA-re(D,RANGEDEFIJEiCT,ZVAEXVAL) 

CC 

WRrre(13,*)  RANGEDEFLBCT,ZVAUXVAL 
CC 

CC  Deflection,  diagonal  shear,  and  direct  diear  analysis  are  performed 
using  the 

CC  ZVAL,XVAL  detonation  locations  determined  above.  Sxxild  any  limit 
states 

CC  be  exceeded  the  procedure  jumps  ahead  and  sets  the  kill  flag  and 
goes  on  to  the 

CC  next  weapon  or  begins  a  new  run  after  recording  the  kill. 

CC 

IF  (WVAL  .LE  1000)  THEN 
DEFLECTION  =  DEF1EC(ZVAUXVAUD JQXL) 

ELSE 

DEFLECTION  =  DEF1EC2(ZVAL»XVAUDEC3CL) 

ENIXF 

CHECK1=XUD 

DCHECK*DEFLECnON/XL 

WRITE(13,*)  DEFLECnOHCHBCKlJOIBCK 

IF  ((CHECKl  .GE  5)  .AND.  (DCHECK  .GE  (0.10)))  THEN 


KDOmr-KDEFLECTrl 
GO  TO  10 
ENIXF 

IF  ((CHEOCl  XT.  5)  .AND.  (DCHBCK  .CE.  (0.06)))  THEN 
KI»FIXCr=KreFLBCT+l 
GO  TO  10 
ENMF 
CC 
CC 
CC 

CC  CHECK  DIAGONAL  SHEAR 
CC 

IF  (WVAL  XE.  500)  THEN 
DIAGSHEAR>SDiAG(ZVAXXVALAFC^) 

ELSE 

DIAGSHEAR=SDIAG2(ZVAL^ALAFC0CL) 

ENIXF 

CC  SCHECK1=^««MAL  DEPTH  TO  SPAN  RATIO  OR  XL/D  >5 
SCHBCKl=((n.5)*SQRTl(FQ*B*Dyi000 
CC  SCHBCK2-XL/D<-5ANDXL/D>=2 

SCHBCK2=((.667)*(10  +  (XL/D))*S(^T(FC)*B*D)/1000 
CC  SCHBCK3=  XL/D<2 

SCHEiCK3=(8*SQRT(FQ*B*Dyi000 

CC 

WRITEda,*)  DIAGSHEAR,SCHBCKl,SCHECK23CHEnG 
IF  ((CHBCKl  .GT.  S)  MfD.  (DIAGSHEAR  .GT.  SCHBCKl))  THQ4 
KSHDIAO=KSHDIAGi-l 
GOTO  10 
ENDIF 

IF  ((CHECKl  .LE.  5)  .AND.  (CHECi?  GE,  2)  .AND. 

&  (DIAGSHEAR  .GT.  SCHECK2))  iHEN 
KSHDIAG=KSHDIAG^-l 
GOTO  10 
ENDIF 

IF  ((CHBCKl  .LT.  2)  .AND.  (DIAGSHEAR  .GT.  SCHECK3))  THEN 
KSHDIAG  =  KSHDIAG^-1 
GOTO  10 
ENDIF 
CC 

CC  CHECK  DIRECT  SHEAR 
CC 

IF  (WVAL  .LE  500)  THEN 
DIRECTSHEAR=SDIR(ZVAXXVALX>,FC^) 

ELSE 

DIRECTSHEAR=SDIR2(ZVAXXVALX>XC^) 

ENDIF 

C 

SCHECK4=((0.51)*FC*B*Dyi000 
WRITE(13,*)  DIRECrSHEAR,SCHECK4 
IF  (DIRECTSHEAR  .GT.  SCHBCK4)  THEN 
KSHDIR=KSHDIR+1 


GOTO  10 
ENDIF 

CC 

GOTO  15 
CC 

CC  COLLECT  THE  NUMBER  OF  KILLS 
CC 

10  CONTINUE 
KILL  =  KILL+  1 
CC 

CC  After  all  nms  are  compkie  the  total  nunber  (Skills  is  divided  by 
the  total 

CC  numbo- of  runs  tunes  four  for  the  four  weapons  analyzed  per  nin,  to 
give  the 

CC  perc^probabilily  of  kill  for  a  single  attadc  with  the  input 
daracteristics. 

CC 

15  CONTINUE 
RUN  =  RUN  +  1 
CC 

25  CONTINUE 
CC 

100  CONTINUE 
K1U>KILL*100 
CC 

PK»  (KILL/RUN) 

WRrra(14,*)  Tirobability  of  Kill  (Pk)  % 

WRrTE(14,*) 

CC 

WRITER  14,*)  'KHIT  =  'JCHIT,’  KBREACH  =^4®REACH 
WRITE(14,*) 

WRITE<14,*)  'KDEFLBCT  =’,KDEFLECT,'  KSHDIAG  “'JCSHDIAG 
WRITE(14,*) 

WR1TE(14,*)'KSHDIR=  'JCSHDIR 
END 
CC 
CC 
CC 

FUNCTION  RANNORM(MU,SIG) 

REAL  MU,SIG 
CALL  SEEIXSEEDVAL) 

1[XJM=SEEDVAL 

RANNORM  =  SIG*GASDEV(1DUM)  +  MU 
END 
C 

FUNCTION  GASDEV(IDUM) 

DATA  ISET/(V 
IF  aSET.BQ.O)  THEN 
1  Vl=2*RANiaDUM>-l. 

V2=2.*RANiaDUM).l. 

R=V1**2+V2**2 


150 


IF(R.GE.1.X30T01 

FAOSQR'IT-2*LOG(RyR) 

GSET=V1*FAC 

GASDEV=V2*FAC 

ISET=1 

ELSE 

GASDEV=GSET 

ISET=0 

ENIXF 

RETURN 

END 

C 

FUNCTION  RANIODUM) 

DIMENSION  R(97) 

PARAMETER  (M1»2S9200M1=714UC1»S4773JIM1«3.SS80247E^ 
PARAMETER  (M2=134456M2=812UC2=2841 1 J»C=7.4373773&6) 
PARAMETER  (M3=243000M3=456UC3=51349) 

DATA  IFF /O/ 

IF  (IDUM1,T.O.ORJFFEQ.O)  THEN 
IFF=1 

DC1=M0D(IC1-IDUM^1) 

IX1=M0D(IA1  •IXl-HCl^ll) 

IX2=M0DaXl>C) 

DC1=M0D(IA1  *D(1+IC1  Ml) 

DO=MODaXlM3) 

DO  U 

D(1=M0D(IA1  *DC14IC1M1) 

IX2=M0D(IA2*DC2+IC2M2) 

R(J)=(ITjOATaXl  >fFlj0AT(IX2)*RM2)*RMl 
11  CONTTNUE 
IDUM=1 
ENDIF 

IX1=M0D(1A1*IX1+IC1M1) 

IX2=M0D(1A2*IX2+IC2M2) 

IX3=M0D(IA3*IX3+iaM3) 

J=l-K97*IX3yM3 
IF(J.GT.97.0RJ.LT.  1  )PAUSE 
RAN1=R(J) 

R(J)=<FLOA'r[lXl  >^-FLOAT(IX2)*RM2)•RMl 
RETURN 
END 
CC 
CC 
CC 

SUBROUTINE  LOCATE(DJUVNGEDEFLECr^AUXVAL) 
DCHECK  =  D*2.5/12 

IF  (ABS(RANGE)  .GT.  ABS(DEFLECI))  THEN 
IF  (  ABS(DEFLECT)  .LT.  (30  -  DCHBCK))  THEN 
ZVAL  =  ABS(RANGE)  -  30 
XVAL  =  0 
ELSE 


ZVAL  =  ABS(RANGE)-30 
XVAL  =  ABS(DEFLECT)  -  (30  -  DCHBCK) 
ENIXF 
ELSE 

IF  (ABS(RANCffi)  .LT.  (30  -  DCHBCK))  THEN 
ZVAL  =  ABS(DEFLECD  -  30 
XVAL  =  0 
ELSE 

ZVAL  =  ABSODEFLECT)  -  30 
XVAL  =  ABS(RANGE)  -  (30  -  DCHBCK) 
ENDIF 
ENDIF 
RETURN 
END 
CC 
CC 
CC 

FUNCTION  DEFLEC(ZAT4=CCy(LL) 

REAL  OOEF(l4)^AUXVAUFC;CUD 
DATACOEF  /-5.969702,  -0.187083, 


& 

-0.186465, 

-0.525068, 

& 

-0.001926, 

0.437493, 

& 

0.006853, 

0.011265, 

& 

-0.002346, 

0.011900, 

& 

-0.002847, 

0.000056, 

& 

-0.011774, 

0.019214/ 

C 

ZVAL  =  Z 
XVAL  =  X 
FC  =  FCC 
XL  =  XLL 
D=  T 

DEFLEC  =  COEF(l)  +  COEF(2)*ZVAL  +  COEF(3)*XVAL  +  COEF(4)*D 
&  +  COEF(5)*FC  +  COEF(6)*XL  +  COEF(7)*ZVAL*XVAL 

&  +  COEF(8)*ZVAL*D  +  COEF(9)*ZVAL*XL+  COEF(10)*XVAL*D 

&  +  COEF(l  1)*XVAL*XL  +  COEF(12)*D*FC  +  COEF(13)*D*XL 

&  +COEF(14)*(D**2) 

C 

PRINT*,  ZVAL,XVAUFC,XUD, DEFLEC 
RETURN 
END 
CC 
CC 
CC 

FUNCTION  SDIR(ZXTJCC,XLL) 

REAL  COEF(14),ZVAUXVAUFC,XUD 
DATA  COEF  /1 18.627133,  -0.640760, 

&  -0.560721,  4.347709, 

&  -0.020032,  -1.021689, 

&  0.014785,  -0.029421, 

&  0.008676,  -0.043227, 


&  0.000106,  -0.000246, 

&  0.000151,  -0.010774/ 

C 

ZVAL  =  Z 
XVAL  =  X 
FC  =  FCC 
XL  =  XLL 
D=  T 

SDIR  =  CX)EF(1)  +  CC>EF(2)*ZVAL  +  C()EF(3)*XVAL  +  CX)EF(4)*D 
&  +  CX)EF(5)*FC  +  COEF(6)*XL  +  CX)EF(7)*ZVAL*XVAL 

&  +  CX)EF(8)*ZVAL*D  +  CX)EF(9)*ZVAL*XL  +  COEF(10)*XVAL*D 

&  +  COEF(l  1)*XVAL*FC  +  COEF(12)*D*FC  -KX)EF(13)*FC*XL 

&  +  aDEF(14)*(ZVAL**2) 

C 

PRINT*,  ZVAL,XVAUFC,XLJ),SDIR 

RETURN 

END 


FUNCTION  SDIAG(Z,X,TJPCC,XLL) 
REAL  COEF(14),ZVAL,XVAUFC,XUD 
DATACOEF  /81.481744,  -0J11996, 


& 

-0.8%116, 

0.981792, 

& 

-0.009560, 

-0.557030, 

& 

0.007351, 

-0.025476, 

& 

0.000083, 

-0.009849, 

& 

0.005363, 

-0.000220, 

& 

0.000075, 

0.035083/ 

C 

ZVAL  =  Z 
XVAL  =  X 
FC  =  FCC 
XL  =  XLL 
D=  T 

SDIAG  =  COEF(l)  +  COEF(2)*ZVAL  +  COEF(3)*XVAL  +  COEF(4)*D 
&  +  COEF(5)*FC  +  OOEF(6)*XL  +  COEF(7)*ZVAL*XVAL 

&  +  COEF(8)*ZVAL*D  +  COEF(9rZVAL*FC  +  COEF(10)*XVAL*D 

&  +  COEF(  1 1  )*XVAL*XL  +  COEF(12)*D*FC  +  CQEF(13)*FC*XL 

&  +COEF(14)*(D**2) 

C 

PRINT*,  ZVAL,XVAUFC,XL,D,SDIAG 
RETURN 
END 
CC 

FUNCTION  DEFl.fiC2(Z,X,T4^CC,XLL) 

REAL  COEF(15),ZVAL,XVAL,FC,XL,D 
DATA  COEF  /-24.193203,  -0.160725, 


& 

-0.193163, 

-2285052, 

& 

0.006236, 

0.651004, 

& 

0.003661, 

0.011528, 

& 

-0.002129, 

0.0I4S30, 

& 

-0.002875, 

0.000091, 

& 

-0.003531, 

-0.000062, 

& 

0.022051/ 

c 

ZVAL  =  Z 
XVAL  =  X 
FC  =  FCC 
XL  =  XLL 
D=  T 

DEFLEC2  =  C0EF(1)  +  C0EF(2)*ZVAL  +  CX)EF(3)*XVAL  +  C0EF(4)*D 
&  +  CXDEF(5)*FC  +  C0EF(6)*XL  +  C0EF(7)*ZVAL*XVAL 

&  +  C0EF(8)*ZVAL*D  +  CX)EF(9)*ZVAL*XL  +  COEF(10)*XVAL*D 

&  +CX)EF(11)*XVAL*XL  +  CX}EF(12)*D*FC  +  CX)EF(13)*D*XL 

&  +CQEF(14)*FC*XL-KX3EF(15)*(D**2) 

C 

PRINT*,  ZVAUCVALJC,XUDJ)EFLEC2 
RETURN 
END 
<X 

FUNCTION  SDIR2(Z,X,T,FOC,XLL) 

REAL  COEF(10),ZVAL,XVAUFC,XLJ> 

DATACOEF  /92.986651,  0.047889, 


8l 

-0.125998, 

5.550065, 

8l 

-0.029305, 

-0.993305, 

& 

0.029550, 

-0.066073, 

& 

-0.063864, 

0.000221/ 

C 

ZVAL  =  Z 
XVAL  =  X 
FC  =  FCC 
XL  =  XLL 
D=  T 

SD1R2  =  COEF(l)  +  COEF(2)*ZVAL  +  COEF(3)*XVAL  +  CQEF(4)*D 
&  +  COEF(5)*FC  +  COEF(6)*XL  +  COEF(7)*ZVAL*XVAL 

&  +  COEF(8)*ZVAL*D  +  COEF(9)*XV.U*D  +  CQEF(10)*FC*XL 

C 

PRINT*,  ZVAUXVAUFC,XUD,SDIR2 
RETURN 
END 
CC 

FUNCTION  SDIAG2(Z.XT,FCC,XLL) 

REAL  COEF(  1 1  ),ZVAL,XVAL,FC,XL,D 
DATACOEF  /436.520515,  0.079125, 


& 

-0.106779, 

5.149210, 

& 

-0.028980, 

-6.187558, 

& 

0.028887, 

-0.065194, 

& 

-0.060722, 

0.000220, 

& 

0.018801/ 

C 

ZVAL  =  Z 


XVAL-X 
FC«FOC 
XL  =  XLL 
D=  T 

SDIAG2  =  C0EF(1)  +  CX)EF(2)*ZVAL  +  C0EF(3)*XVAL  +  C0EF(4)*D 
&  +CX)EF(5)*FC  +  COEF(6)*XL+  CQEF(7)*ZVAL*XVAL 

&  +(XIEF(8)*ZVAL*D  +  CX)EF(9)*XVAL*D  +  COEF(10)*FC*XL 

&  +C0EF(11)*(XL**2) 


PRINT*,  ZVAIJ(VAU5C,XU^SDIAG2 
RETURN 


